Electronic transport in two dimensional graphene 



S. Das Sarma^, Shaffique Adam^^^, E. H. Hwang^, and Enrico Rossi^ 

^ Condensed Matter Theory Center, 
Department of Physics, 
University of Maryland, College Park, 
MD 20742-4111 

'^Center for Nanoscale Science and Technology, 
National Institute of Standards and Technology, 
Gaithersburg, Maryland 20899-6202, 
USA 

(Dated: November 9, 2010) 

We provide a broad review of fundamental electronic properties of two-dimensional 
graphene with the emphasis on density and temperature dependent carrier transport 
in doped or gated graphene structures. A salient feature of our review is a critical com- 
parison between carrier transport in graphene and in two-dimensional semiconductor 
systems (e.g. heterostructures, quantum wells, inversion layers) so that the unique fea- 
tures of graphene electronic properties arising from its gapless, massless, chiral Dirac 
spectrum are highlighted. Experiment and theory as well as quantum and semi-classical 
transport are discussed in a synergistic manner in order to provide a unified and com- 
prehensive perspective. Although the emphasis of the review is on those aspects of 
graphene transport where reasonable consensus exists in the literature, open questions 
are discussed as well. Various physical mechanisms controlling transport are described 
in depth including long-range charged impurity scattering, screening, short-range defect 
scattering, phonon scattering, many-body effects, Klein tunneling, minimum conductiv- 
ity at the Dirac point, electron-hole puddle formation, p-n junctions, localization, per- 
colation, quantum-classical crossover, midgap states, quantum Hall effects, and other 
phenomena. 
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I. INTRODUCTION 
A. Scope 

The experimental discovery of two dimensional (2D) 

gated graphene in 2004 by Novoselov et al. (20041 is a 



seminal event in electronic materials science, ushering in 
a tremendous outburst of scientific activity in the study 
of electronic properties of graphene which continues un- 
abated upto the end of 2009 (with the appearance of more 
than 5000 articles on graphene during the 2005-2009 five- 
year period). The subject has now reached a level so vast 
that no single article can cover the whole topic in any rea- 
sonable manner, and most general reviews are likely to 
become obsolete in a short time due to rapid advances in 
the graphene literature. 

The scope of the current review (written in late 2009 
and early 2010) is transport in gated graphene with the 
emphasis on fundamental physics and conceptual issues. 
Device applications and related topics are not discussed 



(Avouris et al. 




2007 


) 1 


properties ( 


Bunch et al. 



al. 2008a). The 



important subject of graphene materials science, which 
deserves its own separate review, is not discussed at all. 
Details of the band structure properties and related phe- 
nomena are also not covered in any depth, except in the 
context of understanding transport phenomena. What is 
covered in reasonable depth is the basic physics of car- 
rier transport in graphene, critically compared with the 
corresponding well-studied 2D semiconductor transport 
properties, with the emphasis on scattering mechanisms 
and conceptual issues of fundamental importance. In the 
context of 2D transport, it is conceptually useful to com- 
pare and contrast graphene with the much older and well 
established subject of carrier transport in 2D semicon- 
ductor structures (e.g. Si inversion layers in MOSFETs, 
2D GaAs heterostructures and quantum wells). Trans- 
port in 2D semiconductor systems has a number of sim- 
ilarities and key dissimilarities with graphene. One pur- 
pose of this review is to emphasize the key conceptual dif- 
ferences between 2D graphene and 2D semiconductors so 
as to bring out the new fundamental aspects of graphene 
transport which make it a truly novel electronic material 
qualitatively different from the large class of existing and 
well established 2D semiconductor materials. 

Since graphene is a dynamically (and exponentially) 
evolving subject, with new important results appearing 
almost every week, the current review concentrates on 
only those features of graphene carrier transport where 



some qualitative understanding, if not a universal consen- 
sus, has been achieved in the community. As such, some 
active topics, where the subject is in flux, have been left 
out. Given the constraint of the size of this review, depth 
and comprehension have been emphasized over breadth - 
given the huge graphene literature, no single review can 
attempt to provide a broad coverage of the subject at 
this stage. There have already been several reviews of 
graphene physics in the recent literature. We have made 
every effort to minimize overlap between our article and 
these recent reviews. The closest in spirit to our review 



is the one by Castro Neto et al. ( 2009 ) which was written 



2.5 years ago (i.e. more than 3000 graphene publications 
have appeared in the literature since the writing of that 
review) . Our review should be considered complimentary 
to Castro Neto et al. (2009), and we have tried avoiding 



too much repetition of the materials already covered by 
them, concentrating instead on the new results arising 
in the literature following the older review. Although 
some repetition is necessary in order to make our review 
self-contained, we refer the reader to [Castro Neto et al.\ 
(2009) for details on the history of graphene, its band 



structure considerations and the early (2005-2007) ex- 
perimental and theoretical results. Our material empha- 
sizes the more mature phase (2007-2009) of 2D graphene 
physics. 

For further background and review of graphene physics 
beyond the scope of our review, we mention in addition to 



the Rev. Mod. Phys. article by [Castro Neto et aZ.| ( |2009l ), 
the accessible reviews by Andrei Geim and his collabo- 
rators (Geim 2009 Geim and Novoselov 2007), the re- 



cent brief review by (Mucciolo and Lewenkopf 20101, 



as well as two edited volumes of Solid State Communi- 



cations (Das Sarma et al. 2007a Fal'ko et al. 20091, 



where the active graphene researchers have contributed 
individual perspectives. 



B. Background 

Graphene (or more precisely, monolayer graphene - in 
this review, we refer to monolayer graphene simply as 
"graphene") is a single 2D sheet of carbon atoms in a 
honeycomb lattice. As such 2D graphene rolled up in 
the plane is a carbon nanotube, and multilayer graphene 
with weak interlayer tunneling is graphite. Given that 
graphene is simply a single 2D layer of carbon atoms 
peeled off a graphite sample, there was early interest 
in the theory of graphene band structure which was all 
worked out a long time ago. In this review we only con- 
sider graphene monolayers (MLG) and bilayers (BLG), 
which are both of great interest. 
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FIG. 1 (Color online) (a) Graphene honeycomb lattice show- 
ing in different colors the two triangular sublattices. Also 
shown is the graphene Brillouin zone in momentum space. 

(b) 



Adapted from 
otube as a rol 



Castro Neto et al. 



(2009|). (b) Carbon nan- 
(c) Lattice structure 



ed up graphene layer, 
of graphite, graphene multilayer. Adapted from Castro Neto| 
\et al.\ ( [2006| ). (d) Lattice structure of bilayer graphene. 70 
and 71 are respectively the intralayer and interlayer hopping 
parameters t, t±_ used in the text. The interlayer hopping pa- 
rameters 73 and 74 are much smaller than 71 = t±_ and are 
normally neglected. Adapted from |Mucha-Kruczynski et al\ 
( |2010| ) (e) Typical configuration for gated graphene. 
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FIG. 2 ( Color on l ine) ( a) Graphene bandstructure. Ad- 
pated from Wilson] ( 2006 1 . (b) Enlargment of the band struc- 



ture close to the K and K' points showing the Dirac cones. 
Adpated from |Wilson| (|2006|. (c) Model energy dispersion 
E — ftup|k|. (d) Density of states of graphene close to the 
Dirac point. The inset shows the density of states over the 
full electron bandwidth. Adapted from [Castro Neto et al.\ 



(20091 



1. Monolayer graphene 



Graphene monolayers have been rightfully described 



as the "ultimate flatland" ( Geim and MacDonald 2007 ) 



i.e. the most perfect 2D electronic material possible in 
nature, since the system is exactly one atomic monolayer 
thick, and carrier dynamics is necessarily confined in this 
strict 2D layer. The electron hopping in 2D graphene 
honeycomb lattice is quite special since there are two 



equivalent lattice sites (A and B in Fig. LB 1) which give 
rise to the "chirality" in the graphene carrier dynamics. 

The honeycomb structure can be thought of as a tri- 
angular lattice with a basis of two atoms per unit cell, 
with the 2D lattice vectors Aq = (a/2)(3, -y/S) and 
Bo = (a/2)(3,-%/3) (a w 0.142 nm is the carbon- 
carbon distance). K = (27r/(3a), 27r/(3>/3a)) and K' — 
(27r/(3a),-27r/(3V3a)) are the inequivalent corners of 
the Brillouin zone and are called Dirac points. These 
Dirac points are of great importance in the electronic 
transport of graphene, and play a role similar to the role 
of r points in direct band-gap semiconductors such as 
GaAs. Essentially all of the physics to be discussed in 
this review is the physics of graphene carriers (electrons 
and/or holes) close to the Dirac points (i.e. within a 2D 
wavevector g = |q| <C 27r/a of the Dirac points) just as 
all the 2D semiconductor physics we discuss will occur 
around the F point. 

The electronic band dispersion of 2D monolayer 



graphene was calculated by Wallace (1947) and oth- 



ers (McClure 1957 Slonczewski and Weiss 1958) a long 



time ago, within the tight-binding prescription, obtaining 
upto the second-nearest-neighbor hoping term in the cal- 
culation, the following approximate analytic formula for 
the conduction (upper, +, tt*) band, and valence (lower, 
— , tt) band. 



E±{q) ^Zt' ±hvY\q\ 



m'a^ 3ta^ . 
—— ± — sm(30,) 



(1.1) 

with vp = 3ta/2, dg — arctan ^[qx/qy], and where 
t, t' are respectively the nearest-neighbor (i.e. inter- 
sublattice, A-B) and next-nearest-neighbor (i.e intra- 
sublattice, A-A or B-B) hopping amplitudes and i(w 
2.5 eV) » f (« 0.1 eV). 

The almost universally used graphene band dispersion 
at long wavelength puts t' — 0, where the band-structure 
for small q relative to the Dirac point is given by 

E±{q)^±hvFq + 0{q/kf. (1.2) 

Further details on the band structure of 2D graphene 



monolayers can be found in the literature ( 


Castro Neto 


et al.\ 20091 McClure| 1957 Mcclure| 1964 


Reich et al. 


2002 Slonczewski and Weiss 1958 Wallace 1947) and 



will not be discussed here. Instead we provide below. 



a thorough discussion of the implications of Eq. 1.2 for 
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graphene carrier transport. Since much of the fundamen- 
tal interest is in understanding graphene transport in the 
relatively low carrier density regime, complications aris- 
ing from the large q (~ K) aspects of graphene band 
structure lead only to small perturbative corrections to 
graphene transport properties and are, as such, only of 
secondary importance to our review. 

The most important aspect of graphene's energy dis- 
persion (and the one attracting most attention) is its lin- 
ear energy-momentum relationship with the conduction 
and valence bands intersecting at 9 = 0, with no energy 
gap. Graphene is thus a zero band-gap semiconductor 
with a linear rather than quadratic long wavelength en- 
ergy dispersion for both electrons (holes) in the conduc- 
tion (valence) bands. The existence of two Dirac points 
at K and K' where the Dirac cones for electrons and 
holes touch (Fig. ^p) each other in momentum space, 
gives rise to a valley degeneracy = 2 for graphene. 
The presence of any inter-valley scattering between K 
and K' points lifts this valley degeneracy, but such ef- 
fects require the presence of strong lattice scale scatter- 
ing. Intervalley scattering seems to be weak and when 
they can be ignored the presence of a second valley can 
be taken into account symply via the degenercy factor 
— 2. Throughout this introduction we neglect inter- 
vally scattering processes. 

The graphene carrier dispersion E±{q) — hv^q, ex- 
plicitly depends on the constant v-p, sometimes called 
the graphene (Fermi) velocity. In the literature differ- 
ent symbols {vp, vq, j/h) are used to denote this velocity. 
The tight-binding prescription provides a formula for vp 
in terms of the nearest-neighbor hopping t and the lattice 
constant 02 — \/3a: hvp — 3ta/2. The best estimates of 
t ^ 2.5 eV and a = 0.14 nm give vp w 10^ cm/s for the 
empty graphene band i.e. in the absence of any carriers. 
Presence of carriers may lead to a many-body renormal- 
ization of the graphene velocity, which is, however, small 
for MLG, but could, in principle, be substantial for BLG. 

The linear long wavelength Dirac dispersion with 
Fermi velocity roughly 1/300 of the corresponding veloc- 
ity of light, is the most distinguishing feature of graphene 
in addition to its strict 2D nature. It is therefore natural 
to ask about the precise applicability of the linear en- 
ergy dispersion, since it is obviously a long-wavelength 
continuum property of graphene carriers valid only for 
q<.Kfv (0.1 nm)-i. 

There are several ways to estimate the cut-off wavevec- 
tor (or momentum) kc above which the linear contin- 
uum Dirac dispersion approximation breaks down for 
graphene. The easiest is perhaps to estimate the carrier 
energy Ec = hvpkc, and demand that Ec < QAt 1.0 eV, 
so that one can ignore the lattice effects (which lead to 
deviations from pure Dirac-like dispersion). This leads 
to a cutoff wavevector given by kc « 0.25 nm~^. 

The mapping of graphene electronic structure onto the 
massless Dirac theory is deeper than the linear graphene 



carrier energy dispersion. The existence of two equiva- 
lent, but independent, sub-lattices A and B (correspond- 
ing to the two atoms per unit cell), leads to the exis- 
tence of a novel chirality in graphene dynamics where 
the two linear branches of graphene energy dispersion (in- 
tersecting at Dirac points) become independent of each 
other, indicating the existence of a pseudospin quantum 
number analogous to electron spin (but completely in- 
dependent of real spin). Thus graphene carriers have 
a pseudospin index in addition to the spin and orbital 
index. The existence of the chiral pseudospin quantum 
number is a natural byproduct of the basic lattice struc- 
ture of graphene comprising two independent sublattices. 
The long-wavelength low energy effective 2D continuum 
Schrodinger equation for spinless graphene carriers near 
the Dirac point therefore becomes 



ihvpcr ■ V^(r) = £'*(r) 



(1.3) 



where cr — {ax,<yy) is the usual vector of Pauli matri- 
ces (in 2D now), and ^'(r) is a 2D spinor wavefunction. 
Eq.( |1.3[ ) corresponds to the effective low energy Dirac 
Hamiltonian: 



= hvp 







iqy 



iqy 







hvpcr ■ q. (1.4) 



We note that Eq. |1.3| is simply the equation for mass- 
less chiral Dirac Fermions in 2D (except that the spinor 
here refers to the graphene pseudospin rather than real 
spin), although it is arrived at starting purely from the 
tight-binding Schrodinger equation for carbon in a honey- 
comb lattice with two atoms per unit cell. This mapping 
of the low energy, long wavelength electronic structure 
of graphene onto the massless chiral Dirac equation, was 



discussed by Semenoff ( 1984 ) more than 25 years ago 



It is a curious historical fact that although the actual 
experimental discovery of gated graphene (and the be- 
ginning of the frenzy of activities leading to this review 
article) happened only in 2004, some of the key theo- 
retical insights go back a long way in time and are as 
valid today for real graphene as they were for theoreti- 



cal graphene when they were introduced ( 


Gonzalez et al. 


19941 Haldanel 


19881 Ludwig et oLl|1994| 


McClure 


1957 


Semenoff 1984 


Wallace 19471. 



The momentum space pseudo-spinor eigenfunctions for 



Eq. 1.3 can be written as 



*(q,if) 



1 / g-»0,/2 
71 V ±6^'^'/^ 



*(q,^') 



1 / e'^"/'^ 
71 ±e-«^/2 



where the ± signs correspond to the conduction/valence 
bands with E±(q) = ±hvpq. It is easy to show using 
the Dirac equation analogy that the conduction/ valence 
bands come with positive/negative chirality, which is con- 
served, within the constraints of the validity of Eq. |1.3| 
We note that the presence of real spin, ignored so far. 
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would add an extra spinor structure to graphene's wave- 
function (this real spin part of the graphene wavefunc- 
tion is similar to that of ordinary 2D semiconductors). 
The origin of the massless Dirac description of graphene 
lies in the intrinsic coupling of its orbital motion to the 
pseudospin degree of freedom due to the presence of A 
and B sublattices in the underlying quantum mechanical 
description. 



2. Bilayer graphene 

The case of bilayer graphene is interesting in its own 
right, since with two graphene monolayers that are 
weakly coupled by interlayer carbon hopping, it is inter- 
mediate between graphene monolayers and bulk graphite. 

The tight-binding description can be adapted to study 
the bilayer electronic structure assuming specific stack- 
ing of the two layers with respect to each other (which 
controls the interlayer hopping terms). Considering the 
so-called A-B stacking of the two layers (which is the 
3D graphitic stacking), the low energy, long wavelength 
electronic structure of bilayer graphene is described by 



the following energy dispersion relation (Brandt et al. 


1988 


Dresselhaus and Dresselhaus||2002||McCann||2006 


McCann and Fal'ko 


20061 



E±{q)=[V^ 



1/2 



1/2 

(1,5) 



where is the effective interlayer hopping energy (and 
t, v-p are the intra-layer hopping energy and graphene 
Fermi velocity for the monolayer case). We note that 
t_L(~ 0.4 eV) < t(~ 2.5 eV), and we have neglected sev- 
eral additional interlayer hopping terms since they are 
much smaller than t^^. The quantity V with dimensions 
of energy appearing in Eq. |1.5| for bilayer dispersion cor- 
responds to the possibility of a real shift (e.g. by an ap- 
plied external electric field perpendicular to the layers, 
z-direction) in the electrochemical potential between the 
two layers, which would translate into an effective band- 



gap opening near the Dirac point ( 


Castro et al. 


2007 


Ohta et al. 2006 Oostinga et al. 


2008 


Zhang et al. 



2009c I 



Expanding Eq. |1.5| to leading order in momentum, and 
assuming (1/ <C i), we get 

E±{q) = ±[V - 2h\lVq^/tl + h\U^/{2tiV)]. (1.6) 

We conclude that (i) For V ^ bilayer graphene has a 
minimum bandgap of A = 2V —W^ /t\ at g = y/2V/hvp, 
and (ii) for F = 0, bilayer graphene is a gapless semi- 
conductor with a parabolic dispersion relation E± [q) « 
t?v\q^ jt^ = h^q^ /{2m), where m — <j^/(2w|.) for small 
q. The parabolic dispersion (for V = Q) applies only for 
small values of q satisfying hv-pq <C t_L, whereas in the 



opposite limit hv-pq >• i_L, we get a linear band disper- 
sion E±{q) ±hvpq, just hke the monolayer case. We 
note that using the best estimated values for vp and t±, 
the bilayer effective mass is m « (0.03 to 0.05) mo, which 
corresponds to a very small effective mass. 

To better understand the quadratic to linear crossover 
in the effective BLG band dispersion, it is convenient 
to rewrite the BLG band dispersion (for y = 0) in the 
following hyperbolic form 



Eblg = Tmvl ± mu| [l -t- [k/kof] 



211/2 



(1.7) 



where fco = t±/{2hvp) is a characteristic wavevector. In 
this form it is easy to see that -Eblg — ^ k'^{k) for k — 
0(oo) for the effective BLG band dispersion with k k^ 
(k 3> fco) being the parabolic (linear) band dispersion 
regimes, fco ~ 0.3 nm~^ for m ^ 0.03 me. Using the best 
available estimates from band structure calculations, we 
conclude that for carrier densities smaller (larger) than 
5 X 10^^ cm~^, the BLG system should have parabolic 
(linear) dispersion at the Fermi level. 

What about chirality for bilayer graphene? Although 
the bilayer energy dispersion is non-Dirac like and 
parabolic, the system is still chiral due to the A/B sub- 
lattice symmetry giving rise to the conserved pseudospin 
quantum index. The detailed chiral 4-component wave- 
function for the bilayer case including both layer and 
sublattice degrees of freedom can be found in the origi- 



nal literature 


McCann 


2006 


McCann and Fal'ko 


2006 


Nilsson et al. 


2006a|b 


2008 





The possible existence of an external bias-induced 
band-gap and the parabolic dispersion at long wavelength 
distinguish bilayer graphene from monolayer graphene, 
with both possessing chiral carrier dynamics. We note 
that bilayer graphene should be considered a single 2D 
system, quite distinct from "double-layer" graphene. 



Hwang and Das Sarma ( 2009a I which is a composite sys- 



tem consisting of two parallel single layers of graphene, 
separated by a distance in the z-direction. The 2D energy 
dispersion in double-layer graphene is massless Dirac-like 
(as in the monolayer case), and the interlayer separation 
is arbitrary, whereas bilayer graphene has the quadratic 
band dispersion with a fixed inter-layer separation of 
0.3 nm similar to graphite. 



3. 2D Semiconductor structures 

Since one goal of this review is to understand graphene 
electronic properties in the context of extensively stud- 
ied (for more than 40 years) 2D semiconductor systems 
(e.g. Si inversion layers in MOSFETs, GaAs-AlGaAs 
heterostructures, quantum wells, etc.), we summarize in 
this section the basic electronic structure of 2D semi- 
conductor systems which are of relevance in the context 
of graphene physics, without giving much details, which 
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conductor structures, the quantum dynamics is two 
dimensional by virtue of confinement induced by an 
external electric field, and as such, 2D semiconductors 
are quasi-2D systems, and always have an average width 
or thickness (z) 5 nm to 50 nm) in the third direction 
with (z) < Af, where Ap is the 2D Fermi wavelength 
(or equivalently the carrier de Broglie wavelength). The 
condition {z) < Ap defines a 2D electron system. 



FIG. 3 (Color online) (a) Energy band of bilayer graphene 
for V = Q. (b) Enlargment of the energy band close to the 
neutrality point K for different values of V . Adapted from 
Min et al.\ ^mf\ 



can be found in the literature ( Ando ei 1982 Bastard 



1991 Davies 19981 



There are, broadly speaking, four qualitative dif- 
ferences between 2D graphene and 2D semiconductor 
systems. (We note that there are significant quantitative 
and some qualitative differences between different 2D 
semiconductor systems themselves). These differences 
are sufficiently important so as to be emphasized right 
at the outset. 

(i) First, 2D semiconductor systems typically have 
very large (> 1 eV) bandgaps so that 2D electrons and 
2D holes must be studied using completely different 
electron-doped or hole-doped structures. By contrast, 
graphene (except biased graphene bilayers that have 
small band-gaps) is a gapless semiconductor with the 
nature of the carrier system changing at the Dirac 
point from electrons to holes (or vice versa) in a single 
structure. A direct corollary of this gapless (or small 
gap) nature of graphene, is of course, the "always 
metallic" nature of 2D graphene, where the chemical 
potential ("Fermi level") is always in the conduction or 
the valence band. By contrast, the 2D semiconductor 
becomes insulating below a threshold voltage, as the 
Fermi level enters the bandgap. 

(ii) Graphene systems are chiral, whereas 2D semi- 
conductors are non-chiral. Chirality of graphene leads 
to some important consequences for transport behav- 
ior, as we discuss later in this review. (For example, 
2fcp-backscattering is suppressed in MLG at low tem- 
perature.) 

(Hi) Monolayer graphene dispersion is linear, whereas 
2D semiconductors have quadratic energy dispersion. 
This leads to substantial quantitative differences in the 
transport properties of the two systems. 

(iv) Finally, the carrier confinement in 2D graphene 
is ideally two-dimensional, since the graphene layer is 
precisely one atomic monolayer thick. For 2D semi- 



The carrier dispersion of 2D semiconductors is given by 
E{q) = Eo-\- h'^q'^ /{2m*), where Eq is the quantum con- 
finement energy of the lowest quantum confined 2D state 
and q = {qx,qy) is the 2D wavevector. If more than one 
quantum 2D level is occupied by carriers - usually called 
"sub-bands" - the system is no longer, strictly speaking, 
two dimensional, and therefore a 2D semiconductor is no 
longer two dimensional at high enough carrier density 
when higher subbands get populated. 

The 2D effective mass entering m* is known from band- 
structure calculations, and within the effective mass ap- 
proximation m* — 0.07 Too (electrons in GaAs), m* ~ 
0.19 Too (electrons in Si 100 inversion layers), m* = 
0.38 Too (holes in GaAs), and to* = 0.92 TOo (electrons 
in Si 111 inversion layers). In some situations, e.g. Si 
111, the 2D effective mass entering the dispersion rela- 
tion may have anisotropy in the xy plane and a suitably 
averaged to* = ^m^my is usually used. 

The 2D semiconductor wavefunction is non-chiral, and 
is derived from the effective mass approximation to be 



$(r,z) 



(1.8) 



where q and r are 2D wavevector and position, and ^(z) is 
the quantum confinement wavefunction in the z-direction 
for the lowest sub-band. The confinement wavefunc- 
tion defines the width/thickness of the 2D semiconductor 
state with (z) = KCk^lOI^^^- The detailed form for ^(z) 
usually requires a quantum mechanical self-consistent lo- 
cal density approximation calculation using the confine- 
ment potential, and we refer the reader to the extensive 
existing literature for the details on the confined quasi-2D 



subband structure calculations (Ando et al, 1982 Bas- 



tard 


1991 


Davies 


1998 



Stern and Das Sarma 1984) 



Finally, we note that 2D semiconductors may also in 
some situations carry an additional valley quantum num- 
ber similar to graphene. But the valley degeneracy in 
semiconductor structures e.g. Si MOSFET 2D electron 
systems, have nothing whatsoever to do with a pseu- 
dospin chiral index. For Si inversion layers, the valley de- 
generacy {gy =2,4 and 6 respectively for Si 100, 110 and 
111 surfaces) arises from the bulk indirect band struc- 
ture of Si which has 6 equivalent ellipsoidal conduction 
band minima along the 100, 110 and 111 directions about 
85 % to the Brillouin zone edge. The valley degener- 
acy in Si MOSFETs, which is invariably slightly lifted 
(w 0.1 meV), is a well established experimental fact. 
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TABLE I Elementary electronic quantities. Here Ep, D{E), r^, and qtf represent the Fermi energy, the density of states, the 
interaction parameter, and Thomas Fermi wave vector, respectively. Do = D{Ef) is the density of states at the Fermi energy 
and qs = qTr/kp- 





Ef 


D{E) 


Do = D{Ef) 




qTF 


Qs 


MLG 


hvFx/ — 




•Jasa^n 


<P Va^av 


^AiT g s av^^^ 


asavf^^ 






Khvp 2 


Klivp 


Khvp 


BLG/2DEG 




aaav'm 




me^ 3b g„ 










2nh^ 






Kh-^ \/ 47rri 




FIG. 4 (Color online) (a) Heterostructure inversion layer 
quantum well, (b) Band diagram showing the bending of 
the bands at the interface of the semiconductors and the 2 
dimensional subband. 



C. Elementary electronic properties 

We describe, summarize, and critically contrast the el- 
ementary electronic properties of graphene and 2D semi- 
conductor based electron gas systems based on their long 
wave length effective 2D energy dispersion discussed in 
the earlier sections. Except where the context is obvi- 
ous, we use the following abbreviations from now on: 
MLG (monolayer graphene), BLG (bilayer graphene), 
and 2DEG (semiconductor based 2D electron gas sys- 
tems). The valley degeneracy factors are typically = 2 
for graphene and Si 100 based 2DEGs, whereas ffv = 1 (6) 
for 2DEGs in GaAs (Si 111). The spin degeneracy is al- 
ways gs — 2, except at high magnetic fields. The Fermi 
wavevector for all 2D systems is given simply by filling up 
the non-interacting momentum eigenstates upto q = kp. 



\q\<kf 



(2^ 



k-p = 



Airn 



(1.9) 



where n is the 2D carrier density in the system. Unless 
otherwise stated, we will mostly consider electron sys- 
tems (or the conduction band side of MLG and BLG). 
Typical experimental values of n « 10^ cm~^ to 5 x 
10"^^ cm~^ are achievable in graphene and Si MOSFETs, 
wheras in GaAs-based 2DEG systems n 10^ cm~^ to 
5 X 10" cm-2. 



1. Interaction parameter ra 

The interaction parameter, also known as the Wigner- 
Seitz radius or the coupling constant or the effective fine- 
structure constant, is denoted here by r^, which in this 
context is the ratio of the average inter-electron Coulomb 
interaction energy to the Fermi energy. Noting that 
the average Coulomb energy is simply (V) = e^/(K(r)), 
where (r) = (tttt,)"^/^ is the average inter-particle sepa- 
ration in a 2D system with n particles per unit area, and 
K is the background dielectric constant, we obtain rg to 



for MLG and 



for BLG and 2DEG. 



nn) (BLG and 
(GaAs 2DEG). 



A note of caution about the nomenclature is in or- 
der here, particularly since we have kept the degeneracy 
factors gs5v in the definition of the interaction param- 
eter. Putting gsgv — 4, the usual case for MLG, BLG 
and Si 100 2DEG, and g^gv = 2 for GaAs 2DEG, we get 

= e'^/{KhvF) (MLG), = 2me'^/{Kh'' 
Si 100 2DEG), and r, = me^/{Kh^y/¥h 
The traditional definition of the Wigner-Seitz radius for 
a metallic Fermi liquid is the dimensionless ratio of the 
average inter-particle separation to the effective Bohr ra- 
dius Ob = kIt? / [rae^) . This gives for the Wigner-Seitz 
radius r^s ^ me'^ / {kH^ y/¥n) (2DEG and BLG), which 
differs from the definition of the interaction parameter Tg 
by the degeneracy factor gsgv/2. We emphasize that the 
Wigner-Seitz radius from the above definition is meaning- 
less for MLG because the low-energy linear dispersion im- 
plies a zero effective mass (or more correctly, the concept 
of an effective mass for MLG does not apply). For MLG, 
therefore, an alternative definition, widely used in the lit- 
erature defines an "effective fine structure constant (a)" 
as the coupling constant a = /{nhvY), which differs 
from the definition of rg by the factor ^5s5v/2. Putting 
s/g^g^ = 2 for MLG gives the interaction parameter rg 
equal to the effective fine structure constant a, just as 
setting gsgv = 2 for GaAs 2DEG gave the interaction pa- 
rameter equal to the Wigner-Seitz radius. Whether the 
definition of the interaction parameter should or should 
not contain the degeneracy factor is a matter of taste, 
and has been discussed in the literature in the context of 



2D semiconductor systems (Das Sarma et al. 2009). 



A truly significant aspect of the monolayer graphene 
interaction parameter, which follows directly from its 
equivalence with the fine structure constant definition, is 
that it is a carrier density independent constant, unlike 
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Ts parameter for the 2DEG (or BLG), which increases 
with decreasing carrier density as n~^^'^. In particu- 
lar the interaction parameter for MLG is bounded, i.e., 
< rg < 2.2, since 1 < k < oo, and as discussed ear- 
Her Up ~ 10® cm/s is set by the carbon hopping pa- 
rameters and lattice spacing. This is in sharp contrast 
to 2DEG systems where rg ~ 13 (for electrons in GaAs 
with n « 10^ cm~^) and « 50 (for holes in GaAs with 



n ~ 2 X 10^ cm ^) have been reported ( 


Das Sarma et al. 


2005 


Huang et al.\ |2006 


Manfra et al. 


2007 





Monolayer graphene is thus, by comparison, always a 
fairly weakly interacting system, while bilayer graphene 
could become a strongly interacting system at low carrier 
density. We point out, however, that the real low den- 
sity regime in graphene (both MLG and BLG) is dom- 
inated entirely by disorder in currently available sam- 
ples, and therefore a homogeneous carrier density of 
n < lO^'' cm^^ (10^ cm~^) is unlikely to be accessible 
for gated (suspended) samples in the near future. Us- 
ing the BLG effective mass m = 0.03 We, we get the 
interaction parameter for BLG: Tg ~ 68.5/ (k-\/^), where 
n — cm"^. For a comparison, the Tg parame- 

ters for GaAs 2DEG (k = 13, m* = 0.67 TOc), Si 100 on 
Si02 (k = 7.7, m* = 0.19 me,5v = 2) are rg w and 
Tg ^ respectively. 

For the case when the substrate is Si02, k = (Ksi02 + 
l)/2 « 2.5 for MLG and BLG we have Tg w 0.8 and 
Tg w 27.4/(Vn) respecitvely. In vacuum, k = 1 and 
rg « 2.2 for MLG and rg « 68.5/{V^) for BLG. 



2. Thomas-Fermi screening wavevector gxp 

Screening properties of an electron gas depend on the 
density of states, Dq at the Fermi level. The sim- 
ple Thomas-Fermi theory leads to the long wavelength 
Thomas-Fermi screening wavevector 



<Ztf 



-Dn. 



(1.10) 



The density independence of long wavelength screening 
in BLG and 2DEG is the well known consequence of the 
density of states being a constant (independent of en- 
ergy), whereas the property that gxp ~ fcp ~ n^^'^ in 
MLG, is a direct consequence of the MLG density of 
states being linear in energy. 

A key dimcnsionless quantity determining the charged 
impurity scattering limited transport in electronic ma- 
terials is Qs — QTF/kp which controls the dimensionless 
strength of quantum screening. From Table |T] we have 

~ n° for MLG and qs ~ n'^/^ for BLG and 2DEG. Us- 
ing the usual substitutions gsgv = 4(2) for Si 100 (GaAs) 
based 2DEG system, and taking the standard values of 
m and k for graphene-Si02, GaAs-AlGaAs and Si-Si02 



structures, we get (for n = n/lQ^'^ cm ^) 

MLG: qs ~ 3.2, BLG: q^ w 54.8/v^, (1.11a) 
n-GaAs: qs ~ 8/Vh, p-GaAs: qs « 43/\/^(l.llb) 

We point out two important features of the simple 
screening considerations described above: (i)ln. MLG, qs 
being a constant implies that the screened Goulomb inter- 
action has exactly the same behavior as the unscreened 
bare Coulomb interaction. The bare 2D Coulomb in- 
teraction in a background with dielectric constant k 
is given by v{q) = l-Ke^jinq) and the corresponding 
long-wavelength screened interaction is given by M(g) — 
I («:(g-|-gTF))- Putting 5 = fcp in the above equation. 



we get, w(g) ~ (fcp + gTF) 



fcp (l + ^Tp/fep)"^ K 



for MLG. Thus, in MLG, the functional dependence of 
the screened Coulomb scattering on the carrier density 
is exactly the same as unscreened Coulomb scattering, a 
most peculiar phenomenon arising from the Dirac linear 
dispersion, (ii) In BLG (but not MLG, see above) and 
in 2DEG, the effective screening becomes stronger as the 
carrier density decreases since qs = ^tf/^f n~^/^ — >■ 
00 (0) as n — > (00). This counter-intuitive behavior 
of 2D screening, which is true for BLG systems also, 
means that in 2D systems, effects of Coulomb scattering 
on transport properties increases with increasing carrier 
density, and at very high density, the system behaves as 
an unscreened system. This is in sharp contrast to 3D 
metals where the screening effect increases monotonically 
with increasing electron density. 

Finally in the context of graphene, it is useful to give 
a direct comparison between screening in MLG versus 
screening in BLG: 9tf^/9tf^ ~ 16/v^, showing that as 
carrier density decreases, BLG screening becomes much 
stronger than SLG screening. 



3. Plasmons 

Plasmons are self-sustaining normal mode oscillations 
of a carrier system, arising from the long-range nature 
of the inter-particle Coulomb interaction. The plasmon 
modes are defined by the zeros of the corresponding fre- 
quency and wavevector dependent dynamical dielectric 
function. The long wavelength plasma oscillations are 
essentially fixed by the particle number (or current) con- 
servation, and can be obtained from elementary consid- 
erations. We write down the long-wavelength plasmon 
dispersion ujp 

/ 2 \ 1/2 

MLG: c^p(g^O) = f ^^V^7^Vl.l?a) 

BLG & 2DEG: LJp{q ^ 0) = q ■ (1.12b) 

A rather intriguing aspect of MLG plasmon dispersion 
is that it is non-classical (i.e. h appears explicitly in 
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Eq. 1.12 even in the long wavelength limit). This ex- 
plicit quantum nature of long wavelength MLG plas- 
mon is a direct manifestation of its linear Dirac like 
energy-momentum dispersion, which has no classical 



analogy (Das Sarma and Hwang 2009). 



4. Magnetic field effects 

Although magnetic field induced phenomena in 
graphene and 2D semiconductors (e.g. quantum Hall ef- 
fect and fractional quantum Hall effect) are briefly cov- 
ered in section |Vj we want to mention at this point a 
few elementary electronic properties in the presence of 
an external magnetic field perpendicular to the 2D plane 
leading to the Landau orbital quantization of the system. 

(a) Landau Level Energetics: The application of a strong 
perpendicular external magnetic field (B) leads to a com- 
plete quantization of the orbital carrier dynamics of all 
2D systems leading to the following quantized energy lev- 
els En . the so-called Landau levels 



MLG: En = sgn{n)vF^/2ehB\n\, 
with n = 0,±1,±2 
sgn(n) 



(1.13a) 



BLG: En 



V2 



{2\n\ + l){2eBvph) + Am^vj 



-^{2mvFY + 2(2|n| + l){2eBvlh){2mvlY ■ 



2DEG: En = (n+ 1/2) 



with n = 0,±1,±2, • • • . 
eBh^ 



{2eBvlhf , 
(1.13b) 



mc 

with n = 0, 1,2, • 



(1.13c) 

The hallmark of the Dirac nature of graphene is the ex- 



istence of a true zero-energy (n = in Eq. 1.13) Landau 
level, which is equally shared by electrons and holes. The 
experimental verification of this zero energy Landu level 
in graphene is definitive evidence for the long wavelength 



Dirac nature of the system 


(Miller et al. 


2009 


Novoselov 


et al. 2005a 


Zhang et al. 


20051. 



(b) Cyclotron Resonance: 

External radiation induced transitions between Lan- 
dau levels gives rise to the cyclotron resonance in a Lan- 
dau quantized system, which has been extensively stud- 



ied in 2D semiconductor ( Ando et al. 



systems ( jHenriksen et al. 



2010 2008 



1982 ) and graphene 



Jiang et al. 2007). 



The cyclotron resonance frequency in MLG and 2DEG is 
given by 



MLG: Wc 
2DEG: ujr 



VF^2ehB (\^ 
eB 



1 - 



(1.14a) 
(1.14b) 



For BLG, the cyclotron frequency should smoothly in- 
terpolate between the formula for MLG for very large 



that of the 2DEG for small n so that E^ ^ 2mvp (where 
m ss 0.033 is the approximate B = effective mass of the 
bilayer parabolic band dispersion). Experimental BLG 



cyclotron resonance studies (Henriksen et al. 2010) in- 
dicate the crossover from the quadratic band dispersion 
(i.e. 2DEG-like) for smaller q to the linear band disper- 
sion (i.e. MLG-like) at larger q seems to happen at lower 
values of q than that implied by simple band theory con- 
siderations. 

A particularly interesting and important feature of 
cyclotron resonance in graphene is that it is affected 
by electron-electron interaction effects unlike the usual 
parabolic 2DEG, where the existence of Kohn's theo- 
rem prevents the long wavelength cyclotron frequency 
from being renormalized by electron- electron interac- 
tions (I Ando et al.l 119821 IKohnl 119611). For further dis- 



cussion of this important topic, we refer the reader to the 



recent literature on the subject (Henriksen et al. 2010 



Shizuya 2010) 



(c) Zeeman splitting: 

In graphene, the spin splitting can be large since the 
Lande g-factor in graphene is the same {g — 2) as in 
vacuum. The Zeeman splitting in an external magnetic 
field is given by {^J-B■ Bohr magneton) E^, = gfJ^BB — 
0.12 B[T] meV, for 5 = 2 (MLG, BLG, Si 2DEG) and 
E^ = -0.03 B[T] meV for g = -0.44 (GaAs 2DEG). We 
note that the relative value of E^/E-p is rather small in 
graphene, E^/E^ ~ 0.01(B[T]/\/t) ^ 0.01 for B = 10 T 
and n — 10^^ cm"^. Thus the spin splitting is only 
1 percent even at high fields. Of course, the polariza- 
tion effect is stronger at low carrier densities, since Ep is 
smaller. 



n, so that E^ in Eq. 1.13 is much larger than 2mvp, to 



D. Intrinsic and extrinsic graphene 

It is important to distinguish between intrinsic and ex- 
trinsic graphene because gapless graphene (either MLG 
or BLG) has a charge neutrality point (CNP), i.e. the 
Dirac point, where its character changes from being 
electron-like to being hole-like. Such a distinction is not 
meaningful for a 2DEG (or BLG with a large gap) since 
the intrinsic system is simply an undoped system with no 
carriers (and as such is uninteresting from the electronic 
transport properties perspective). 

In monolayer and bilayer graphene, the ability to gate 
(or dope) the system by putting carriers into the conduc- 
tion or valence band by tuning an external gate voltage 
enables one to pass through the CNP where the chem- 
ical potential (Ep) resides precisely at the Dirac point. 
This system, with no free carriers at T = 0, and Ep pre- 
cisely at the Dirac point is called intrinsic graphene with 
a completely filled (empty) valence (conduction) band. 
Any infinitesimal doping (or, for that matter, any finite 
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TABLE II Electronic quantities for monolayer graphene. Note: The graphene Fermi velocity {vf — 10* cm/s) and the 
degeneracy factor g — g^g^ — 4, i.e. the usual spin degeneracy (gs = 2) and a valley degeneracy ((?„ — 2), are used in this table. 
Here h — 71/(10^" cm~^), and B, q, and a are measured in T, cm~^, and e^/h = 38.74 fiS (or h/e^ — 25.8 kSl) respectively. 



Quantity 



Fermi wave vector {kp) 
Thomas Fermi wave vector [qxp) 
Interaction parameter (r^) 
DOS at Ef {Do = D{Ep)) 
Fermi energy (Ef) 
Zeeman splitting (Ez) 
Cyclotron frequency (cuc) 
Landau level energy (E„) 

Plasma frequency (cjp(q)) 
Mobility (/i) 
Scattering time (r) 
Level broadening (P) 



Scale values 



1.77 X lO^V^ [cin-^] 
1.55 X 10^V^/k [cm~^] 
2.19/k 

1.71 X lO^Vfi [meV-^cm-2] 
11.65Vfi [meV] 
0.125 [meV] 
5.51 X IO^^VB [s^^] 
sgn(/) 36.29y/B|T[ [meV], I = 0, ±1, ±2 

5.80 X IQ-'^^Vfiq/K [meV] 
2.42 X 10" a/h [cm^Vs] 
2.83 X 10"^" a/Vfi [s] 
11.63\/fi/(j [mcV] 



temperature) makes the system "extrinsic" with elec- 
trons (holes) present in the conduction (valence) band. 
( Miiller et al. 2009 ) . Although the intrinsic system is a 



set of measure zero (since Ep has to be precisely at the 
Dirac point), the routine experimental ability to tune the 
system from being electron-like to to being hole-like by 
changing the external gate voltage, manifestly establishes 
that one must be going through the intrinsic system at 
the CNF. If there is an insulating regime in between, 
as there would be for a gapped system, then intrinsic 
graphene is not being accessed. 

Although it is not often emphasized, the great achieve- 



ment of Novoselov et al. { 2004 ) in producing 2D graphene 



in the laboratory is not just fabricati ng (jNovoselov et'of. 



2005b) and identifying (Ferrari 2007 Ferrari ei aL 2006) 



stable monolayers of graphene flakes on substrates, but 
also establishing its transport properties by gating the 
graphene device using an external gate which allows one 
to simply tune an external gate voltage and thereby con- 
tinuously controlling the 2D graphene carrier density as 
well as their nature (electron or hole). If all that could 
be done in the laboratory was to produce beautiful 2D 
graphene flakes, with no hope of doping or gating them 
with carriers, the subject of graphene would be many or- 
ders of magnitude smaller and less interesting. What led 
to the exponential growth in graphene literature is the 
discovery of gatable and density tunable 2D graphene in 
2004. 

Taking into account the quantum capacitance in 
graphene the doping induced by the external gate voltage 



Rossier et al. 2007): 



Vg is given by the relation ( Fang et al. 2007 Fernandez 



(1.15) 



of the electron charge and n g = f ( ^2^ ) . The sec- 
ond term on the r.h.s. of ( 1.15[ ) is analogous to the 
term due to the so-called quantum capacitance in reg- 
ular 2DEG. Notice that in graphene, due to the linear 
dispersion, contrary to parabolic 2D electron liquids, the 
quantum capacitance depends on Vg. For a background 
dielectric constant k sa 4 and gate voltages larger than 
few millivolts the second term on the r.h.s. of (1.15) 



where C is the gate capacitance, e the absolute value 



can be neglected for thicknesses of the dielectric larger 
than few angstroms. In current experiments on exfoliated 
graphene on Si02 the oxide is 300 nm thick and there- 
fore quantum-capacitance effects are completely negligi- 
ble. In this case, a simple capacitance model connects 
the 2D carrier density (n) with the applied external gate 
voltage n « CVg, where C « 7.2 x 10^° cm^VV for 
graphene on Si02 with roughly 300 nm thickness. This 
approximate value of the constant C seems to be pretty 
accurate, and the following scaling should provide n for 
different dielectrics 

n [IQio cm-2] = 7.2 x ^i^— K [V], (1.16) 
^ ^ 300 3.9 ^ ^ ^ ^ ^ 

where t is the thickness of the dielectric (i.e. the distance 
from the gate to the graphene layer) and k is the dielectric 
constant of the insulating substrate. 

It is best, therefore, to think of 2D graphene on Si02 
(see Fig. 1(e)) as a metal-oxide-graphene-field-effect- 
transistor (MOGFET) similar to the well known Si MOS- 
FET structure, with Si replaced by graphene where the 
carriers reside. In fact, this analogy between graphene 
and Si 100 inversion layer is operationally quite effec- 
tive: Both have the degeneracy factor gsg^ = 4 and both 
typically have Si02 as the gate oxide layer. The quali- 
tative and crucial difference is, of course, that graphene 
carriers are chiral, massless, with linear dispersion and 
with no band gap, so that the gate allows one to go 
directly from being n type to a p type carrier system 
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through the charge neutral Dirac point. Thus a graphene 
MOGFET is not a transistor at all (at least for MLG), 
since the system never becomes insulating at any gate 



voltage (Avouris et al. 2007). 



We will distinguish between extrinsic (i.e. doped) 
graphene with free carriers and intrinsic (i.e. undoped) 
graphene with the chemical potential precisely at the 
Dirac point. All experimental systems (since they are 
always at T 7^ 0) are necessarily extrinsic, but intrinsic 
graphene is of theoretical importance since it is a critical 
point. In particular, intrinsic graphene is a non-Fermi liq- 



uid in the presence of electron-electron interactions ( Das 



Sarma et al. 2007b), whereas extrinsic graphene is a 
Fermi liquid. Since the non-Fermi liquid fixed point for 
intrinsic graphene is unstable to the presence of any fi- 
nite carrier density, the non-Fermi liquid nature of this 
fixed point is unlikely to have any experimental impli- 
cation. But it is important to keep this non-Fermi liq- 
uid nature of intrinsic graphene in mind when discussing 
graphene's electronic properties. We also mention (see 
Sec. 



IV) 



that disorder, particularly long-ranged disor- 
der induced by random charged impurities present in the 
environment, is a relevant strong perturbation affecting 
the critical Dirac point, since the system breaks up into 
spatially random electron-hole puddles, thus masking its 
zero-density intrinsic nature. 



E. "Other topics" 

There are several topics that are of active current re- 
search that we could not cover in any depth in this review 
article. Some of these remain controversial and others are 
still poorly understood. Yet these subjects are of impor- 
tance, both in terms of fundamental physics and for the 
application of graphene for useful devices. Here we sketch 
the status of these important subjects at the time of writ- 
ing this article. For example, there have recently emerged 
several novel methods of fabricating graphene including 
chemical vapor deposition on nickel ( |Kim et al.\ |2009[ ) 



and copper (Li et al. 2009), as well as directly unzip- 



ping carbon nanotubes (Kosynkin et al. 2009 Sinitskii 



et al. 2009 ) and other chemical methods ( Jiao et al. 
2009). As of early 2010, all of these other fabrication 
processes are just in their infancy. The notable excep- 
tion is "epitaxial graphene" manufactured by heating SiC 
wafers, causing the Si atoms to desorb, resulting in sev- 



eral 


graphene layers at the surface (Berger et al. 2004 


2006 


|de Heer et aZ.||2010||Emtsev et al. 


2009 IFirst et al. 


2010 


) that are believed to be very weakly coupled and of 



very good quality ( 


Hass et al. 


2008 


Miller et al. 


2009 


Orlita et al. 


2008 


Rutter et al. 


20071. We note that 



graphene can be used as a component of more compli- 



cated structures by 


exploiting its spin 


Cho et al. 


(2007 


I; 


Han et al. 


( 


2009 


); 


Hill et al. 


(2006 


); Huertas- Hernando 


et al. 


( 


2006 


2009) 


Jozsa et al. 


( 


2009 


); 


Tombros et al. 



(2007) or valley Rycerz et al. (2007b) degeneracy or by 



patterning gates with a periodic super-potential (Brey 



and Fertig 2009a Park ei aZ. 2008). Graphene can also 



be made to superconduct by coupling it to superconduct- 



ing leads and through the proximity effect (Beenakker 



2006 2008 Du et al. 



2008a Heersche et al. 2007) or 



other novel proposals (Feigel'man et al. 2008 Lutchyn 



et al. 2008). This review could not cover any of these 



topics in any reasonable depth. 



1. Optical conductivity 



It was pointed out as early as 1994 by |Ludwig et al.\ 
( 1994 ) that if one examined the conductivity of Dirac 



Fermions in linear response theory, keeping a finite fre- 
quency, i.e. <y{uj) while taking the limit of zero temper- 
ature [T — > 0) and vanishing disorder (F — >• 0), then 
one obtained a universal and frequency independent op- 
tical conductivity (i.e. electrical conductivity at finite 
frequency) 



o'(w) = gsg-. 



Tie 
'~8h 



(1.17a) 



Ludwig et al. (1994) also noted that this result did not 
commute with the d.c. conductivity where one first took 
the limit — > and then F — > 0, in which case one 
obtained 



9sgv 



■nh 



(1.17b) 



These T — Q results apply to intrinsic graphene where 
E-p is precisely at the Dirac point. The crossover be- 
tween these two theoretical intrinsic limits remains an 



open problem ( Katsnelson 2006 Ostrovsky et al. 2006 1 . 



The optical conductivity (Eq. 1.17a) has been mea- 



sured experimentally both by infrared spectroscopy (Li 



et al.\ 2008 1 and by measuring the absorption of sus- 



pended graphene sheets (Nair et al. 2008). In the IR 



measurements, <j{uj) is close to the predicted univer- 
sal value for a range of frequencies 4000 cm~^ < uj < 
6500 cm^^. While in the absorption experiment, the 
attenuation of visible light through multilayer graphene 
scales as na per layer. The authors claimed that this was 
an accurate measurement of the fine structure constant a, 
and is a direct consequence of having cr(w) being a univer- 
sal and frequency independent constant. In some sense, it 
is quite remarkable that disorder and electron-electron in- 
teractions do not significantly alter the value of the opti- 
cal conductivity. This has attracted considerable theoret- 



ical interest ( Gusynin and Sharapov , 2006 Herbut et al. 



2008| |Katsnelson[ 
MacDonald' '2009 



2008; Kuzmenko et al, 2008; Min and 



Mishchenkol \M)T\ ,2009, |Peres anc 



Stauber, 2008; Peres et al., 2008;|Sheehy and Schmalian 



2009[ [Stauber et al. 2008a|), where it has been argued 
that it is a fortuitous cancellation of higher order terms 
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that explains the insensitivity of a{uj) to interaction ef- 
fects. We refer the reader to these works for detailed 
discussion of how interaction effects and disorder change 
cr{ui) from the universal value, although, a consensus is 
yet to emerge on whether these effects could be observed 
experimentally or how accurate a{uj) is for a measure of 



the fine structure constant ( Gusynin et al. 2009 Mak 



et al. 2008). 



2. Graphene nanoribbons 

It was realized in the very first graphene trans- 
port experiments that the finite minimum conductivity 



(Eq. 1.17) would be an obstacle for making a useful tran- 
sistor since there is no "off' state. One way to cir- 
cumvent this problem is to have a quasi- ID geometry 
that confines the graphene electrons in a strip of (large) 
length L and a finite (small) width W . The confinement 



gap typically scales as 1/W^ (Wakabayashi et oL) 19991, 



however this depends on the imposed boundary condi- 
tions. This is quite similar to carbon nanotubes (since 
a nanotube is just a nanoribbon with periodic boundary 
conditions). The nomenclature in graphene is slightly 
different from carbon nanotubes, where a zigzag-edge 
nanoribbon is similar to an armchair nanotube in that 
it is always metallic within the tight-binding approxima- 
tion. Similarly, an armchair nanoribbon is similar to a 
zigzag nanotube in that it can be either metallic or semi- 
conducting depending on the width. Early theoretical 
calculations ( Son et al. 2006a|b Yang et al. 2007 1 used 
a density functional theory to calculate the bandgap of 
armchair graphene nanoribbons and found that just like 
carbon nanotubes, the energy gaps come in three fami- 
lies that are all semiconducting (unlike the tight-binding 
calculation, which gives one of the families as metallic). 
Brey and Fertig ( 2006 ) showed that simply quantizing the 



Dirac Hamiltonian (the low energy effective theory) gave 
quantitatively similar results for the energy gaps as the 
tightbinding calculation, while Son et al. ( |2006a I showed 
that the density functional results could be obtained from 
the tightbinding model with some added edge disorder. 
By considering arbitrary boundary conditions, |Akhmerov| 
and Beenakker ( 2008 ) demonstrated that the behavior of 



the zigzag edge is the most generic for graphene nanorib- 
bons. These theoretical works gave a simple way to un- 
derstand the gap in graphene nanoribbons. 

However, the first experiments on graphene nanorib- 
bons (Han et al. 2007[ ) presented quite unexpected re- 
sults. As shown in Fig. [5] the transport gap for nar- 
row ribbons is much larger than that predicted by the- 
ory (with the gap diverging at widths « 15 nm), while 
wider ribbons have a much smaller gap than expected. 
Surprisingly, the gap showed no dependence on the ori- 
entation (i.e. zig-zag or armchair direction) as required 
by the theory. These discrepancies have prompted sev- 




FIG. 5 (Color online) (a) Graphene nanoribbon energy gaps 
as a function of width, adapted from Han et al. 



(20071. 



Four devices (P1-P4) were orientated parallel to each other 
with varying width, while two devices (D1-D2) were oriented 
along different crystallographic directions with uniform width. 
The dashed line is a fit to a phenomenological model with 
Eg = A/(W — W*) where A and W* are fit parameters. The 
inset shows that contrary to predictions, the energy gaps have 
no dependence on crystallographic direction. The dashed 
lines are the same fits as in the main panel, (b) Evidence for 
a percolation metal-insulator transition in graphene nanorib- 
bons, adapted from Adam et al. { 2008a!). Main panel shows 



graphene ribbon conductance as a function of gate voltage. 
Solid lines are a fit to percolation theory, where electrons and 
holes have different percolation thresholds (seen as separate 
critical gate voltages Vc). The inset shows the same data in 
a linear scale, where even by eye the transition from high- 
density Boltzmann behavior to the low-density percolation 
transport is visible. 



eral studies (Abanin and Levitov 



2008 



Adam et al. 



2008a', 'Areshkin et aL','2 007[|Basu a Z.','2008','Biel et al. 



2009a b; Chen et al, 2007 



Blanterl |2009| [SoTT et al 



|Dietl et al , 2009; Mar tin and 



[2007| IStampfer et'al.] |2009 



Todd et al 20091. In particular, Sols et al (2007) ar- 



gued that fabrication of the nanoribbons gave rise to 
very rough edges breaking the nanoribbon into a series 
of quantum dots. Coulomb blockade of charge transfer 



between the dots (Ponomarenko et al. 2008) explains 



the larger gaps for smaller ribbon widths. In a similar 



spirit, Martin and Blanter (20091 showed that edge dis- 



order qualitatively changed the picture from that of the 
disorder-free picture presented earlier, giving a localiza- 
tion length comparable to the sample width. For larger 
ribbons, Adam et al. (2008a) argued that charged im- 



purities in the vicinity of the graphene would give rise 
to inhomogeneous puddles so that the transport would 
be governed by percolation (as shown in Fig. [s] points 
are experimental data, and the solid lines, for both elec- 
trons and holes, show fits to cr ~ {V — Vc)" , where v 
is close to 4/3, the theoretically expected value for per- 
colation in 2D systems). The large gap for small rib- 
bon widths would then be explained by a dimensional 
crossover as the ribbon width became comparable to 
the puddle size. A numerical study including the effect 
of quantum localization and edge disorder was done by 
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Mucciolo et al. ( 2009 ) who found that a few atomic layers 
of edge roughness was sufficient to induce transport gaps 
appear that are approximately inversely proportional to 
the nanoribbon width. Two very recent and detailed ex- 



periments (Gallagher et al. 2010 Han et al. 20101 seem 



to suggest that a combination of these pictures might be 
at play (e.g. transport trough quantum dots that are 
created by the charged impurity potential), although as 
yet, a complete theoretical understanding remains elu- 
sive. The phenomenon of the measured transport gap 
being much smaller than the theoretical bandgap seems 
to be a generic feature in graphene, occuring not only in 
nanoribbons but also in biased bilayer graphene where 
the gap measured in transport experiments appears to 
be substantially smaller than the theoretically calculated 



and a corresponding carrier density of ~ 10^^ cm~^. To 
avoid the warping, it was proposed that one should use 
a top gate with the opposite polarity, but at the time 
of writing, this has yet to be demonstrated experimen- 
tally. Despite the limited variation in carrier density, sus- 
pended graphene has achieved a carrier mobility of more 
than 2 00, 000 cm^/Vs ( [Bolotin al| |2008a|b| |Du et al. 



2008b ) . Very recently, suspended graphene bilayers were 



demonstrated experimentally ( Feldman et al. 2009 ) 



bandgap (Oostinga et al. 2008) or even the measured 



optical gap (Mak et al. 2009 Zhang et al. 2009c) 



4. Many-body effects in graphene 

The topic of many-body effects in graphene is itself a 
large subject, and one that we could not cover in this 
transport review. As already discussed earlier, for intrin- 
sic graphene, the many-body ground state is not even 



3. Suspended graphene 



a Fermi liquid (Das Sarma et al. 2007b), an indication 
of the strong role played by interaction effects. Exper- 
imentally, one can observe the signature of many body 



effects in the compressibility (Martin et al. 20081 and 



Since the substrate affects both the morphology of using ARPES (j Bostwick et aL[ |2007[ |Zhou eVaL[|2007 |. 



graphene (Ishigami et al. 2007 Meyer et al. 2007 Stol 



yarova et al. 2007 ) as well as provides a source of impu- 



rities, it became clear that one needed to find a way to 



Away from the Dirac point, where graphene behaves as 
a normal Fermi liquid, the calculation of the electron- 
electron and electron-phonon contribution to the quasi- 



particle self-energy was studied by several groups (Bar- 



neath the graphene so that the graphene is suspended 
over a trench that is approximately 100 nm deep. As 
a historical note, we mention that suspended graphene 
without electrical contacts was made earlier by |Meyer| 
et al. (2007). Quite surprisingly, the suspended samples 



of the underlying substrate. The making of "suspended 
graphene" or "substrate-free" graphene was an important 


las et al.\ 2007| Calandra and Mauril |2007 Carbotte 


et al. 


|2010; Hwang and Das Sarma, 2008; Hwang et al. 


experimental milestone ( Bolotin et al. 2008a|b Du et al. 


2007c|dnPark et aL 2007i,2009,|Polini et aZ.,|2007 2008a 


2008b ) where after exfoliating graphene and making elec- 


Tse and Das Sarma 2007), and show reasonable agree- 


trical contact, one then etches away the substrate under- 


ment with experiments ( Bostwick et al. 2007 Brar et al. 



2010). For both bilayer graphene (Min et al. 2008b) and 



for double layer graphene (Min et al. 2008a I an insta- 



bility towards an exitonic condensate has been proposed. 
In general, monolayer graphene is a weakly interacting 
system since the coupling constant [r^ < 2) is never 
as prepared did not show much difference from unsus- large (jMuller et "aT] |2009[). In principle, bilayer graphene 



pended graphene, until after current annealing ( Barreiro 



et al. 2009 Moser et al. 2007). This suggested that 



most of impurities limiting the transport properties of 
graphene were stuck to the graphene sheet and not buried 
in the substrate. After removing these impurities by 
driving a large current through the sheet, the suspended 
graphene samples showed both ballistic and diffusive car- 
rier transport properties. Away from the charge neu- 
trality point, suspended graphene showed near-ballistic 
transport over hundreds of nm, which prompted much 



could have arbitrarily large coupling at low carrier den- 
sity where disorder effects are also important. We refer 
the reader to these original works for details on this vast 
and interesting subject. 



5. Topological insulators 



There is a deep connection between graphene and topo- 



theoretical interest ( 


Adam and Das Sarma, 2008b 


Fogler 


et al. 


2008a 


Miiller et al. 2009 


Stauber et al. 




2008b 


1. 



logical insulators (Kane and Mele 2005a Sinitsyn et al. 



One problem with suspended graphene is that only a 
small gate voltage iVg « 5 V) could be applied before the 
graphene buckles due to the electrostatic attraction be- 
tween the charges in the gate and on the graphene sheet, 
and binds to the bottom of the trench that was etched 
out of the substrate. This is in contrast to graphene on 
a substrate that can support as much as 14 ~ 100 V 



2006). Graphene has a Dirac cone where the "spin" 
degree of freedom is actually related to the sublattices 
in real space, whereas it is the real electron spin that 
provides the Dirac structure in the topological insula- 



tors (Hasan and Kane 2010) on the surface of BiSb and 
BiTe (iGhen et aLl |2009dl iHsieh et a/.l|2008p. Graphene 



is a weak topological insulator because it has two Dirac 
cones (by contrast, a strong topological insulator is char- 
acterized by a single Dirac cone on each surface); but 
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in practice, the two cones in graphene are mostly decou- 
pled and it behaves like two copies of a single Dirac cone. 
Therefore many of the results presented in this review, 
although intended for graphene, should also be relevant 
for the single Dirac cone on the surface of a topological 
insulator. In particular, we expect the interface trans- 
port properties of topological insulators to be similar to 
the physics described in this review as long as the bulk 
is a true gapped insulator. 



F. 2D nature of graphene 

As the concluding section of the Introduction, we ask: 
what precisely is meant when an electronic system is cat- 
egorized as 2D and how can one ensure that a specific 
sample/system is 2D from the perspective of electronic 
transport phenomena. 

The question is not simply academic, since 2D does 
not necessarily mean a thin film (unless the film is liter- 
ally one atomic monolayer thick as in graphene, and even 
then, one must consider the possibility of the electronic 
wavefunction extending somewhat into the third direc- 
tion) . Also the definition of what constitutes a 2D may 
depend on the physical properties or phenomena that one 
is considering. For example, for the purpose of quantum 
localization phenomena, the system dimensionality is de- 
termined by the width of the system being smaller than 
the phase coherence length (or the Thouless length). 
Since could be very large at low temperature, metal 
films and wires can respectively be considered 2D and 
ID for localization studies at ultra-low temperature. For 
our purpose, however, dimensionality is defined by the 
3D electronic wavefunction being "free" plane-wave like 
(i.e. carrying a conserved 2D wavevector) in a 2D plane, 
while being a quantized bound state in the third dimen- 
sion. This ensures that the system is quantum mechani- 
cally 2D. 

Considering a thin film of infinite (i.e. very large) di- 
mension in the xy-plane, and a finite thickness w in the 
z-direction, where w could be the typical confinement 
width of a potential well creating the film, the system is 
considered 2D if Ap = > w, For graphene, we have 

Af « (350/v^) nm, where n ~ n/{10^^ cm ^), and since 
w « 0.1 nm to 0.2 nm (the monolayer atomic thickness), 
the condition Af 3> w is always satisfied, even for un- 
physically large n = 10^'' cm~^. 

Conversely, it is essentially impossible to create 2D 
electronic systems from thin metal films since the very 
high electron density of metals, provides Ap ~ 0.1 nm, so 
that even for a thickness w w 1 nm (the thinnest metal 
film that one can make), Ap < w, making them effectively 
3D. By virtue of the much lower carrier densities in semi- 
conductors, the condition \p > w can be easily satisfied 
for w = 5 nm to 50 nm for n = 10^ cm~^ to 10^^ cm~^, 
making it possible for 2D semiconductor systems to be 



readily available since confinement potentials with width 
ss 10 nm can be implemented by external gate voltage or 
band structure engineering. 



We now briefly address the question of the experimen- 
tal verification of the 2D nature of a particular system 
or sample. The classic technique is to show that the or- 
bital electronic dynamics is sensitive only to a magnetic 
field perpendicular to the 2D plane (i.e. Bz) (Practi- 
cally, there could be complications if the spin proper- 
ties of the system affect the relevant dynamics, since the 
Zeeman splitting is proportional to the total magnetic 
field). Therefore, if either the magnetoresistance oscilla- 
tions (Shubnikov-de Hass effect) or cyclotron resonance 
properties depend only on B^, then the 2D nature is es- 
tablished. 



Both of these are true in graphene. The most defini- 
tive evidence for 2D nature, however, is the observation 
of the quantum Hall effect, which is a quintessentially 
2D phenomenon. Any system manifesting an unambigu- 
ous quantized Hall plateau is 2D in nature, and therefore 
the observation of the quantum Hall effect in graphene 



in 2005 by Novoselov et al. (2005a) and Zhang et al. 



(2005) absolutely clinched its 2D nature. In fact, the 



quantum Hall effect in graphene persists to room temper- 



ature ( Novoselov et al. 2007 ) , indicating that graphene 



remains a strict 2D electronic material even at room tem- 
perature. 



Finally, we remark on the strict 2D nature of graphene 
from a structural viewpoint. The existence of finite 2D 
flakes of graphene with crystalline order at finite temper- 
ature does not in any way violate the Hohenberg-Mermin- 
Wagner-Coleman theorem which rules out the breaking 
of a continuous symmetry in two dimensions. This is be- 
cause the theorem only asserts a slow power law decay of 
the crystalline (i.e. positional order) correlation with dis- 
tance, and hence, very large flat 2D crystalline flakes of 
graphene (or for that matter, of any material) are man- 
ifestly allowed by this theorem. In fact, a 2D Wigner 
crystal, i.e. a 2D hexagonal classical crystal of electrons 
in a very low-density limit, was experimentally observed 



more than thirty years ago (Grimes and Adams 1979) 



on the surface of liquid He4 (where the electrons were 
bound by their image force). A simple back of the enve- 
lope calculation shows that the size of the graphene flake 
has to be unphysically large for this theorem to have any 



effect on its crystalline nature ( Thompson-Flagg et al. 



2009). There is nothing mysterious or remarkable about 



having finite 2D crystals with quasi-long-range positional 
order at finite temperatures, which is what we have in 2D 
graphene flakes. 
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II. QUANTUM TRANSPORT 
A. Introduction 

The phrase "quantum transport" usually refers to the 
charge current induced in an electron gas in response to 
a vanishing external electric field in the regime where 



quantum interference effects are important (Akkermans 



and Montambaux 2007 Rammer 1988). This is relevant 



at low temperatures where the electrons are coherent 
and interference effects are not washed out by dephas- 
ing. Theoretically, this corresponds to the systematic 
application of diagrammatic perturbation theory or field 
theoretic techniques to study how quantum interference 
changes the conductivity. For diffusive transport in two 
dimensions (including graphene), to lowest order in this 
perturbation theory, interference can be neglected, and 
one recovers the Einstein relation ctq — e^D{Ey)'D, where 
D{Ef) is the density of states at Ep and V = VpT/2 is 
the diffusion constant. This corresponds to the classical 
motion of electrons in a diffusive random walk scattering 
independently off the different impurities. 

Since the impurity potential is typically calculated us- 
ing the quantum-mechanical Born approximation, this 
leading order contribution to the electrical conductivity 
is known as the semi-classical transport theory and is 
the main subject of Sec. |ni.A| below. Higher orders in 
perturbation theory give "quantum corrections" to this 
semi-classical result i.e. a = (Tq + Sa, where Sa << a. 
In some cases these corrections can be divergent; a result 
that simultaneously implies a formal breakdown of the 
perturbation theory itself, while suggesting a phase tran- 
sition to a non-perturbatively accessible ground state. 

For example, it is widely accepted that quantum inter- 
ference between forward and backward electron trajec- 
tories is the microscopic mechanism responsible for the 



Anderson metal-insulator transition (Abrahams et al. 



1979). For this reason the leading quantum correction 
to the conductivity is called "weak localization" and is 
interpreted as the precursor to Anderson localization. 

Weak localization is measured experimentally by us- 
ing a magnetic field to break the symmetry between the 
forward and backward trajectories causing a change in 
the resistance. In this case the zero field conductivity 
a{B = 0) = (jQ + 6(7 includes the quantum corrections 
while (t{B > B*) — (Tq has only the semiclassical con- 
tribution; {B* is the approximately the magnetic field 
necessary to thread the area of the sample with one flux 
quantum.) 

The second hallmark of quantum transport is meso- 
scopic conductance fluctuations. If one performed the 
low temperature magnetotransport measurement dis- 
cussed above, one would notice fluctuations in the mag- 
netoresistance that would look like random noise. How- 
ever, unlike noise, these traces are reproducible and are 
called magneto-fingerprints. These magneto-fingerprints 



depend on the positions of the random impurities as seen 
by the electrons. Annealing the sample relocates the im- 
purities and changes the fingerprint. The remarkable fea- 
ture of these conductance fluctuations is that their mag- 
nitude is universal (depending only on the global sym- 
metry of the system) and notwithstanding the caveats 
discussed below, are completely independent of any mi- 
croscopic parameters such as material properties or type 
of disorder. 

While the general theory for weak localization and 
universal conductance fluctuations is now well estab- 



lished (Lee and Ramakrishnan 1985), in Section II. C. 3 



we discuss its application to graphene. 

The discussion so far has concerned diffusive transport, 
in what follows, we also consider the ballistic properties 
of non-interacting electrons in graphene. Early studies on 
the quantum mechanical properties of the Dirac Hamil- 
tonian revealed a peculiar feature - Dirac carriers could 



not be confined by electrostatic potentials (Klein 1929). 

An electron facing such a barrier would transmute into 
a hole and propagate through the barrier. In Section [lI.B| 
we study this Klein tunneling of Dirac carriers and dis- 
cuss how this formalism can be used to obtain graphene's 
ballistic universal minimum conductivity. There is no 
analog of this type of quantum-limited transport regime 
in two dimensional semiconductors. The "metallic na- 
ture" of graphene gives rise to several interesting and 
unique properties that we explore in this section includ- 
ing the absence of Anderson localization for Dirac elec- 
trons and a metal-insulator transition induced by atomi- 
cally sharp disorder (such as dislocations). We note that 
many of the results in this section can be also obtained us- 



ing field-theoretic methods ( Altland 


2006 


Fradkin 


1986 


Fritz et al. |2008 


Ludwig et al.\ 1994 Ostrovsky et al. 


2006 Ryu et al. 


2007bj iSchuessler et al. 2009). 



B. Ballistic transport 

1. Klein tunneling 

In classical mechanics, a potential barrier whose height 
is greater than the energy of a particle will confine that 
particle. In quantum mechanics, the notion of quantum 
tunneling describes the process whereby the wavefunc- 
tion of a non-relativistic particle can leak out into the 
classically forbidden region. However, the transmission 
through such a potential barrier decreases exponentially 
with the height and width of the barrier. For Dirac par- 
ticles, the transmission probability depends only weakly 
on the barrier height, approaching unity with increas- 
ing barrier height (Katsnelson et al. 2006). One can 



understand this effect by realizing that the Dirac Hamil- 
tonian allows for both positive energy states (called elec- 
trons) and negative energy states (called holes) . Whereas 
a positive potential barrier is repulsive for electrons, it is 
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attractive for holes (and vice versa). For any potential 
barrier one needs to match the electron states outside the 
barrier with the hole states inside the barrier. And since 
the larger the barrier, the greater the mode matching be- 
tween electron and hole states, the greater the transmis- 
sion. For an infinite barrier, the transmission becomes 



perfect. This is called Klein tunneling (Klein 19291 



By solving the transmission and reflection coefficients 



for both the graphene p-n junction ( Cheianov and Fal'ko 



2006b[ [Low and Appenzeller[ |2009[ ) and the p-n-p junc- 
tion ( Katsnelson et al. 2006 ) , it was found that for 



graphene the transmission at an angle normal to the bar- 
rier was always perfect (although there could be some re- 
flection at other angles). This can be understood in terms 
of pseudospin conservation. At normal incidence, the in- 
coming electron state and the reflected electron state are 
of opposite chirality resulting in vanishing probability for 
reflection. 

At finite angles of incidence, the transmission depends 
on how sharp the barrier is. In the limit of a per- 
fectly sharp step, the transmission probability is deter- 
mined only by pseudospin conservation and given by 
'7step(^) = cos^O. For a smooth variation in the elec- 
trostatic potential that defines the p-n junction (char- 
acterized by a length-scale ^), the transmission proba- 
bility was shown by Cheianov and Fal'ko (2006b) to be 
T^{e) = exp[-7r(A:FC)sin^6')]. This implies that for both 
sharp and smooth potential barriers, a wavepacket of 
Dirac Fermions will coUimate in a direction perpendicular 
to the p-n junction. One can estimate the conductance 
of a single p-n junction (of width W) to be 



4e2 



dO , , fcpC>>i 2e2 



-W. (2.1) 



Although the conductance of smooth p-n junctions are 
smaller by a factor of ^/kpS, compared to sharp ones, 
this result suggests that the presence of p-n junctions 
would make a small contribution to overall resistivity of a 



graphene sample (see also Sec. IV. D below), i.e. graphene 
p-n junctions are essentially transparent. 

The experimental realization of p-n junctions came 



shortly after the theoretical predictions (Huard et al. 


2007 


Lemme et al. 


2007 


Ozyilmaz et aLj|2007,,Williams 


et al. 


2007). At zero magnetic field the effect of creating 



a p-n junction was to modestly change the device resis- 
tance. More dramatic was the change at high magnetic 
field in the quantum Hall regime (see Sec. [v| below). 

More detailed calculations of the zero field conductance 
of the p-n junction were performed taking into account 
the effect of non-linear electronic screening. This tends to 
make the p-n junction sharper, and f or r., <g; 1, increas- 



ing the conductance by a factor (IZhang and Fogler 



tion of p-n junctions to the total resistance. The effect 



studied how the p-n junction resistance changed from its 
ballistic value in the absence of disorder to the diffusive 
limit with strong disorder. More recently, [Rossi et al.\ 
(2010a) used a microscopic model of charged impurities 



to calculate the screened disorder potential and solved 
for the conductance of such a disordered n-p-n junction 
numerically. The broad oscillations visible in Fig. [6] arise 




-10 12 

An [10^^ cm'^] 

tg 



■ (b) 






■ 








• 











5 10 15 20 

Aspect Ratio W/L 



FIG. 6 (Color online), (a) Disorder averaged resistance as a 
function of top gate voltage for a fixed back gate density nbg = 
5 X lO^^cm"^ and several values of the impurity density (from 
bottom to top Tiimp ~ 0,1,2.5,5,10, and 15 x lO^'^cm"^). 
Results were obtained using 10'^ disorder realizations for a 
square samples of size W = L = IQQ nm in presence of a 
top gate placed in the middle of the sample 10 nm above 
the graphene layer, 30 nm long and of width W. The charge 
impurities were assumed at a distance d = 1 nm and the 
uniform dielectric constant k was taken equal to 2.5. Adapted 
from (Rossi et al. 2010a I. (b) Solid line is Eq. 2.4 for armchair 
boundary conditions showing the aspect ratio dependence of 



the Dirac point ballistic conductivity ( Tworzydlo et al. 2006 1 . 



For W ^ L the theory approaches the univ ersal value Ae^ /nh . 



Circles show experimental data taken from Miao et al. ( 2007 1 



and squares show the data from Danneau et al. (20081. Inset: 
Illustration of the configuration used to calculate graphene's 
universal minimum conductivity. For Vg > 0, one has a p- 
p-p junction, while for Vg < 0, one has a p-n-p junction. 
The figure illustrates the ballistic universal conductivity that 
occurs at the transition between the p-p-p and p-n-p junctions 



from resonant tunneling of the few modes with smallest 
transverse momentum. These results demonstrate that 
the signatures of the Klein tunneling are observable for 
impurity densities as high as 10^^ cm~^ and would not 
be washed away by disorder as long as the impurity lim- 
ited mean-free-path is greater than length of the mid- 
dle region of the opposite polarity. This implies that at 
zero magnetic field the effects of Klein tunneling are best 
seen with a very narrow top-gate. Indeed, recent ex- 



periments have succeeded in using an "air-bridge" ( Gor- 



bachev et al. 2008 Liu et al. 2008 ) or very narrow top 



gates ( Slander et al. 2009 Young and Kim 20091. The 



[2008), and thereby further reducing the overall contribu- et al. (2010a 



observed oscillations in the conductivity about the semi- 
classical value are in good agreement with theory of |Rossi| 



of disorder was examined by Fogler et al. (2008b) who 



There is a strong similarity between the physics of 
phase-coherent ballistic trajectories of electrons and that 
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of light waves that is often exploited (Cheianov et al. leads and V{x) 



2007a Ji et al. 2003 Shytov et al. 2008). In particular, of Vc 



Vg for 



< X < L. For the case 



Liang et al. (2001 ) demonstrated that one could construct 



a Fabry-Perot resonator of electrons in a carbon nan- 
otube. This relies on the interference between electron 
paths in the different valleys K and K' . The same physics 
has been seen in "ballistic" graphene, where the device 
geometry is constructed such that the source and drain 
electrodes are closer than the typical electronic mean- 



free-path (Cho and Fuhrer 2009 Miao et al. 2007) 



2. Universal quantum limited conductivity 

An important development in the understanding of 
graphene transport is to use the formalism of Klein tun- 
neling to address the question of graphene's minimum 



conductivity (Katsnelson 2006 Tworzydlo et al. 2006). 



This of course considers non-interacting electrons at zero 
temperature, and in the limit of no disorder. As shown 
in inset of Fig. [6] (b) , the insight is to consider the source- 
graphene-drain configuration as the n-p-n or n-n-n junc- 
tion i.e. the leads are heavily electron doped, while the 
graphene sheet in the middle could be electron doped, 
hole doped or be pinned at the Dirac point with zero 
doping. Since there is no disorder, the electronic mean 
free path is much longer than the distance between the 
source and drain {£ ^ L). It is this situation we have in 
mind when we talk of graphene's "ballistic conductivity" . 

For a non-Dirac metal, at finite carrier density, the ab- 
sence of scattering would imply that the semi-classical 
electrical conductivity is infinite, since there is nothing 
to impede the electron motion. However, the conduc- 
tance would then vanish as the carrier density is tuned 
to zero. This metal-insulator transition will be discussed 
in more details later in the context of two dimensional 
semiconductors . 

The situation is quite different for graphene. From 
studying the Klein tunneling problem, we already know 
that both the n-p-n junction and the n-n-n junction have 
finite transmission coefficients. The interesting question 
is: what is the tunnelling at the precise point where the 
junction changes from a n-p-n junction to the n-n-n junc- 
tion? The conductivity at this transition point would 
then be the quantum limited (ballistic) conductivity of 
graphene at the Dirac point. 

The solution is obtained by finding the transmission 
probabilities and obtaining the corresponding ballistic 
conductivity. This is analogous to the quantum mechan- 
ics exercise of computing the transmission through a po- 
tential barrier, but now for relativistic electrons. Using 
the non-interacting Dirac equation 



[hvY-cr ■ k + eV{x)] ^{v) = e^'(r), 



(2.2) 



with the boundary conditions corresponding to: V{x < 
0) = V{x > L) = Voo to represent the heavily doped 



oo and at the Dirac point {Vg 



0), 



the transmission probability (i.e. the square of the 
transmission amplitude) is given by purely evanescent 



modes (Tworzydlo et al. 2006| ) 



1 

cosh(9„i) 



(2.3) 



This is in contrast to the non-relativistic electrons (i.e. 
with the usual parabolic dispersion), where for fixed 
Qn, the analogous calculation gives vanishing transmis- 
sion probability T„ 1/V^oo- The remaining subtle 
point is determining the form transverse wave-vector 
Qn. While it is clear that qn ^ nW~^ for large n, the 
choice of the boundary condition changes the precise re- 
lation e.g. Qn = mr/W for metallic armchair edges and 
q„ = (n ± l/3)7r/H^ for semiconducting armchair edges. 
FoUowing pfworzydlo et al] ( |2006[ ) , we use twisted bound- 
ary conditions ^'(j/ = 0) = o'x'^iy — 0) and '^{y — 
W) = -(Jx'ifiy = W) which gives ^ {n + l/2)Tr/W. 
This boundary condition is equivalent to having massless 
Dirac Fermions inside the strip of width W , but infinitely 
massive Dirac Fermions outside of the strip, thereby con- 



fining the electrons (Ryu et al. 2007a) 



The Landauer conductivity is then given by 



T 2 



W 

Ae' 



(2.4) 



n=0 



h :^^W cosh^ [T:{n + 1/2)L/W] tt/i ' 



E 



Since at the Dirac point (zero energy) there is no energy 
scale in the problem, the conductivity (if finite) can only 
depend on the aspect ratio L/W. The remarkable fact 



is that for W ^ L, the sum in Eq. 2.4 converges to a 
finite and universal value - giving for ballistic minimum 
conductivity a^ain — 4e^/7r/i. This result also agrees with 
that obtained using linear response theory in the limit of 
vanishing disorder suggesting that the quantum mechan- 
ical transport through evanescent modes between source 
and drain (or equivalently, the transport across two p-n 
junctions with heavily doped leads), is at the heart of 
the physics behind the universal minimum conductivity 
in graphene. 

Miao et al] ( |2007[ ) and [Danneau et al] ( |2008[ ) have 



probed this ballistic limit experimentally using the two- 
probe geometry. Their results, shown in Fig. [6](b), are 
in good agreement with the theoretical predictions. Al- 
though, it is not clear what role contact r esistance ([Blake 



2om 



Datta 



eTai] [20091 [Blanter and Martin', fSOOTl ICayssol et aT. 



Giovannetti et al.\ 



2009 Huard et al. 



2008; Golizadeh-Mojarad anc 



in these 2-probe measurements. 



2008||Lee et a/.[[2008b[ ) played 
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3. Shot noise 



Shot noise is a type of fluctuation in electrical cur- 
rent caused by the discreteness of charge carriers and 
from the randomness in their arrival times at the detec- 
tor or drain electrode. It probes any temporal correla- 
tion of the electrons carrying the current, quite distinct 
from "thermal noise" (or Johnson-Nyquist noise) which 
probes their fluctuation in energy. Shot noise is quanti- 
fled by the dimensionless Fano factor T , deflned as the 
ratio between noise power spectrum and the average con- 



C. Quantum interference effects 

1. Weak antilocalization 

Over the past 50 years, there has been much progress 
towards understanding the physics of Anderson localiza- 



ductance. Scattering theory gives (Biittiker 1990) 



(2.5) 



Some well known limits include = 1 for "Poisson noise" 
when T„ ^ 1 (e.g. in a tunnel junction), and = 1/3 



for disordered metals (Beenakker and Biittiker 1992 1. 



For graphene at the Dirac point, we can use Eq. |2.3| 
to get T 1/3 for » L ([Tworzydlo eTal] |2006|). 



tion (for a recent review, see Evers and Mirlin (2008)). 
Single particle Hamiltonians are classified according to 
their global symmetry. Since the Dirac Hamiltonian (for 
a single valley) Ji — hvpcr-h is invariant under the trans- 
formation H — ay'H*ay (analogous to spin- rotation sym- 
metry (SRS) in pseudospin space) it is in the All class 
(also called the symplectic Wigner-Dyson class). The 
more familiar physical realization of the symplectic class 
is the usual disordered electron gas with strong spin-orbit 
coupling. 

2m 

V„k,/3k' = 14_k' - «V^k-k'(k' X k) • a^p, (2.6) 

where a and (3 are (real) spin indices, cr^^ a vector of 
Pauli matrices. Notice that this Hamiltonian is also in- 
variant under SRS, % 



ay'H*ay. We have 



One should emphasize that obtaining the same numerical 
value for the Fano factor T = 1/3 for "ballistic" quan- 
tum transport in graphene as that of diffusive transport 
in disordered metals could be nothing more than a coin- 



S{q 



2ttvt 



lysoyso 



(5(q 



(2.7) 



cidence ( Dragomirova et al. 2009). Cheianov and Fal'ko rection to the conductivity is positive 



It was shown by Hikami et al. { 1980 ) that when the clas- 
sical conductivity is large (ctq ^ e^//i) , the quantum cor- 



( |2006b[ ) found that the shot noise of a single p-n junction 
1/2 that is numerically quite close to 1/3. 



was T = 1 ^ 

Since several different mechanisms all give J- ~ 1/3, 
this makes shot noise a complicated probe of the un- 
derlying physical mechanism. Recent numerical stud- 
ies by 



6a 



nh 



\n{L/e). 



(2.8) 



Equivalently, one can define a one parameter scaling 



function (Abrahams et al. 1979) 



Lewenkopf et al. (2008); San- Jose et al. (20071; 



Sonin 2008[ 2009|) treated the role of disorder to exam- 

1/3 in ballistic graphene 
Within the 



/3(a) 



ine the crossover from the 
to the diffusive regime (see also Sec. IV. B) 
crossover, or away from the Dirac point, the Fano fac- 
tor is no longer universal and shows disorder dependent 
deviations. The experimental situation is less clear. Dan- 
neau et al. (20081 measured the Fano factor decrease 



diner 
d\nL' 



(2.9) 



from « 1/3 with increasing carrier density to claim 
agreement with the ballistic theory. While [DiCarlo et al.\ 
( 2008 ) found that was mostly insensitive to carrier type 



and density, temperature, aspect ratio and the presence 
of a p-n junction, suggesting diffusive transport in the 
dirty limit. 

Since shot noise is, in principle, an independent probe 
of the nature of the carrier dynamics, it could be used as 
a separate test of the quantum-limited transport regime. 
However, in practice, the coincidence in the numerical 
value of the Fano factor with that of diffusive transport 
regime makes this prospect far more challenging. 



where for the symplectic class, it follows from Eq. |2.8| that 
/3((t) = l/(7rcr) for large a. To have /3 > means that the 
conductivity increases as one goes to larger system sizes 
or adds more disorder. This is quite different from the 
usual case of an Anderson transition where a negative /3- 
function means that for those same changes, the system 
becomes more insulating. 

Since perturbation theory only gives the result for 
f3{a :s> e'^/h), the real question becomes what happens 
to the /? function at small a. If the (3 function crosses 
zero and becomes negative as cr —> 0, the system exhibits 
the usual Anderson metal-insulator transition. Numeri- 



cal studies of the Hamiltonian (Eq. 2.6), show that for 
the spin-orbit system, the /3-function vanishes at a* « 1.4 
and below this value, the quantum correction to the clas- 
sical conductivity is negative resulting in an insulator at 
zero temperature, cr* is an unstable fixed point for the 
symplectic symmetry class. 



19 



However, as we have seen in Sec. ILB.2[ graphene has a 
minimum ballistic conductivity (T,nin = 4e^/7r/i and does 
not become insulating in the limit of vanishing disorder. 
This makes graphene different from the spin-orbit Hamil- 
tonian discussed above, and the question of what happens 
with increasing disorder becomes interesting. 



Bardarson et aL (2007) studied the Dirac Hamiltonian 



(Eq. 2.2 1 with the addition of a Gaussian correlated dis- 
order term U{r), where 



(C/(r)C/(r'))=ifo^^exp 



— r — r 



2^2 



(2.10) 



One should think of Kq as parameterizing the strength 
of the disorder and ^ as its correlation length. If the the- 
ory of one-parameter scaling holds for graphene, then it 
should be possible to rescale the length L* = fo{Ko)L, 
where /o is a scaling function inversely proportional 
to the effective electronic mean free path. Their nu- 
merical results are shown in Fig. [7] and demonstrate 
that (i) graphene does exhibit one-parameter scaling 
(i.e. there exists a /3-function) and (ii) the /3-function 
is always positive unlike the spin-orbit case. Therefore, 
Dirac Fermions evade Anderson localization and are al- 



ways metallic. Similar conclusions were obtained by No- 



mura et aL| (|2007); |San- Jose et a/.| (|2007[); |Titov| (|2007[); 
Tworzydlo et aL||2008[ )." 



The inset of Fig. [7] shows an explicit computation of 
the /3-function comparing Dirac Fermions with the spin- 
orbit model. The difference between these two classes of 
the All symmetry class has been attributed to a topo- 
logical term (i.e. two possible choices for the action of 
the field theory describing these Hamiltonians). Since 
it allows for only two possibilities, it has been called a 
Z2 topological symmetry ( Evers and Mirliii] |2008[ [Kane 



and Mele 



2005b I. The topological term has no effect 



at (J » e /h, but is responsible for the differences at 
a « e'^/h and determines the presence or absence of a 
metal-insulator transition. Nomura et al. (20071 present 
an illustrative visualization of the differences between 
Dirac Fermions and the spin-orbit symplectic class shown 
in Fig. [Sj By imposing a twist boundary condition in the 
wave-functions such that '^{x = 0) = exp[i(/)]^'(a; = L) 
and ^(y = 0) = ^'(y = W)^ one can examine the single 
particle spectrum as a function of the twist angle 0. For 
4> = Q and </> = tt, the phase difference is real and eigen- 
values come in Krammers degenerate pairs. For other 
values of 0, this degeneracy is lifted. As seen in the fig- 
ure, for massless Dirac Fermions all the energy states are 
connected by a continuous variation in the boundary con- 
ditions. This precludes creating a localized state, which 
would require the energy variation with boundary con- 
dition (also called Thouless Energy) be smaller than the 
level spacing. Since this is a topological effect, |Nomura| 
et al. ( 2007 ) argued that this line of reasoning should be 
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FIG. 7 (Color online) Demonstration of one parameter scal- 
ing at the Dirac point ( Bardarson et al\ 20071. Main panel 
shows the conductivity as a function of L'/^, where L* = 
f(Ko)L is the scaled length and ^ is the correlation length 
of the disorder potential. Notice that j3 — din a / d\n L > 
for any disorder strength. The inset shows explicit compar- 
ison (Nomura et al. 20071 of the /3 function for the Dirac 



Fermion model and for the symplectic (All) symmetry class. 



(a) massless Dirac model 



(b) random SO model 





FIG. 8 (Color online) Picture proposed by Nomura et al. 



(20071 to understand the difference in topological structure 



between the massless Dirac model and the random spin-orbit 
symmetry class. 



robust to disorder. 



The situation for the spin-orbit case is very different. 
The same Krammer's pairs that are degenerate at (f> — 0, 
reconnect at (j) — n. In this case, there is nothing to 
prevent localization if the disorder would push the Kram- 
mer's pairs past the mobility edge. Similar considerations 
regarding the 2'2-symmetry hold also for topological in- 
sulators where the metallic surface state should remain 
robust against localization in the presence of disorder. 
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2. Crossover from the symplectic universality class 

As is already apparent in the preceding discussion, 
each Dirac cone is described by the Dirac Hamihonian 
H — a p. The effective SRS H = ay'H*ay is preserved in 
each cone, and for most purposes graphene can be viewed 
as two degenerate copies of the All symplectic symmetry 
class. However, as was first pointed out by |Suzuura and| 
Ando (2002a), a material defect such as a missing atom 



would couple the two Dirac cones (and since each cone is 
located in a different "valley" , this type of interaction is 
called intervalley scattering). One can appreciate intu- 
itively why such scattering is expected to be small. The 
two valleys at points K and K' in the Brillouin zone are 
separated by a large momentum vector that is inversely 
proportional to the spacing between two neighboring car- 
bon atoms. This means that the potential responsible for 
such intervalley coupling would have to vary appreciably 
on the scale of 0.12 nm in order to couple the K and K' 
points. 

We note that while such defects and the corresponding 
coupling between the valleys are commonly observed in 



STM studies on epitaxial graphene (Rutter et al. 2007), 



they are virtually absent in all similar studies in exfoli- 



ated graphene (Ishigami et al.. 2007 Stolyarova et al. 



2007 Zhang et al. 2008 ) . In the presence of such atomi- 



cally sharp disorder, Suzuura and Ando ( 2002a I proposed 



a model for the two- valley Hamiltonian that captures the 
effects of intervalley scattering. The particular form of 
the scattering potential is not important, and below in 
Sec. |II.C3| we will discuss a generalized Hamiltonian that 
includes all non-magnetic (static) impurities consistent 
with the honeycomb symmetry and is characterized by 



five independent parameters (Aleiner and Efetov 2006 



McCann et al. 2006). Here, the purpose is simply to 



emphasize the qualitative difference between two types 
of disorder: long-range (i.e. diagonal) disorder that 
preserves the effective SRS, and a short-range potential 
that breaks this symmetry. 

We note that with the intervalley term U^^, the Hamil- 
tonian belongs to the Wigner-Dyson orthogonal symme- 
try class, while as was discussed in Sec. |H.CT] including 
only diagonal disorder , one is in the Wigner-Dyson 
symplectic class. 

A peculiar feature of this crossover is that it is gov- 
erned by the concentration of short-range impurities thus 
questioning the notion that the "universality class" is de- 
termined only by the global symmetries of the Hamilto- 
nian and not by microscopic details. However, a similar 



crossover was observed in Miller et al. (2003) where the 



strength of spin-orbit interaction was tuned by carrier 
density, moving from weak localization at low density 
and weak spin-orbit interaction, to weak antilocalization 
at high density and strong spin-orbit interaction. 

From symmetry considerations one should expect that 
without atomically sharp defects, graphene would exhibit 
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FIG. 9 (Color online) Diagrammatic representation for (a) 
diffuson, (b) Cooperon, and dressed Hikami boxe s [(c) and 
(d)]. Figure adapted from Kharitonov and Efetov (2008) 



weak antilocalization (where 5a > 0) and no Ander- 
son localization (see Sec. ILC.l). However, with inter- 
valley scattering, graphene should have weak localiza- 
tion {da < 0) and be insulating at zero temperature. 



These conclusions were verified by Suzuura and An"do 
(2002a) from a microscopic Hamiltonian by calculating 



the Cooperon and obtaining the corrections to the con- 



ductivity from the bare Hikami box (see Sec. H.C.S below 
for a more complete discussion). 

For the case of no intervalley scattering U = , the 
resulting Cooperon is 



A 



1 



(vFTQy 



(2.11) 



with area A = LW , Q = k„ k^g and e*^'^''" "''"-''^ « 
— 1, giving (5ctlr = (2e^/7r^S) ln[L0/^]. As expected for 
the symplectic class, without intervalley scattering, the 
quantum correction to the conductivity is positive. 

With intervalley scattering U — calculating the 
same diagrams gives 



1 



(2.12) 



with current ja — —jp and ^ctsr 



This negative Sa is consistent with the orthogonal sym- 
metry class. The explicit microscopic calculation demon- 
strates the crossover from weak antilocalization to weak 
localization induced by atomically sharp microscopic de- 
fects providing intervalley coupling. 

This crossover was recently observed experimen- 



tally (Tikhonenko et al. 20091. They noted empirically 



that for their samples, the scattering associated with 
short-range defects is stronger at high carrier density. In 
fact, this is what one expects from the microscopic theory 



discussed in Sec. [in] below. Due to the unique screening 
properties of graphene, long-range scatterers dominate 
transport at low carrier density while short-range scat- 
terers dominate at high-density. Assuming that these 
short-range defects are also the dominant source of in- 
tervalley scattering, then one would expect to have weak 
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FIG. 10 (Color online) Left panel: Schematic of different 
magnetoresistance regimes: (i) For str ong in tervalley scatter- 
ing (i.e. Ti <^ and <^ r^), Eq. 2.18 gives weak local- 
ization or 5a ~ p{B) — p(0) < 0. This is similar to quan- 
tum transport in the usual 2DEG; (ii) For weak intervalley 
scattering (i.e. Ti ^ and Tz ^ T0), one has weak antilo- 
calization, characteristic of the symplectic symmetry class; 
(iii) For n 3> r^, but Tz <ti t^, Eq. 2.18 gives a regime 
of suppressed weak localization. Right Panel: Experimen- 



tal realization of these three regimes. Graphene magnetocon- 
ductance is shown for carrier density (from bottom to top) 
n = 2.2 X 10^° cm-^ 1.1 x 10^^ cm'^ and 2.3 x 10^^ cm'^ 
at T = 14 K. The lowest carrier density (bottom curve) has a 
small contribution from short-range disorder and shows weak 
anti-localization (i.e. the zero-field conductivity is larger than 
at finite field. a{B — 0) = ao + Sa and a{B > B*) — cto, with 
Sa > 0. B* is the phase-breaking field). In contrast, the 
highest density data (top curve) has a larger contribution of 
intervalley scatteri ng and shows weak localiz ation i.e. Sa < 0. 
Figure taken from Tikhonenko et al. (20091. 



localization at high carrier density (due to the large inter- 
valley scattering), and weak antilocalization at low car- 
rier density where transport is dominated by "atomically 
smooth" defects like charged impurities in the substrate. 
This is precisely what was seen experimentally. Figure 
|10| shows a comparison of the magnetoconductance at 
three different carrier densities. At the lowest carrier 
density, the data show the weak antilocalization charac- 
teristic of the symplectic symmetry class, while at high 
density, one finds weak localization signaling a crossover 
to the orthogonal universality class. 

A second crossover away from the symplectic universal- 



ity class was examined by Morpurgo and Guinea (2006). 



As discussed earlier, a magnetic field breaks time reversal 
symmetry and destroys the leading quantum corrections 
to the conductivity Sa{B > B*) — 0. This can also be 
understood as a crossover from the symplectic (or orthog- 
onal) universality class to the unitary class. The unitary 
class is defined by the absence of time reversal symmetry 
and hence vanishing contribution from the Cooperonj^ 



^ The sign of the quantum correction in relation to the global sym- 
metry of the Hamiltonian can also be obtained from Random 
Matrix Theory ( [Beenakker] |1997| l Sa/ao = (1 - 2/l3)/4, where 
/3 = 1,2,4 for the orthogonal, unitary and symplectic Wigner- 
Dyson symmetry classes. 



Similar to short-range impurities inducing a crossover 
from symplectic to orthogonal classes, [Morpurgo and 
Guinea ( 2006 1 asked if there were other kinds of disor- 



der that could act as pseudo-magnetic fields and induce 
a crossover to the unitary symmetry class leading to the 
experimental signature of a suppression of weak antilo- 
calization. This was in part motivated by the first ex- 
periments on graphene quantum transport showing that 
the weak localization correction was an order of magni- 
tude smaller than expected (Morozov et al. 2006). The 



authors argued that topological lattice defects (Ebbe- 



sen and Takada 1995) (e.g. pentagons and heptagons) 



and non-planarity of graphene (commonly referred to as 
"ripples" ) would generate terms in the Hamiltonian that 
looked like a vector potential and correspond to a pseudo- 
magnetic field. 

In addition, experiments on both suspended 
graphene (Meyer et al. 2007) and on a substrate 
showed that graphene is not perfectly co-planar. It 
is noteworthy, however, that experiments on Si02 
substrate showed that these ripples were correlated with 
the height fluctuations of the substrate and varied by 



less than 1 nm (Ishigami et al. 2007); while graphene 



on mica was even smoother with variations of less than 
0.03 nm (Lui et al. 2009). On the other hand, one could 



deliberately induce lattice defects (Chen et al. 2009c) 



or create controlled ripples by straining graphene before 
cooling and exploiting graphene's negative thermal 



expansion coefficient ( [Balandin et al. 
^0^91). 



2008 Bao et al. 



Just like a real magnetic field, these terms would break 
the TRS in a single valley (while preserving the TRS of 
the combined system). If T[ ^ r^, the two valleys are 
decoupled and these defects would cause a crossover to 
the unitary symmetry class and the resulting Cooperon 
(Fig. [9| would vanish. For example, considering the 
case of lattice defects, the disorder Hamiltonian would 
be given by = (l/4)[(Ta; (g) (t^]V {dyU^{r) - d^Uy{r)), 
where u(r) is the lattice strain vector induced by the de- 
fect. One notices that this term in the Hamiltonian has 
the form of an effective magnetic field +B in the K valley 
and — i? in the K' valley (Morpurgo and Guinea 20061. 



In the absence of intervalley coupling, this would sup- 
press weak antilocalization when the effective magnetic 



field |-B| is larger than the field B* discussed in Sec. II. A 



3. Magnetoresistance and mesoscopic conductance fluctuations 

As already discussed, at low energies and in the ab- 
sence of disorder, graphene is described by two decou- 
pled Dirac cones located at points K and K' in the Bril- 
louin zone. Within each cone, one has a pseudospin space 
corresponding to wavefunction amplitudes on the A and 
B sublattice of the honeycomb lattice. The two valley 
Hamiltonian is then the outer product of two SU(2) spin 
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spaces KK' (E) AB. The most generic Hamiltonian in this 
space of 4 X 4 Hermitian matrices can be parameterized 



by the generators of the group U (4) ( Aleiner and Efetov 



2006 Altland 2006 McCann et al. 2006 ) 



E 

s.l—x,y,z 



E,A,u,z(r) , (2.13) 



where S = (E:r,Sy,I]^) = (a^ (g) ax,a-z (8) cry, I2 (8) a^) 
is the algebra of the sublattice SU{2) space (recall that 
the outer product is in the space KK' ® AB, and the S 
operator is diagonal in the KK' space). Similarly A = 
{a^^^cTz, o-y(8(7z, ^ I2) forms the algebra of the valley- 
spin space (being diagonal in the AB space). 

The Hamiltonian of Eq. 2.13 can be understood in sim- 
ple terms. The first term is just two decoupled Dirac 
cones and is equivalent to the disorder free case discussed 
earlier, but written here in a slightly modified basis. The 
second term is identical to hl^^ and as discussed previ- 
ously represents any long-range diagonal disorder. The 
last term parameterized by the nine scattering potentials 
Usi{r) represent all possible types of disorder allowed by 
the symmetry of the honeycomb lattice. For example, 
a vacancy would contribute to all terms (including uq) 
except Uxz and Uyz; while bond disorder contributes to 



all terms except Uzz (Aleiner and Efetovl 2006) 



The "diagonal" term uo{r) is the dominant scattering 
mechanism for current graphene experiments and origi- 
nates from long ranged Coulomb impurities and is dis- 
cussed in more detail in Sec. |III| Due to the peculiar 
screening properties of graphene, such long-range disor- 
der cannot be treated using the Gaussian white noise ap- 
proximation. To circumvent this problem (both for the 
long-range Uq and short-range Ugi terms) we simply note 
that for each kind of disorder, there would be a corre- 
sponding scattering time {tq, Tsi} that could, in principle, 
have very different dependence on carrier density. 

For the special case of Gaussian white noise, i.e. where 
{usi{r)us'i'{r')) = u'^siSs^s'Si,i'S{r - r'), we have Ht~'^ = 



ttD{Ef)u^i.) Moreover, one could assume that after dis- 
order averaging, the system is isotropic in the xy plane. 
Denoting {x, y} =_L, the total scattering time is given by 



2r 



±z 



2r. 



z± 



-At 



_L_L- 



(2.14) 

These five scattering times could be viewed as indepen- 



dent microscopic parameters entering the theory ( Aleiner 



and Efetov 2006 1, or one could further classify scattering 
times as being either "intervalley" t^^ — AtT] + 2r7^ 

1 _ A^-i I 0,-l 



or "intravalley" ^ = 4tj^^ -|- 2r^^ 



'_L_L ' ^'z± 

A small contri- 



bution from trigonal warping (a distortion to the Dirac 
cone at the energy scale of the inverse lattice spac- 
ing) could be modeled by the perturbative term ~ 
Y,x{'^p)-Az'^x{^p)'^x which acts as an additional source 



of intravalley scattering ( McCann et al. 2006 ) 



The transport properties of the Hamiltonian (Eq. 2.13 ) 
are obtained by calculating the two particle propagator. 
In general, both the classical contribution (diffusons) and 
quantum corrections (Cooperons) will be 4 x 4 matrices 
defined in terms of the retarded (R) and advanced (A) 
Green's functions G^'^ as (see also Fig. ^ 

2?(w, r, r') - (g^(e + lo, r, r') (g> g^{e, r', r)) (2.15) 
C(t^,r,r') = (g^(e + w,r,r')(»e-^(e,r,r')). 



As will be discussed in Sec. HI. A the scattering rate 
is dominated by the diagonal disorder r w tq. Since 
both this term and the Dirac part of Eq. |2.13| is invariant 
under the valley SU{2), one can classify the diffussons 
and Cooperons as "singlets" and "triplets" in the AB- 
sublattice SU{2) space. Moreover, one finds that for 
both the diffussons and Cooperons, only the valley "sin- 
glets" are gapless, and one can completely ignore the 
valley triplets whose energy gap scales as Tq~^. Con- 
sidering only the sublattice singlet {j = 0) and triplet 



(j = x,y,z), one finds ( 


Fal'ko et al. 


2007 Kechedzhi 


et aLl|2008||Kharitonov and Efetov 


2008 


McCann et al. 


2006) 



-iu - 2%'^o ( V 



1 2 2eA 

-lu — -WpTo V H 



2eA 



V^{y,y') = 5{r-v'), 
C^(r,r') = (5(r-r'), 



(2.16) 



with r" = (singlet); = F** = r,- ^ -f ^, and F^ = 2t^ ^ (triplet). This equation captures all the differences in 
the quantum corrections to the conductivity between graphene and usual 2DEGs. 



The magnetoresistance and conductance fluctuations 
properties in graphene follow from this result. The quan- 
tum correction to the conductivity 5a ~ — , 



where Nf is the number of gapless triplet Cooperon 
modes and is the number of gapless singlet Cooperon 
modes. In this context, gapped modes do not have a di- 



23 



vergent quantum correction and can be neglected. Simi- 
larly, the conductance fluctuations are given by {[SG]^) = 
2DEGj where Ncd counts the total number of 
gapless Cooperons and diffusons modes, and {[SG]^)2ueg 
is the conductance variance for a conventional 2D elec- 



tron gas. For a quasi- ID geometry, 
Yg(2e^/ft,)^ (Lee and Ramakrishnan 1985). 



SGY) 



2DEG 



We can immediately identify several interesting 
regimes that are shown schematically in Fig. [TOj 



(i) Strong intervalley scattering. Even with strong 
short-range disorder (i.e. ^ and <C t^), both 
the singlet Cooperon and the singlet diffuson 
remain gapless since r° = 0. Contributions from all 
triplet Cooperons and difussons vanish. In this situation, 
the quantum corrections to the conductivity in graphene 
is very similar to the regular 2DEG. The Hamiltonian 
is in the orthogonal symmetry class discussed earlier 
and one has weak localization (6a < 0). Similarly, for 
conductance fluctuations (typically measured at large 
magnetic fields), one would have the same result as the 
non-relativistic electron gas. 



(a) Weak short-range disorder. For ^ and ^ r^, 
all sublattice Cooperons and diffusons remain gapless 
at zero magnetic field. One then has Sa > or weak 
anti-localization (symplectic symmetry). This regime 



was observed in experiments on epitaxial graphene (Wu 



et al. 2007). The diffuson contribution to the conduc- 
tance fluctuations is enhanced by a factor of 4 compared 
with conventional metals. 



(Hi) Suppressed localization regime. In the case that 
there is strong short-range scattering <^ t^, but weak 
intervalley scattering r^. The Cooperons and 

will be gapped, but will remain and cancel the 
effect of the singlet C°. In this case one would have 
the suppressed weak localization that was presumably 
seen in the first graphene quantum transport experi- 



ments (Morozov et al. 2006). 



Although the discussion above captures the main 
physics, for completeness we reproduce the results of 



calculating the dressed Hikami boxes in Fig. [9] ( Aleiner 


and Efetov, 2006; Kcchedzhi et al. |2008| |Kharitonov and 


Efetov 


2008, McCann et al. 


2006 1 and using known re- 


suits (" 


Lee and Ramakrishnan 


1985). The quantum cor- 



rection to the conductance is 



2e^D f d^q 



(C^ + Cy + C'~ C") , (2.17) 



and for the magnetoresistance 
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where B^^i^^ = {^(^1 ^De)T^j ., and F{x) = lna;-h V'(l/2 + 
1/a;), with ij) as the digamma function. The function 
F(x) is the same as for 2DEGs ^Lee and Ramakrishnan 



1985), however the presence of three terms in Eq. 2.18 
unique to graphene. The universal conductance fluctua- 
tions are 



{[5Gf) = Y.^ 
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(2.19) 



with only diffusions contributing for B > B^ « B* . 

In this section we have assumed Gaussian white noise 
correlations to calculate the Green's functions. Since we 
know that this approximation fails for the semi-classical 
contribution arising from Coulomb disorder, why can we 
use it successfully for the quantum transport? It turns 
out that the quasi-universal nature of weak localization 
and conductance fluctuations means that the exact na- 
ture of the disorder potential will not change the result. 
Several numerical calculations using long-range Coulomb 



potential have checked this assumption (Yan and Ting 



2008 ) . Many of the symmetry arguments discussed here 



apply to confined geometries like quantum dots (Wurm 



et al. 2009 ) . Finally, the diagrammatic perturbation 
theory discussed here applies only away from the Dirac 
point. As discussed in Sec. |II.C.Tj numerically calculated 
weak (anti) localization corrections remain as expected 



even at the Dirac point. However, Rycerz et al. (2007a 
found enhanced conductance fluctuations at the Dirac 
point, a possible consequence of being in the ballistic to 
diffusive crossover regime. 

As for the experimental situation, in addition to the 



observation of suppressed localization (Morozov et al. 



2006) and anti-localization (Wu et al. 



2007), Horsell 



et al. ( 2009 ) made a systematic study of several samples 
fltting the data to Eq. 2.18 to extract t^, and Tz. Their 
data show a mixture of localization, anti-localization and 
saturation behavior. An interesting feature is that the 
intervalley scattering length Li — [DTi)^^"^ is strongly 
correlated to the sample width (i.e. Li « W/2) imply- 
ing that the edges are the dominant source of intervalley 
scattering. This feature has been corroborated by Ra- 
man studies that show a strong D-peak at the edges, but 



not in the bulk (Chen et al. 2009c Graf et al.^ 20071. 



The most important finding of Horsell et al. (2009) is 
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that the contribution from short-range scattering [tz) is 
much larger than one would expect (indeed, comparable 
to tq). The authors showed that any predicted micro- 
scopic mechanisms such as ripples or trigonal warping 
that might contribute to (but not r^) were all negligi- 
ble and could not explain such a large t~^. It remains an 
open question as to why r"^ 3> in the experiments. 



4. Ultraviolet logarithmic corrections 



In addition to these considerations, interaction effects 
could also affect quantum transport (e.g. the Altshuler- 
Aronov phenomena) , particularly in the presence of dis- 
order. Although such interaction effects are probably 
relatively small in monolayer graphene, they may not be 
negligible. Interaction effects may certainly be impor- 
tant in determining graphene transport properties near 



the charge neutral Dirac point ( 


Bistritzer and MacDon- 


aid 20091 IFritz et a?.[|2008||Kashuba 


2008 


Miiller et al. 


2008). 



III. TRANSPORT AT HIGH CARRIER DENSITY 



The semiclassical Boltzmann transport theory treats 
the impurities within the first Born approximation. In 
a diagrammatic perturbation theory, this is the leading 

order term in an expansion of riimp 0. Typically A. Boltzmann transport theory 

for other conventional metals and semiconductors, one 
makes a better approximation by trying to include more 
diagrams that capture multiple scattering off the same 
impurity. For example, in the Self-Consistent Born Ap- 
proximation (SCBA) one replaces bare Green's functions 
with dressed ones to obtain a self-consistent equation for 



the self energy (Bruus and Flensberg, 2004). 



In practice, for graphene, one often finds that attempts 
to go beyond the semiclassical Boltzmann transport the- 
ory described in Sec. |III.A| fare far worse than the sim- 
ple theory. The theoretical underpinnings for the failure 



of SCBA was pointed out by Aleiner and Efetov ( 2006 ) 
where they argued that the SCBA (a standard technique 
for weakly disordered metals and superconductors) , is not 
justified for the Dirac Hamiltonian. They demonstrated 
this by calculating terms to fourth order in perturbation 
theory, showing that SCBA neglects most terms of equal 
order. This could have severe consequences. For exam- 
ple, considering only diagonal disorder, the SCBA breaks 
time reversal symmetry. To further illustrate their point, 



Aleiner and Efetov ( 2006 ) argued that for the full disor- 



der Hamiltonian (Eq. 2.13), considering three impurity 
scattering, there are 54 terms to that order, and only 6 
are captured by the SCBA. 

These terms provide a new divergence in the dia- 
grammatic perturbation series that is distinct from the 
weak localization discussed in Sec. III. C. 31 Unlike weak 
localization that for 2D systems diverges as the size 
{Sa ~ ln[L/£]), this additional divergence occurs at all 
length scales, and was called "ultraviolet logarithmic cor- 
rections". The consequences of this divergence include 
the logarithmic renormalization of the bare disorder pa- 
rameters which was studied using the Renormalization 
Group (RG) in [Foster and Aleiner] ( |2008[ ). For the ex- 
perimentally relevant case of strong diagonal disorder, 
the renormalization does not change the physics. How- 
ever, when all disorder couplings (i.e. intervalley and in- 
travalley) are comparable e.g. relevant for graphene after 
ion irradiation, the system could flow to various strong 
coupling fixed points depending on the symmetry of the 
disorder potential. 



In this section we review graphene transport for large 
carrier densities {n ^ rii, Ui being the impurity density) , 
where the system is homogeneous. We discuss in detail 
the microscopic transport properties at high carrier den- 
sity using the semiclassical Boltzmann transport theory. 

It was predicted that the graphene conductivity lim- 
ited by the short ranged scatterers (i.e., (5-range disor- 
der) is independent of the carrier density, because of the 



linear-in-energy density of states ( Shon and Ando 1998 1 . 
However, the experiments (Fig. 11) show that the con- 
ductivity increases linearly in the carrier density concen- 
tration. To explain this linear-in density dependence of 
experimental conductivity, the long range Coulomb dis- 
order was introduced (|Ando| ]2006| [Cheianov and Fal'ko 
2006a; Hwang et al 



mura and MacDonald^ |^2006 



2007a; Katsncls on et aLj |2009; 'No- 
2007j iJrushin and Schlie^ 



mann 



2008). The long range Coulomb disorder also suc- 



cessfully explains several recent transport experiments. 



Tan et al. ( 2007a ) have found the correlation of the sam- 
ple mobility with the shift of the Dirac point and min- 
imum conductivity plateau width, showing qualitative 
and semi-quantitative agreement with the calculations 
with long range Coulomb disorder (Fig. [TTJb)). [Chen 
et al. (2008a) investigated the effect of Coulomb scat 



terers on graphene conductivity by intentionally adding 
potassium ions to graphene in ultrahigh vacuum, observ- 
ing qualitatively the prediction of the transport theory 
limited by Coulomb disorder (Fig. [TTJc)). Jang et al. 



(2008) tuned graphene's fine structure constant by de- 
positing ice on the top of graphene and observed an en- 
hancement in mobility which is predicted in the Boltz- 
mann theory with Coulomb disorder (Fig. 11 ^d) and 
Chen et al.\ |2009a|b[ and |Kim e^"aH |2009[ ). The role of 
remote impurity scattering was further confirmed in the 
observation of drastic improvement of mobility by reduc- 
ing carrier scattering in suspended graphene through cur- 



rent annealing (Bolotin et al. 2008b Du et al. 2008b I. 



Recent measurement of the ratio of transport scatter- 
ing time to the quantum scattering time by |Hong et al.\ 
(2009b) also strongly supports the long range Coulomb 



25 




FIG. 11 (color online) (a) The measured conductivity a of 
graphene as a function of gate voltage Vg (or carrier density). 
The conductivity increases linearly with the density. Adapted 
from Novoselov et aL ( 2005a I . (b) cr as a function of Vg for five 



different samples For clarity, curves are vertically displaced. 
The inset shows the detailed view of the density-dependent 
conductivity near the Dirac point for the data in the main 
panel. Adapted from Tan et ffl/!| ( 2007a[ ). (c) a versus Vg for 
the pristine sample and three different doping concentrations. 
Adapted from Chen et al. (2008a I. (d) cr as a function Vg for 
pristine graphene (circles) and after deposition of 6 monolay- 
ers of ice (triangles). Inset: Optical microscope image of the 
device. Adapted from Jang et al. (20081. 



disorder as the main scattering mechanism in graphene 
(Fig. [12]) . 

The conductivity a (or mobihty ^ = a/ne) is calcu- 
lated in the presence of randomly distributed Coulomb 
impurity charges with the electron-impurity interaction 
being screened by the 2D electron gas in the random 
phase approximation (RPA). Even though the screened 
Coulomb scattering is the most important scattering 
mechanism, there are additional scattering mechanisms 
(e.g., neutral point defects) unrelated to the charged im- 
purity scattering. The Boltzmann formalism can treat 
both effects, where zero range scatterers are treated with 
an effective point defect density of n^. Phonon scattering 
effects, important at higher temperatures, are treated in 
the next section. We also discuss other scattering mech- 
anisms which could contribute to graphene transport. 

We start by assuming the system to be a homogeneous 
2D carrier system of electrons (or holes) with a carrier 
density n induced by the external gate voltage Vg. When 
the external electric field is weak and the displacement 
of the distribution function from thermal equilibrium is 
small, we may write the distribution function to the low- 
est order in the applied electric field (E) /k = /(ek)+^/k, 
where ek is the carrier energy and /(ek) is the equilibrium 
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FIG. 12 (color online) The ratio of transport scattering time 
(rt) to quantum scattering time (r^) as a function of density 
for different samples. Dashed (solid) lines indicate the theo- 
retical calculations with Coulomb disorder (5-rangc disorder) 
(Hwang and Das Sarma 2008b I. Adapted from |Hong et al. 
( 2009b| ). 



Fermi distribution function and (5/k is proportional to the 
field. When the relaxation time approximation is valid, 
we have (5/k = —^^^^j^eEi-vi^^ij^, where Vk = de^/dk is 
the velocity of carrier and T(ek) is the relaxation time or 
the transport scattering time, and is given by 



1 



2n 



r(ek) 



d^k' 



X |(14,k'W>ni-cosekk']^(ek-ek'), (3-1) 

where ^kk' is the scattering angle between the scattering 
in- and out- wave vectors k and k', n';j^\z) is the con- 
centration of the a-th kind of impurity, and z represents 
the coordinate of normal direction to the 2D plane. In 



Eq. (3.1) (Vk,k'(2:)) is the matrix element of the scat- 
tering potential associated with impurity disorder in the 
system environment. Within Boltzmann transport the- 
ory by averaging over energy we obtain the conductivity 



J deD{e)wlT{e) 



dl 
de 



(3.2) 



and the corresponding temperature dependent resistiv- 
ity is given by p(T) = l/cr(T). Note that /(e^) = 
{1 -|- exp[(efc — p))\/kBT}^^ where the finite temperature 
chemical potential, ^{T), is determined self-consistently 
to conserve the total number of electrons. At T = 0, /(e) 
is a step function at the Fermi energy Ep = ^(T — 0), 
and we then recover the usual conductivity formula: 



'-D{Ef)t{Ef). 



(3.3) 



where vp '^s the carrier velocity at the Fermi energy. 
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B. Impurity scattering 

The matrix element of the scattering potential is de- 
termined by the configuration of the 2D systems and the 
spatial distribution of the impurities. In general, im- 
purities are located in the environment of the 2D sys- 
tems. But for simplicity we consider the impurities are 
distributed completely at random in the plane parallel 
to the 2D systems located at z = d. The location 'd' is 
a single parameter modeling the impurity configuration. 
Then the matrix element of the scattering potential of 
randomly distributed screened impurity charge centers is 
given by 



dzn|'^^(z)|(Vu,k'W)l'=n. 



e{q) 



Fiq) (3.4) 



where q 



Ik - k' 



7kk' 



the number of impu- 



rities per unit area, F{q) is the form factor associated 
with the carrier wave function of the 2D system, and 
Vi{q) = 27re^/(K(7)e~''^ is the Fourier transform of the 
2D Coulomb potential in an effective background lattice 



dielectric constant n. The form factor F{q) in Eq. (3.4) 
comes from the overlap of the wave function. In 2D semi- 
conductor systems it is related to the quasi-2D nature of 
systems, i.e. finite width of the 2D systems. The real 
functional form depends on the details of the quantum 
structures (i.e. heterostructures, square well, etc.). F{q) 
becomes unity in the two dimensional limit (i.e. (5-layer). 
However, in graphene the form factor is related to the 
chirality, not to the quantum structure since graphene is 
strictly a 2D layer. In Eq. (lO), e{q) = e{q,T) is the 



2D finite temperature static RPA dielectric (screening) 
function appropriate and is given by 



:iq,T) = l + v,iqMq,T), 



(3.5) 



where Tl{q, T) is the irreducible finite-temperature polar- 
izability function and Vc{q) is the Coulomb interaction. 
For short-ranged disorder we have 



dzn^t\z)\{y^^M'{z))? =ndVSF{q), 



(3.6) 



where Ud the 2D impurity density and Vq a constant 
short-range (i.e. a (5-function in real space) potential 
strength. 

One can also consider the effect on carrier transport 
by scattering from cluster of correlated charged impuri- 
ties (Katsnelson et al. 20091, as was originally done for 



2D semiconductors by Kawamura and Das Sarma (Das 
Sarma and Kodiyalam[|199"8l|Kawamura and Das Sarma[ 



1996). Without any detailed knowledge about the clus- 
tering correlations, however, this is little more than ar- 
bitrary data fitting. 

Because the screening effect is known to be of vital 



importance for charged impurities (Ando 2006 Hwang 



et al. 2007a), we first provide the static polarizability 



function. It is known that the screening has to be consid- 
ered to explain the density and temperature dependence 



of the conductivity of 2D semiconductor systems (Das 



Sarma and Hwang 1999 20051 and the screening prop- 



erty in graphene exhibits significantly different behavior 



( Hwang and Das Sarma 2007 ) from that in conventional 



2D metals. Also significant temperature dependence of 
the scattering time t may arise from the screening func- 
tion in Eq. (3.5). Thus, before we discuss the details of 
conductivity we first review screening in graphen and in 
2D semiconductor systems. 



1. Screening and polarizability 
(a) graphene: 

The polarizability is given by the bare bubble diagram 



(Ando 


2006 


Hwang and Das Sarma 


2007 


Wunsch et al. 


2006 


) 



n(g,r) 



.9 /sk — /s'k' 



A f-^ £sk — Es'k' 

K.SS 



^^s.'(k,k'). 



(3.7) 



where s — ±1 indicate the conduction (+1) and valence 
(— 1) bands, respectively, k' — k -I- q, £sk — s/iWF|k|, 
F,,-(k,k') = (l + cos0)/2, and /,k = [exp{/3(e,k - m)} + 
l]-i with /3 = l/ksT. 

After performing the summation over ss' it is useful 
to rewrite the polarizability as the sum of intraband and 
interband polarizaibility \l{q, T) = H+(g, T) + H" (g, T). 
At T = 0, the intraband (H+) and interband (H~) polar- 
izability becomes ( Ando 2006[ [Hwang and Das Sarma[ 
2007| 



n+(9) = 



1 - 



1 - 



1 _ 



1 



4kj 
g2 



fi'iq) 



■nq 
8k^' 



q < 2kF 

sin-i q>2kF 
(3.8a) 

(3.8b) 



where H^ — IV^/Do, where Dq = gEp/^Trffvp is the 
DOS at Fermi level. Intraband H+ (interband H^) po- 
larizability decreases (increases) linearly as q increases 
and these two effects exactly cancel out up to g = 2fci?, 
which gives rise to the total static polarizability being 
constant for q < 2kF sls in the 2DEG (Stern 1967), i.e. 
H(g) = n+{q)+Il-{q) = D{Ef) for q < 2kF. In Fig. 
[13] we show the calculated graphene static polarizability 
as a function of wave vector. In the large momentum 
transfer regime, q > 2kF, the static screening increases 
linearly with q due to the interband transition. In a nor- 
mal 2D system the static polarizability falls off rapidly 
for q > 2kF with a cusp at q 



2k F (Stern 1967). The 



linear increase of the static polarizability with q gives 
rise to an enhancement of the effective dielectric constant 
K* iq — > c») = nll+gsgvi^fs/^) in graphene. Note that in 
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a normal 2D system k* — >■ k as g — >■ oo. Thus, the effec- 
tive interaction in 2D graphene decreases at short wave 
lengths due to interband polarization effects. This large 
wave vector screening behavior is typical of an insulator. 
Thus, 2D graphene screening is a combination of "metal- 
lic" screening (due to n"*") and "insulating" screening 
(due to ) , leading, overall, to rather strange screening 
property, all of which can be traced back to the zero-gap 
chiral relativistic nature of graphene. 

It is interesting to note that the non-analytic behav- 
ior of graphene polarizability at q = 2kF occurs in the 
second derivative, (PJl{q)/dq^ oc \/q^ — ikp, i.e., the 
total polarizability, as well as its first derivative are con- 
tinuous at q = 2kF- This leads to an oscillatory decay 
of the screened potential in the real space ( Friedel os- 
dilation) which goes as 4>{r) ~ cos{2kFr)/r^ (Cheianov 



and Fal'ko 2006a Wunsch et al. 20061. This is in con 



trast to the behavior of a 2DEG, where Friedel oscilla- 
tions scale like 4>{r) ~ cos(2fci?r)/r^. The polarizability 
also determines the RKKY interaction between two mag- 
netic impurities as well as the induced spin density due to 
a magnetic impurity, both quantities being proportional 
to the Fourier transform of tt{q)- Like for the screened 
potential, the induced spin density decreases as for 
large distances. Again, this contrasts with the behav- 
ior found in a 2DEG. For the particular case of intrinsic 
graphene the Fourier transform of interband polarizabil- 
ity (n~(g)) diverges [even though n(r) formally scales 
as r~^, its magnitude does not converge], which means 
that intrinsic graphene is susceptible to ferromagnetic or- 
dering in the presence of magnetic impurities due to the 



divergent RKKY coupling (Brey et al. 2007 ) . 

Since the explicit temperature dependence of screen- 
ing gives rise to significant temperature dependence of 
the conductivity, we consider the properties of the polar- 
izability at finite temperatures. The asymptotic form of 
polarizability is given by 



fl((7, T»Tj^) 

n(g,T«rF) 
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(3.9a) 
(3.9b) 



where Tp — Ep/ks is the Fermi temperature. In addi- 
tion, the finite temperature Thomas-Fermi wave vector 
in the g — > long wavelength limit is given by (Ando 



2006 Hwang and Das Sarma 2009b) 



(?,(r>TF) w 81n(2)r,fcF 



T 
Tf 



qsiT <^TF)~Ar,kF 



t:WT 



6 \Tf 



(3.10a) 



(3.10b) 



The screening wave vector increases linearly with tem- 
perature at high temperatures (T ^ Tf), but becomes 
a constant with a small quadratic correction at low tem- 



peratures {T Tf). In Fig. [Ts] we show the finite 
temperature polarizability Il{q,T). 

(b) bilayer graphene: 

For bilayer graphene we have the polarizability of 
Eq. (EtI) with Esk = skV2TO and i^,.s'(k,k') = (1 -I- 



ss' cos26')/2 due to the chirality of bilayer graphene. At 
T = the polarizability of bilayer graphene ( Hwang and 



Das Sarma 2008a) is given by 



niq) - Do [f{q) - 9{q)e{q - 2fcj.)] , 



(3.11) 



where Dq — gsgym/2Tili^ is the BLG density of states at 
the Fermi level and 



1 



5(9) - -v/4 + g4-log 



i + yrTW4" 



;3.12a) 



(3.12b) 



where q — q/kp- 

In Fig. [13] the wave vector dependent BLG polariz- 
ability is shown. For MLG intraband and interband ef- 
fects in polarizability exactly cancel out up to 9 = 2^;^, 
which gives rise to the total static polarizability being 
constant for q < 2kF. However, for BLG the cancella- 
tion of two polarizability functions is not exact because of 
the enhanced backscattering, so the total polarizability 
increases as q approaches 2k f, which means screening 
increases as q increases. Thus BLG, in spite of having 
the same parabolic carrier energy dispersion of 2DEG 
systems, does not have a constant Thomas- Fermi screen- 
ing up to g = 2kF (Borghi et al. 2009 Hwang and Das 



Sarma 2008a) as exists in MLG and 2DEG. In the large 
momentum transfer regime, q > 2kF, the BLG polariz- 
ability approaches a constant value, i.e., H(g) — >• Nq log 4, 
because the interband transition dominates over the in- 
traband contribution in the large wave vector limit. For 
q > 2k F the static polarizability falls off rapidly (~ 1/q^) 
for 2DEG (|Sternl|1967| and for MLG it increases linearly 



with q (see Sec. |III.B.l[ ). 

The long-wavelength (g — >■ 0) Thomas-Fermi screening 
can be expressed as qrF = gsgv'me'^ / kH^ , which is the 
same form as a regular 2D system and independent of 
electron concentration. The screening at q = 2kF is given 
by qs{2kF) = qrFKh^ [^5 - log {(1 -|- \/5)/2}] . Screen- 
ing at g = 2fci;- is about 75 % larger than normal 2D 
TF screening, which indicates that in bilayer graphene 
the scattering by screened Coulomb potential is much 
reduced due to the enhanced screening. 

A qualitative difference between MLG and BLG po- 
larizability functions is at g = 2fci?. Due to the sup- 
pression of 2kF backward scattering in MLG, the total 
polarizability as well as its first derivative are continu- 
ous. In BLG, however, the large angle scattering is en- 
hanced due to chirality [i.e. the overlap factor F^s' in 
Eq. (3.7|], which gives rise to the singular behavior of 
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FIG. 13 (color online) Polarizability n(g, T) in units of the density of states at the Fermi level Da. (a) and (b) show n(g,T) 
of monolayer graphene (a) as a function of wave vector for different temperatures and (b) as a function of temperature for 
different wave vectors, (c) and (d) show BLG polarizability. (e) and (f) show the 2DEG polarizability. 



polarizability at q = 2k p- Even though the BLG polar- 
izability is continuous at g = 2kp, it has a sharp cusp 
and its derivative is discontinuous at 2kp. As q ^ 2kp, 
dli{q)/dq oc 1/ \/q^ — 4:kp. This behavior is exactly the 
same as that of the regular 2DEG, which also has a cusp 
at q = 2kp. The strong cusp in BLG Il{q) a,t q — 2kp 
leads to Friedel oscillations in contrast to the MLG be- 
havior. The leading oscillation term in the screened po- 
tential at large distances can be calculated as 



4:qTpkp sm{2kpr) 



K. {2kf 



CqTFY 



{2kprf 



(3.13) 



where C — -y/S — log[(l -I- -\/5)/2], which is similar to the 
2DEG except for the additional constant C {C ~ \ for 
2DEG), but different from MLG where Friedel oscilla- 
tions scale as 0(r) ~ cos{2kpr)/r^ (Cheianov and Fal'ko 
2006a[ [mmsch aL[|2506| . 



The enhanced singular behavior of the BLG screening 
function at g = 2kp has other interesting consequences 



related to Kohn anomaly ( Kohn 1959 ) and RKKY inter- 



action. For intrinsic BLG the Fourier transform of T]-{q) 
simply becomes a ^-function, which indicates that the lo- 
calized magnetic moments are not correlated by the long 
range interaction and there is no net magnetic moment. 
For extrinsic BLG, the oscillatory term in RKKY interac- 
tion is restored due to the singularity of polarizability at 
q ~ 2kp, and the oscillating behavior dominates at large 
kpr. At large distances 2kpr ^ 1, the dominant oscil- 



lating term in n(r) is given by n(r) cx sm{2kpr) / {kpr)'^ . 
This is the same RKKY interaction as in a regular 2DEG. 

In Fig. [13] the wave vector dependent BLG polariz- 
ability is shown for different temperatures. Note that at 
(7 = 0, n(0,T) — Np for all temperatures. For small q, 
n((7, T) increases as g^. The asymptotic form of polariz- 
ability becomes 



I].{q,T-»Tp) 
U{q,T<$:Tp) 



1 



q 



dkj, T 



(3.14a) 
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More interestingly, the polarizability at g = is temper- 
ature independent, i.e., the finite temperature Thomas- 
Fermi wave vector is constant for all temperatures. 



qs{T) = qxF- 



(3.15) 



In BLG polarizability at g = two temperature effects 
from the intraband and the interband transition exactly 
cancel out, which gives rise to the total static polarizabil- 
ity at (J = being constant for all temperatures. 

(c) 2D semiconductor systems: 

The polarizability of ordinary 2D system was first cal- 
culated by Stern and all details can be found in the liter- 



ature (Ando et al. 1982 Stern 19671. Here we provide 



the 2D polarizability for comparison with graphene. The 
2D polarizability can be calculated with e^k = fi.^k^/2m 
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and Fss' = Sss' /2 because of the non-chiral property of 
the ordinary 2D systems. Il{q) at T = becomes ( Stern 
1967]) 



n(g) = Do 



1-^1- {2kF/qre{q-2kF) 



where Do = gsgvm/2TTh^ is the 2D density of states at 
the Fermi level. Since the polarizability is a constant for 
q < 2kp both the long wave length TF screening and 
2k p screening are same, which is given by — qxF — 
gsgvme'^/Kh^. 

The asymptotic form for the regular 2D polarizability 
are given by 



n((? = 0,T< Ti.) 
n(g,T»Tj.) 



1 - e-^-/^, 



Tf 
T 



1 



6kl T 



(3.17a) 
(3.17b) 



For g = 0, in the T ^ Tp limit, we get the usual Debye 
screening for the regular 2D electron gas system 



qs{T > Tp) w qpp 



Tp 
T 



(3.18) 



A comparison of Eq. (3.181 with Eqs. (3.10) and (3.15) 



shows that the high-temperature Debye screening be- 
haviors are different in all three systems just as the 
low-temperature screening behaviors, i.e., the high tem- 
perature screening wave vector qg in semiconductor 2D 
systems decreases linearly with temperature while q^ in 
MLG increases linearly with temperature and qs in BLG 
is independent of temperature. 

In Fig. [13] we show the corresponding parabolic 2D po- 
larizability normalized by the density of states at Fermi 
level, Dq = gm/lT?2'K. Note that the temperature depen- 
dence of 2D polarizability at g = 2^^- is much stronger 
than that of graphene polarizability. Since in normal 2D 
systems the 2kp scattering event is most important for 
the electrical resistivity, the temperature dependence of 
polarizability at q ^ 2kp completely dominates at low 
temperatures (T <C Tp). It is known that the strong 
temperature dependence of the polarizability function at 
q = 2kp leads to the anomalously strong temperature de- 



pendent resistivity in ordinary 2D systems (Das Sarma 



and Hwang 1999 Stern 1980). 



In the next section the temperature-dependent conduc- 
tivities are provided due to the scattering by screened 
Coulomb impurities using the temperature dependent 
screening properties of this section. 



2. Conductivity 

(a) Single layer graphene: 

The eigenstates of single layer graphene are given by 



the plane wave ■(/'sk(r) 



exp(ik • r)Fsk, where A is 



the area of the system, s = ±1 indicate the conduction 



(-1-1) and valence (—1) bands, respectively, and fJ, 



1 

V2 



sk 



(e*®'', s) with = tan(fcy/fca;) being the polar angle of 
the momentum k. The corresponding energy of graphene 
for 2D wave vector k is given by est = s/iwi? |k|, and the 
(3.16) density of states (DOS) is given by D{e) = g\e\/{2TTh'^Vp), 
where g — gsgv is the total degeneracy {gs — 2, g^ — 2 be- 
ing the spin and valley degeneracies, respectively). The 
corresponding form factor F{q) in the matrix elements of 



Eqs. (3.4) and (3.6) arising from the sublattice symmetry 



(overlap of wave function) (Ando 2006 Auslender and 



Katsnelson 2007) becomes F{q) = (1 + cos9)/2, where 



g = |k — k'l, = 0kk'- The matrix element of the scatter- 
ing potential of randomly distributed screened impurity 
charge centers in graphene is given by 



sk,skVl — 



Vt{q) 



eiq) 



1 + cos 



(3.19) 



and the matrix element of the short-ranged disorder is 



|(V;k.sk')l' = Vo'(l + cos0)/2, 



(3.20) 



where Vq is the strength of the short-ranged disorder 
potential measured in eVm^. The factor (1 — cos 6*) in 



Eq. (3.1) weights the amount of backward scattering of 



the electron by the impurity. In normal parabolic 2D sys- 
tems ( Ando et al. 1982 ) the factor (1 — cos 6) favors large 



angle scattering events. However, in graphene the large 
angle scattering is suppressed due to the wave function 
overlap factor (1 -t- cos 9), which arises from the sublat- 
tice symmetry peculiar to graphene. The energy depen- 
dent scattering time in graphene thus gets weighted by 
an angular contribution factor of (1 — cos0)(l -I- cos 9), 
which suppresses both small-angle scattering and large- 
angle scattering contributions in the scattering rate. 

Assuming random distribution of charged centers with 
density rii, the scattering time r at T = is given by 



(Adam et al. 2007 Hwang and Das Sarma 2008b) 



TT d 

4 — 

To I 2 dr„ 



where Tq 
x'^)f{x) with 



2y/TrniVp/^/n, and g{x) 



■ cosh ^ ^ for x < 1 
= cos^^ - for a; > 1 



(3.21) 

+ (1- 



(3.22) 



Since is independent of the carrier density the scat- 



tering time is simply given by r cx \/n. With Eq. 3.3 
we find the density dependence of graphene conductivity 
a{n) (X n because D{Ep) oc ^/n. For graphene on Si02 
substrate the interaction parameter « 0.8, then the 
conductivity is given by (j{n) « ^i^;^:- (Adam et al. 
2007) On the other hand the corresponding energy de- 



pendent scattering time of short-ranged disorder is 



1 



ridVS Ep 
h A{hvpy 



(3.23) 



30 



100 



-2- 







-500 500 
Density dO'^cm"^) 





10 



n/n. 




FIG. 14 (color online) (a) Calculated graphene conductivity 
as a function of carrier density (rii is an impurity density) 
limited by Coulomb scattering with experimental data. Solid 
lines (from bottom to top) show the minimum conductivity 
of 4e^//i, theory for d = 0, and d = 0.2 nm. The inset shows 
the results in a linear scale assuming that the impurity shifts 
by d = 0.2 nm for positive voltage bias. Adapted from |Hwang| 



e£ al. (2007a I. (b) Graphene conductivity calculated using a 
combination of short and long range disorder. In the calcu- 
lation, Ud/ni = 0, 0.01, 0.02 (top to bottom) are used. In 
inset the graphene mobility as a function of dielectric con- 
stant (k) of substrate is shown for different carrier densities 
n — 0.1, 1, 5x10^^ cm~^ (from top to bottom) in the presence 
of both long ranged charged impurity {ui — 2 x 10^^ cm~^) 
and short-ranged neutral impurity (n^ = 0.4 x lO^'' cm~^). 
Vb = 10 eVnm'^ is used in the calculation. 



Thus, the density dependence of scattering time due to 
the short-range di sorde r scattering is given by Tin) oc 
7T,^i/2, With Eq. 3.3 we find the conductivity to be 
independent of density for short-range scattering, i.e., 
(j{n) oc nP , in contrast to charged impurity scattering 
which produces a conductivity linear in n. 

In Fig. [T4{^a) the calculated graphene conductivity lim- 
ited by screened charged impurities is shown along with 



the experimental data (Chen et al. 2008a Tan et al. 



2007a). In order to get quantitative agreement with ex- 
periment, the screening effect must be included. The ef- 
fect of remote scatterers which are located at a distance 
d from the interface is also shown. The main effect of 
remote impurity scattering is that the conductivity devi- 
ates from the hnear behavior with density and increases 



with both the distance d and n/rii (Hwang et al. 2007a). 

For very high mobility samples, a sub-linear conduc- 
tivity, instead of the linear behavior with density, is 



found in experiments (Chen et al. 2008a Tan et al. 



2007a). Such high quality samples presumably have a 
small charge impurity concentration rii and it is therefore 
likely that short-range disorder plays a more dominant 



role. Fig. 14 b) shows the graphene conductivity calcu- 
lated including both charge impurity and short range dis- 
order for different values of nd/rii. For small Ud/ni the 
conductivity is linear in density, which is seen in most 
experiments, and for large nd/rii the total conductivity 
shows the sub-linear behavior. This high-density flat- 
tening of the graphene conductivity is a non-universal 
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FIG. 15 (color online) Hall mobility as a function of temper- 
ature for different hole densities in (a) monolayer graphene, 
(b) bilayer graphene, and (c) trilayer graphene. The symbols 
are the measured data and the lines are fits. Adapted from 



Zhu et al. (20091. 



crossover behavior arising from the competition between 
two kinds of scatterers. In general this crossover occurs 
when two scattering potentials are equivalent, that is, 
UiV^ ~ ndV^ . In the inset of Fig. 14 |b) the mobility in 



the presence of both charged impurities and short-ranged 
impurities is shown as a function of k. As the scattering 
limited by the short-ranged impurity dominates over that 
by the long-ranged impurity (e.g. ridV^ ^ ^iV?) the 
mobility is no longer linearly dependent on the charged 
impurity and approaches its limiting value 



e [hvpY 
Ah n 



1 



(3.24) 



The limiting mobility depends only on neutral impu- 
rity concentration Ud and carrier density, i.e. long-range 
Coulomb scattering is irrelevant in this high density limit. 

The temperature dependent conductivity of graphene 
arising from screening and the energy averaging defined 
in Eq. (3.2 1 is given at low temperatures (T ^ Tp) 
a{T)/ao w 1 - CiiT/TpY, where ctq = e^vj.D{EF)To/2 
and Ci is a positive constant depending only on 
(Hwang and Das Sarma[ 2009b I. The conductivity de- 



creases quadratically as the temperature increases and 
shows typical metallic temperature dependence. On the 
other hand, at high temperatures {T/Tp ^ 1) it be- 
comes a{T)/ao ~ C2{T/Tp)^, where C2 is a positive 
constant. The temperature dependent conductivity in- 
creases as the temperature increases in the high tempera- 
ture regime, characteristic of an insulating system. (Ts{T) 
of graphene due to the short-range disorder (with scat- 
tering strength Vo) is given by cFs{T) = jq^^^: where 
a,o - e^vj,D{EF)Ts/2 with r, = ^Epvf/{hvp)^. In 
the low temperature limit the temperature dependence 
of conductivity is exponentially suppressed, but the high 
temperature limit of the conductivity approaches aso/2 
as T — 00, i.e., the resistivity at high temperatures in- 
creases up to a factor of two compared with the low tem- 
perature limit resistivity. 

Recently, the temperature dependence of resistivity of 
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graphene has been investigated experimentally (Bolotin 



et al. 


|2008b 


et al. 


2009). 



2008b| pim et al.[ [2008bl |Tan et al] |2007b| |Zhu' 



In Fig. 15 a) the graphene mobility is shown 



ior at high density is observed as explained with screened 
Coulomb impurities. However, it is not obvious whether 
the temperature dependent correction is quadratic be- 
cause phonon scattering also gives rise to a temperature 
dependence (see Sec. IIILCI). 



(b) bilayer graphene: 

The eigenstates of bilayer graphene can be written as 
V'sk — e''""(e~^*^'', s)/V2 and the corresponding energy is 
given by e<;k = sh'^k'^ /2m, where 6]^ = ta.rr^ (ky/kx) and 
s = ±1 denote the band index. The corresponding form 



factor F{q) of Eqs. (3.4) and (3.6) in the matrix elements 
arising from the sublattice symmetry of bilayer graphene 
becomes F{q) = (1 + cos20)/2, where q = |k — k'|, 9 = 
Then the matrix element of the scattering potential 



7kk' 



of randomly distributed screened impurity charge centers 
in graphene is given by (Adam and Das Sarmal, 



2006b 2008) 



2008a 



.Katsnelsonl [20071 IKoshino and Ando||2006j |Nilsson et al. 



|(V^sMk')l' = |«.(g)/£(g)|2(l + cos20)/2. (3.25) 



The matrix element of the short-ranged disorder is 
given by |(14k,sk')P = V^{1 + cos26')/2, and the cor- 
responding energy dependent scattering time becomes 
T~^{ek) = ndV^m/ft' . The density dependent conduc- 
tivity is given by cr(n) ^ -n? in the weak screening limit 
(90 = qxpf^kF 1) or for the unscreened Coulomb 
disorder, and in the strong screening limit (gg ^ 1) 
cr(n) ^ 71. In general for screened Coulomb disorder 
cr(n) ~ 'I 



(Das Sarma et al.. 2010) where a is den- 



sity dependent and varies slowly changing from 1 at low 
density to 2 at high density. Increasing temperature, 
in general, suppresses screening, leading to a slight en- 
hancement of the exponent a. For short range disorder 
cr(n) ~ n. 
Fig. 



Fig. 17 



16 shows the experiment of BLG conductivity. In 
'a) the density dependent conductivities both for 
screened Coulomb disorder and for short range disorder 
are shown. For screened Coulomb disorder the conduc- 
tivity shows super-linear behavior, which indicates that 
pure Coulomb disorder which dominates mostly in MLG 
transport can not explain the density dependent con- 



ductivity as seen experimentally (see Fig. 16) (Morozov 



et al. 2008 Xiao et al. 2010). The density dependence 



of conductivity with both disorders is approximately lin- 
ear over a wide density range, which indicates that BLG 
carrier transport is controlled by two distinct and inde- 
pendent physical scattering mechanisms, i.e. screened 
Coulomb disorder due to random charged impurities in 
the environment and a short-range disorder. The weaker 
scattering rate of screened Coulomb disorder for BLG 
than for MLG is induced by the stronger BLG screening 
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FIG. 16 (color online) The measured conductivity of bilayer 
graphene as a function of gate voltage Vg (or carrier density) . 
The measured conductivity increases linearly with the density. 
Adapted fromlMorozov et aZ.|(|2008|. 
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FIG. 17 (color online) (a) Density dependence of bilayer 
graphene conductivity with two scattering sources: screened 
long-range Coulomb disorder and short-ranged neutral disor- 
der, (b) Density dependence of BLG conductivity for differ- 
ent temperatures, T=0 K, 50 K, 100 K, 150 K, 200 K, 300 K 
(from bottom to top). Top inset shows cr as a function of T 
in presence of short range disorder. Bottom inset shows a as 
a function of T in presence of screened Coulomb disorder for 
different densities n = [5, 10, 30] X lO^'^ cm"'^ (from bottom 



to top). Adapted from Das Sarma et al. (20101 



than MLG screening, rendering the effect of Coulomb 
scattering relatively less important in BLG (compared 
with MLG). 

The temperature dependent conductivity due to 



screened Coulomb disorder (Adam and Stiles 



2010 



Das Sarma et al. |2010[ [Hwang and Das Sarma[ 
Lv and Wan[ |2010[) is given by a(T)/ao « 1 - Cq{T/Tf) 



2010 



at low temperatures, where Cq = 41og2/(C-|- 1/qo) with 
<Zo = QTr/^kp, and a{T) w ai{T/Tp)'^ at high tempera- 
tures. When the dimensionless temperature is very small 
[T/Tp ^ 1) a linear-in T metallic T dependence arise 
from the temperature dependence of the screened charge 
impurity scattering, i.e. the thermal suppression of the 



2fci;--peak associated with back-scattering (see Fig. 13). 
For the short-ranged scattering the temperature depen- 
dence only comes from the energy averaging and the con- 
ductivity becomes <j{T) = cr(0)[l-|-iln(H-e~^/*)], where 
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t = T/Tp. At low temperatures the conductivity is ex- 
ponentially suppressed, but at high temperatures it in- 
creases linearly. 

Fig. [iTt^b) shows the finite temperature BLG conduc- 
tivity as a function of n. The temperature dependence is 
very weak at higher densities as observed in recent exper- 
iments (Morozov et al. 2008). At low densities, where 



Sarachik 2004 1 . More recent perspectives can be found 



Das Sarma and Hwang (2005); Spivak et al. (2010). 



T/Tp is not too small, there is a strong insulating-type 
T dependence arising from the thermal excitation of car- 
riers (which is exponentially suppressed at higher den- 
sities) and energy averaging, as observed experimentally 
dMorozov e^"aH [20081 ). Note that for BLG Tp = 4.23n K, 
^). In bottom inset the conductiv- 



where ; 



n/(10 



tit: 



ity due to screened Coulomb disorder is shown as a func- 
tion of temperature for different densities. At low tem- 
peratures [T/Tp <^ 1) the conductivity decreases linearly 
with temperature, but <j{T) increases quadratically in 
high temperature limit. By contrast, for the short-range 
disorder a always increases with T, as shown in the upper 
inset of Fig. [T7|^b) . Thus for bilayer graphene the metal- 
lic behavior due to screening effects is expected at very 
low temperatures for low mobility samples, in which the 
screened Coulomb disorder dominates. In Fig. [Tsj^b) the 
temperature dependence of mobility for bilayer graphene 
is shown. As we expect the metallic behavior shows up 
at very low temperatures (T < 100 K). 

We conclude this section by emphasizing the similarity 
and the difference between BLG and MLG transport at 
high densities from the perspective of Boltzmann trans- 
port theory considerations. In the MLG the linear den- 
sity dependent conductivity arises entirely from Coulomb 
disorder. However, in the BLG the existence of short- 
range disorder scattering must be included to explain 
the linearity because the Coulomb disorder gives rise to 
a higher power density dependence in conductivity. The 
importance of short-range scattering in BLG compared 
with MLG is understandable based on BLG screening 
being much stronger than MLG screening leading to the 
relative importance of short-range scattering in BLG. 

(c) 2D semiconductor systems: 

Transport properties of 2D semiconductor based 
parabolic 2D systems (e.g. Si MOSFETs, GaAs het- 
erostructures and quantum wells, SiGe-based 2D struc- 
tures) have been studied extensively over the last 



forty years ( 


Abrahams et al. 


2001 


Kravchenko and Sarachik 




2004 


). 



Ando et al. 1982 



More recently, 2D 
transport properties have attracted much attention be- 
cause of the experimental observation of an appar- 
ent metallic behavior in high-mobility low-density elec- 
tron inversion layer in Si metal-oxide-semiconductor- 



field-effect transistor (MOSFET) structures (Kravchenko 
et al. 1994). However, in this review we do not make 
any attempt at reviewing the whole 2D MIT literature. 
Early comprehensive reviews of 2D MIT can be found in 



Our goal in this review is to provide a direct comparison 
of the transport properties of 2D semiconductor systems 
with those of MLG and BLG, emphasizing similarities 
and differences. 

It is well known that the long-range charged impurity 
scattering and the short-range surface-roughness scatter- 
ing dominate, respectively, in the low and the high car- 
rier density regimes of transport in 2D semiconductor 
systems. In Fig. [18] the experimental mobility of Si- 
MOSFETs is shown as a function of density. As den- 
sity increases, the measured mobility first increases at 
low densities and after reaching the maximum mobility 
it decreases at high densities. This behavior is typical for 
all 2D semiconductor systems, even though the mobility 
of GaAs systems decreases very slowly at high densities. 
This mobility behavior in density can be explained with 



mainly two scattering mechanisms as shown in Fig. 18 'b). 
In the low temperature region phonons do not play much 
of a role in resistive scattering. At low carrier densities 
long-range Coulomb scattering by unintentional random 
charged impurities invariably present in the environment 
of 2D semiconductor systems dominates the 2D mobility 
(Ando et al. 1982). However, at high densities as more 



carriers are pushed to the interface the surface roughness 
scattering becomes more significant. Thus transport in 
2D semiconductor systems is limited by the same mech- 
anisms as in graphene even though at high densities the 
unknown short-range disorder in graphene is replaced by 
the surface roughness scattering in 2D semiconductor sys- 
tems. The crucial difference between 2D transport and 
graphene transport is the existence of the insulating be- 
havior of 2D semiconductor systems at very low densities 
which arises from the gapped nature of 2D semiconduc- 
tors. However, the high density 2D semiconductor trans- 
port is not qualitatively different from graphene trans- 
port since charged impurity scattering dominates carrier 
transport in both cases. 

The experimentally measured conductivity and mo- 
bility for three different systems as a function of den- 
sity are shown in Figs. [T8][l9l At high densities, the 
conductivity depends on the density approximately as 
a (X n" with 1 < a < 2, where a{n) depends weakly 
on the density for a given system but varies strongly 
from one system (e.g., Si-MOSFET) to another (e.g. 
GaAs). At high densities, before surface roughness scat- 
tering sets in the conductivity is consistent with screened 
charged impurity scattering for all three systems. As n 
decreases, a{n) starts decreasing faster with decreasing 
density and the experimental conductivity exponent a 
becomes strongly density dependent with its value in- 
creasing substantially, and the conductivity vanishes as 
the density further decreases. To explain this behav- 
ior a density-inhomogeneity-driven percolation transition 



the literature (Abrahams et al. 2001 Kravchenko and was proposed (Das Sarma et al. 2005), i.e. the density- 
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FIG. 18 (Color online) (a) Experimental mobility /i as a 
function of density for two Si-MOSFET samples at a temper- 
ature T=0.25 K. (b) Calculated mobility with two different 
scatterings, i.e. charge impurities and surface-roughness scat- 
terings, (c) Measured conductivity (j(n) for Si-MOSFET as a 
function of electron density n for two different samples. The 
solid lines are fits to the data of the form (j{n) oc A{n — UpY . 
The upper and lower insets show the exponent p and critical 
density Up, respectively, as a function of temperature. Solid 
lines a re a guides for the eyes. Adapted from [Tracy et al.\ 
(20091). 



dependent conductivity vanishes as cr{n) cx (n — np)P 
with the exponent p = 1.2 being consistent with a perco- 
lation transition. At the lowest density, linear screening 
in a homogeneous electron gas fails qualitatively in ex- 
plaining the cr(ri) behavior whereas it gives quantitatively 
accurate results at high densities. As has been found 
from direct numerical simulations (Efros 1988 Nixon 



and Davies 1990 Shi and Xie 2002 ) homogeneous linear 



screening of charged impurities breaks down at low car- 
rier densities with the 2D system developing strong inho- 
mogeneities leading to a percolation transition at n < rip. 
Nonlinear screening dominates transport in this inhomo- 
geneous low carrier density regime. For n < Hp, the 
system is an insulator containing isolated puddles of elec- 
trons with no metallic conducting path spanning through 
the whole system. By contrast, graphene, being gapless, 
goes from being an electron metal to a hole metal, i.e., 
the conductivity is always finite for all densities, as the 
chemical potential passes through the puddle region. 

Except for being an insulator at very low densities the 
transport behavior of 2D semiconductor systems is not 
qualitatively different from graphene transport because 
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FIG. 19 (a),(b) Experimentally measured (symbols) and cal- 
culated (lines) conductivity of two different n-GaAs samples. 
The high density conductivity limited by the charged impu- 
rities fit well to the experimental data. Adapted from |Das| 



Sarma et al. (20051. (c) Mobility of p-GaAs 2D system vs. 
density at fixed temperature T =47 mK. (d) The correspond- 
ing conductivity vs. density (solid squares) along with the fit 
generated assuming a percolation transition. The das hed line 
in (c) indicates the /.i ~ p"'^ behavior. Adapted from Manfra 
et al.\\2m7\. 



both systems are governed by the charged impurities. 
To understand the p (or a) behavior at high density we 
start with the Drude-Boltzmann semiclassical formula, 
Eq. (3.2), for 2D transport limited by screened charged 



impurity scattering (Das Sarma and Hwang 1999 ). How- 



ever, due to the finite extent in the z-direction of the real 
2D semiconductor system the Coulomb potential has a 
form factor depending on the details of the 2D structure. 
For comparison with graphene we consider the simplest 
case of 2D limit, i.e. ^-layer. For (5-layer 2D systems with 
parabolic band the scattering times at T = for charged 
impurity centers with impurity density rii located at the 
2D systems is calculated by 



aqo 



To 



(3.26) 



where ^ = 27r/i^( 



- -""mV^)^9o. <?0 = qTEf^kp {qTF IS a 
2D Thomas-Fermi wave vector), and f{x) is given in Eq. 
(3.22). Then, the density dependence of conductivity 
can be expressed as <T(n) oc n" with 1 < a < 2. In 
the strong screening limit (go ^ 1) the scattering time 
becomes oc (/g oc n~^, then the conductivity behaves 
as (j{n) oc n^. In the weak screening limit rx rx nP 
and cr{n) oc n. These conductivity behaviors are common 
for 2D systems with parabolic bands and are qualitatively 
similar to graphene where a oc n behavior is observed. 
However, due to the complicated impurity configuration 
(spatial distribution of impurity centers) and finite width 
effects of real 2D semiconductor systems the exponent 
a varies with systems. In general, modulation doped 
GaAs systems have larger a than Si-MOSFETs due to 
the configuration of impurity centers. 
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temperature behaviors of 2D conductivity are given by: 
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FIG. 20 (Color online) (a) Experimental resistivity p of Si- 
MOSFET as a function of temperature at 2D electron den- 
sities (from top to bottom) n = [1.07, 1.10, 1.13, 1.20, 1.26, 
1.32, 1.38, 1.44, 1.50, 1.56, 1.62, and 1.68] x lO" cm^ Inset 
(b) shows p forn = [1.56, 1.62, and 1.68] x 10"cm^ (c) The- 
oretically calculated temperature and density-dependent re- 
sistivity for sample A for densities n — [1.26, 1.32, 1.38, 1.44, 
1.50, 1.56, 1.6 2, and 1.68] x lO" cm^ (from top to bottom). 

( |2009| ). (e) Experimental p{T) for 
^). The density ranges 
"cm~^. 



Adapted from Tracy et 



2.3 X 10^ 



cm 



n-GaAs (where n, 

from 16xl 0"cm~'' to 1.06x 10^"cm"^. Adapted from [Lilly 
\et al. \ ( 2003| ). (f) Temperature dependence of the resistivity 
for p-GaAs systems for densities ranginR from 9.0 x 10 ^ cm~^ 
to 2.9x10^ cm-2. Adapted from iManfra et aZ. I (12007 1. 



An interesting transport property of 2D Semiconductor 
systems is the remarkable observation of the extremely 
strong anomalous metallic (i.e., dp/dT > 0) tempera- 
ture dependence of the resistivity p{T) in the density 
range just above a critical carrier density Uc where dp/dT 



changes its sign at low temperatures (see Fig. 201, which 
is not seen in graphene. Note that the experimentally 
measured p{T) of graphene shows very weak metallic 
behavior at high density due to the weak temperature 



dependence of screening function. It is suggested (Das 



Sarma and Hwang 1999) that the anomalously strong 
metallic temperature dependence discovered in 2D semi- 
conductor systems arises from the physical mechanism of 
temperature, density, and wave vector dependent screen- 
ing of charged impurity scattering in 2D semiconduc- 
tor structures, leading to a strongly temperature depen- 
dent effective quenched disorder controlling p{T, n) at low 
temperatures and densities. Interaction effects also lead 



to a linear-T conductivity in 2D semiconductors (Zala 



et al. 2001 1 



With temperature dependent screening function 



e{q,T) in Eq. (3.51, the asymptotic low (Das Sarma and 



Hwang[ 20031 and high (Das Sarma and Hwang 2004) 



cr(t < 1) 
a{t > 1) 



•2D 



-2D 



1 - Ci (T/Tf)] , (3.27a) 
T/Tf + (3V^go/4)yi>7r]3-27b) 



where t = T/Tf, a'o = <T = 0), Ci = 2go/(l + go), 
and cr\^ — {e^ /h){n/ni)TTqQ. Here an ideal 2D elec- 
tron gas with zero thickness is considered in order to 
compare with the 2D graphene sheet which also has a 
zero thickness. It is important to include the tempera- 
ture dependent polarizability of Fig. [13] in the calcula- 
tion in order to get strong temperature dependent re- 
sistivity. Since the most dominant scattering occurs at 
q = 2kF and the temperature dependence of screening 
function at 2k f is strong, the calculated 2D resistivity 
shows the strong anomalous linear T metallic behav- 
ior, which is observed in many different semiconductor 
systems (e.g. 



Si-MOSFET (Kravchenko et al.\ 19941, p- 



GaAs (Manfra et al. 2007 Noh et al. 20031, n-GaAs 



(ILilly et aZ.l|2003|, SiGe (Senz et aLl|2002p, AlAs (|Pa 



padakis and Shayegan 1998)). In addition, for the ob- 



servation of a large temperature-induced change in resis- 
tivity it is required to have a comparatively large change 
in the value of the dimensionless temperature t = T/Tp 
and the strong screening condition, qxF ^ 2fci?, which 
explains why the Si MOS 2D electron system exhibits 
substantially stronger metallic behavior than the GaAs 
2D electron system, as is experimentally observed, since 
(gTF/2fc_F)sj ~ W{qTF/'2kF)GaAs at similar density. 

Before concluding this basic transport theory section 
of this review we point out the key qualitative similarities 
and differences in the transport theory of all systems (i.e. 
graphene, bilayer graphene, and 2D semiconductor based 
parabohc 2D systems). First, the graphene conductivity 
is qualitatively similar to that of 2D semiconductor sys- 
tems in the sense that the conductivity at high density 
of both systems follows the power law in terms of den- 
sity, cr{n,) ~ n" . The formal Boltzmann theory for the 
scattering times is the same in all systems except for the 
different angular factor arising from chiral properties of 
graphene. This angular part does not play a role in the 
density dependence of conductivity, but significantly af- 
fects the temperature dependence of conductivity. The 
explicit differences in the density of states D{e) and the 
dielectric function e{q, T) also lead to different temper- 
ature dependent conductivities in these systems. The 
most important qualitative difference between graphene 
and semiconductor 2D systems occurs at low carrier den- 
sities, in which semiconductor 2D systems become insula- 
tors, but graphene conductivity is finite for all densities. 



C. Phonon Scattering in Graphene 

In this section we review the phonon scattering lim- 
ited carrier transport in graphene. Lattice vibrations are 
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inevitable sources of scattering and can dominate trans- 
port near room temperature. It is an intrinsic scatter- 
ing source of the system, i.e., it limits mobility at finite 
temperatures when all extrinsic scattering sources are re- 
moved. In general, three different types of phonon scat- 
tering are considered, i.e., intravalley acoustic (optical) 
phonon scattering which induce the electronic transition 
within a single valley via a acoustic (optical) phonons, 
and intervalley scattering, i.e. electronic transition be- 
tween different valleys. 

The intravalley acoustic phonon scattering is induced 
by low energy phonons and is considered an elastic pro- 
cess. The temperature dependent phonon-limited resis- 
tivity (Hwang and Das Sarma 2008a| [Stauber et al. 



2007) was found to be hnear (i.e. pph oc T) for T > Tbg 
where Tbq is the Bloch-Griineisen (BG) temperature 
(Kawamura and Das Sarma 19921, and pph{T) ^ 
for T < Tbg- The acoustic phonon scattering gives 
a quantitatively small contribution in graphene even at 
room temperature due to the high Fermi temperature of 
graphene in contrast to 2D semiconductors where room- 
temperature transport is dominated by phonon scatter- 
ing ( Kawamura and Das Sarma[ 1990 1992). The in- 
travalley optical phonon scattering is induced by optical 
phonons of low momentum (q « 0) and very high en- 
ergy {ujop ~ 200 meV in graphene ) and is negligible. 
The intervalley scattering can be induced by the emission 
and absorption of high momentum, high energy acoustic 
or optical phonons. In graphene intervalley scattering 
may be important at high temperatures because of rel- 
atively low phonon energy (« 70 meV, the out-of-plane 



acoustic (ZA) phonon mode at the K point) (Maultzsch 



et al. 2004 Mounet and Marzari 20051. Even though 



the effects of inter-valley phonon scattering can explain 
a crossover (Figs. 21 and 22) seen in experiments in 



the 150 K to 250 K range, more work is needed to val- 
idate the model of combined acoustic phonon and ZA 
phonon scattering contributing to the temperature de- 
pendent graphene resistivity. 

The remote interface polar optical phonons in the sub- 



strate (i.e. Si02) have recently been considered (Chen 
et al] |2008b| [Fratini and Guinea[ |2008[ ). Even though 



Before we discuss the phonon transport theory we men- 
tion the phonon contribution obtained from the experi- 
mental resistivity data. Since the experimentally mea- 
sured resistivity in the current graphene samples is com- 
pletely dominated by extrinsic scattering (impurity scat- 
tering described in Sec. III.B ) even at room temper- 



atures the experimental extraction of the pure phonon 
contribution to graphene resistivity is not unique. In 
particular, the impurity contribution to resistivity also 
has a temperature dependence arising from Fermi statis- 
tics and screening which, although weak, cannot be ne- 
glected in extracting the phonon contribution (particu- 
larly since the total phonon contribution itself is much 
smaller than the total extrinsic contribution). In addi- 
tion, the experimental phonon contribution is obtained 
assuming Matthiessen's rule, i.e. ptot = Pph + Pi where 
ptot is the total resistivity contributed by impurities and 
defects (pi) and phonons {ppt), which is not valid at room 
temperature (Hwang and Das Sarma 2008a). Thus, two 



different groups ( Chen et al] 2008b Morozov et al. 2008 ) 



(Figs. 21 and 22 ) have obtained totally different behavior 



of phonon contribution of resistivity. In [Morozov et al] 
(2008) it is found that the temperature dependence is a 



rather high power (T^) at room temperatures, and the 
phonon contribution is independent of carrier density. In 
Chen et al. ( 2008b[ ) the extracted phonon contribution is 
strongly density dependent and is fitted with both linear 
T from acoustic phonons and Bose-Einstein distribution. 
Therefore the phonon contribution, as determined by a 
simple subtraction, could have large errors due to the 
dominance of extrinsic scattering. 

In this section we describe transport only due to the 
longitudinal acoustic (LA) phonons since either the cou- 
pling to other graphene lattice modes is too weak or the 
energy scales of these (optical) phonon modes are far too 
high for them to provide an effective scattering channel 
in the temperature range (5 K to 500 K) of our inter- 
est. Since graphene is a non-polar material the most im- 
portant scattering arises from the deformation potential 
due to quasi-static deformation of the lattice. Within the 



Boltzmann transport theory ( Kawamura and Das Sarma 



1990 1992) the relaxation time due to deformation po- 



these modes are known to be not very important in Si- tential coupled acoustic phonon mode is given by 



MOSFETs (Hess and Vogl 1979 Moore and Ferry 1980) 



their role in graphene transport seems to be important 



(Chen et al. 2008b DaSilva et al. 2010). Another pos- 



sibility considered in Morozov et al. (2006) is that the 



1 



^(l-cos0kk')W^kk'Y 



fie') 



lie) 



(3.28) 



thermal fluctuations (ripplons) of the mechanical rip- 
ples invariably present in graphene samples contribute 
to the graphene resistivity. In addition [Mariani and von| 
Oppen ([2008) investigated the role of the flexural (out- 
of-plane) phonons of free standing graphene membranes 
which arise from the rotation and reflection symmetries. 
Flexural phonons make a contribution to the resistivity 
at low temperatures with an anomalous temperature de- 
pendence oc T^/^ InT. 



where 0kk' is the scattering angle between k and k', e = 
/lUi? |k|, and Wkk' is the transition probability from the 
state with momentum k to the state with momentum k' 
and is given by 

^kk' = ^El^(q)l'^(^'^') (3.29) 



where C(q) is the matrix element for scattering by acous- 
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FIG. 22 (color online) Temperature dependent resistivity 
for four different MLG samples (symbols). The solid curve 
is the best fit by using a combination of T and functions. 
The inset shows T dependence of maximum resistivity at the 
neutrality point for MLG and BLG (circles and squares, re- 
spectively). Adapted from Morozov et al. (20081. 



FIG. 21 (color online) Temperature-dependent resistivity of 
graphene on Si02. Resistivity of two graphene samples as a 
function of temperature for different gate voltages. Dashed 
lines are fits to the linear T-dependence with Eq. ( |3.32[ ). 



Adapted from |Chen et aL| ( |2008bp 



tic phonons, and A(e,e') is given by 

A{e,e') = Ngd{e-e'+uj^) + {Ng + l)5{e-e'~ujq), (3.30) 

where Wq = hvpiiq is the acoustic phonon energy with 
being the phonon velocity and Nq the phonon occupation 
number Nq = l/(exp(/3a;q) — 1). The first (second) term 
is Eq. (3.30) corresponds to the absorption (emission) 
of an acoustic phonon of wave vector q = k — k'. The 
matrix element C(q) is independent of the phonon oc- 
cupation numbers. The matrix element |C(q)p for the 
deformation potential is given by 



D^hq 

2ApraVph 



1 



(3.31) 



where D is the deformation potential coupling constant, 
Pm is the graphene mass density, and A is the area of the 
sample. 

The scattering of electrons by acoustic phonons may be 
considered quasi-elastic since hu^ <^ Ep , where Ep is the 
Fermi energy. There are two transport regimes, which ap- 
ply to the temperature regimes T <C Tbg and T ^ Tbg, 
depending on whether the phonon system is degener- 
ate (Bloch-Griineisen, BG) or non-degenerate (equiparti- 
tion, EP). The characteristic temperature Tbg is defined 
as ksTBG = 'i.hkpVph, which is given, in graphene, by 
Tbg = '2vphkp /kB ~ bA^/n K with density measured 



in unit of n = 10 



12^ 



The relaxation time in the EP 



regime is calculated to be ( 


Hwang and Das Sarma 


2008a 


Stauber et al. 


2007 Vasko and Ryzhii 


2007 


) 



Z?2 



-(e) 



'Jp PrnVph 



-kuT. 



(3.32) 



Thus, in the non-degenerate EP regime (fiwq ^ ksT) the 
scattering rate [1/T(e)] depends linearly on the temper- 
ature. At low temperatures {Tbg ^ T <^ Ep/kB) the 
calculated conductivity is independent of electron den- 
sity. Therefore the electronic mobility in graphene is in- 
versely proportional to the carrier density, i.e. jj, ex. 1/n. 
The linear temperature dependence of the scattering time 



has been reported for nanotubes (Kane et al. 1998) and 



graphites (Pietronero et al. 1980 Suzuura and Ando 



2002b Woods and Mahan 2000). 



In BG regime the scattering rate is strongly reduced by 
the thermal occupation factors because the phonon pop- 
ulation decreases exponentially, and the phonon emission 
is prohibited by the sharp Fermi distribution. Then, in 
the low temperature limit T <^ Tbg the scattering time 
becomes ( Hwang and Das Sarma[ [2'008a ) 



1 



111 4!C(4) 

TT Ep kp 2pmVph (flVph) 



,{kBT)\ 



(3.33) 



Thus, the temperature dependent resistivity in BG 
regime becomes p oc T^. Even though the resistivity 
in EP regime is density independent, Eq. (3.33) indi- 



cates that the calculated resistivity in BG regime is in- 
versely proportional to the density, i.e. pbg oc n~^^'^ 
since p cx [D{Ep){t)]~^ . More experimental and theo- 
retical work would be needed for a precise quantitative 
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FIG. 23 (color online) (a) Acoustic phonon-limited mobility 
of n-GaAs 2D system as a function of density for two different 
temperatures, (b) Calculated n-GaAs mobility as a function 
of temperature for different impurity densities. At low tem- 
peratures (T < lA") the mobility is completely limited by 
impurity scattering. Adapted from |Hwang and Das Sarma| 
(2008b|). 



understanding of phonon scattering effect on graphene 
resistivity. 



D. Intrinsic mobility 

Based on the results of previous sections one can ex- 
tract the possible (hypothetical) intrinsic mobility of 2D 
systems when all extrinsic impurities are removed. In 



Fig. 23 the acoustic phonon-limited mobility is shown 
for 2D n-GaAs system. For lower temperatures, /i(T) 
increases by a large factor (/i cx T^^ for deformation- 
potential scattering and \x oc for piezoelectric scat- 
tering) since one is in the Bloch-Griineisen regime where 



phonon occupancy is suppressed exponentially (Kawa- 



mura and Das Sarma 1990 1992 1. Thus the intrinsic 



mobility of semiconductor systems is extremely high at 
low temperatures (T < Tbg)- For currently available 
semiconductor samples the mobility below Tbg is com- 
pletely limited by extrinsic impurity scattering in 2D sys- 
tems. Above the BG regime (or T > 4 K) the mobility 
is dominated by phonons. In this limit the mobility lim- 
ited by phonon scattering is much lower than that for 
charged impurity scattering. Therefore it will be impos- 
sible to raise 2D mobility (for T > 4 K) by removing 
the extrinsic impurities since acoustic phonon scattering 
sets the intrinsic limit at these higher temperatures (for 



T > 100 K, optical phonons become dominant) (Pfeiffer 



et all 1989 



In Fig. |24[ the acoustic phonon-limited graphene mo- 
bility, n = (enp)~^, is shown as functions of temperature 
and carrier density, which is given by /x > lO^'^ / D^riT 
cm^ /Vs where D is measured in eV, the temperature T in 
K, and n carrier density measured in units of 10^^ cm~^. 
Thus, the acoustic phonon scattering limited graphene 
mobility is inversely proportion to T and n for T > Tbg- 
Also with the generally accepted values in the literature 
for the graphene sound velocity and deformation coupling 




10 



o 



' 50 100 150 200 250 300 
T(K) 



10^ 



T=300 K 



1 2 3 4 
n(10''cm-2) 



FIG. 24 (color online) Calculated graphene mobility limited 
by the acoustic phonon with the deformation potential cou- 
pling constant D — 19 eV (a) as a function of temperature 
and (b) as a function of density. Adapted from |Hwang and| 
Das Sarma|(|2008a|). 



( [Chen et al] |2008b| (i.e., Vph = 2 x 10^ cm/s, and de- 
formation potential D = 19 eV) fi could reach values as 
high as 10*^ cm^/Vs for lower carrier densities (n < 10^^ 
cm 



) at T = 300i^ (Hwang and Das Sarma 



2008a 



Shishir and Ferry 2009). For larger (smaller) values of 
D, fi would be smaller (larger) by a factor of D^. It 
may be important to emphasize here that we know of no 
other system where the intrinsic room-temperature car- 
rier mobility could reach a value as high as 10^ cm^/Vs, 
which is also consistent with the experimental conclusion 



by( 


Chen et al. 


2008b 


Hong et al. 


2009a 


Morozov et al.\ 


2008 


\. This would, however, require the elimination of 



all extrinsic scattering, and first steps in this direction 
have been taken in fabricating suspended graphene sam- 



ples (Bolotin et al. 2008a Du et al. 2008b |. Finally we 



point out the crucial difference between graphene and 2D 
GaAs in phonon limited mobility. In the 2D GaAS sys- 
tem the acoustic phonon scattering is important below 
T — 100 K and polar optical phonon scattering becomes 
exponentially more important for T > 100 K whereas in 
graphene a resistivity linear in T is observed upto very 
high temperatures (w 1000 K) since the relevant opti- 
cal phonons have very high energy (w 2000 K) and are 
simply irrelevant for carrier transport. 



E. Other scattering mechanisms 

1. Midgap states 

The Boltzmann transport theory developed in 
Sec. |III.A| considered the limit of weak scattering. One 
can ask about the opposite limit of very strong scattering. 
The unitarity of the wavefunctions implies that a poten- 
tial scatterer can only cause a phase-shift in the outgoing 
wave. Standard treatment of s-wave elastic scattering 
gives the scattering time 



h 

Tk 



TTD{Ek\ 



■ sin^((5fc 



(3.34) 
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FIG. 25 (Color online) Left panel: Raman spectra (wave- 
length 633 nm) for (a) pristine graphene and (b) graphene 
irradiated by 500 eV Ne^ ions that are known to cause vacan- 
cies in the graphene lattice. Right panel: Increasing the num- 
ber of vacancies by ion irradiation caused a transition from 
the pristine graphene (where Coulomb scattering dominates) 
to the lower curves where scattering from vacancies dominate. 
Also shown is a fit to Eq. 3.36 from Stauber el al. (20071 that 



describes scattering off vacancies that have midgap states. 
Figures taken from Chen et al. (2009c I. 



where the conductivity is then given by the Einstein re- 
lation a = {2e^ / h)v-pk-pTi^^ . 

To model the disorder potential induced by a va- 
cancy, [Hentschel and Guinea ( |2007 ) assumed a circu- 
larly symmetric potential with 1/(0 < r < R') = oo, 
V{R' < r < R) = const, and V{R > r) = 0. This 
corresponds to a circular void of radius i?', and appro- 
priate boundary conditions are chosen to allow for zero 
energy states (also called mid-gap states). By matching 
the wavefunctions of incoming and outgoing waves, the 



scattering phase shift can be calculated as ( Guinea 2008 
Hentschel and Guinea| |2007[ ) 



arctan 



MkR') \ fe^p 
Yo{kR') 



1 



2 ln(A:i?') 



(3.35) 



where Jq{x) \Yi^{x)] is the zeroth order Bessel function of 
the first (second) kind. Expanding for small carrier den- 



sity, one then finds for the conductivity (Stauber et al. 



20071 



2p/ n 
nh rid 



In^ikpR'), 



(3.36) 



which other than the logarithmic factor, mimics the be- 
havior of charged impurities, and is linear in carrier den- 
sity. 

, ^, . , „„„„ 



(2009c 



In recent experimental work, Chen et al. 
radiated graphene with He and Ne ions to deliberately 
create large vacancies in the graphene sheet. They fur- 
ther demonstrated that these vacancies induced by ion 
irradiation gave rise to a strong D-peak in the Raman 
spectra, inferring that the absence of such a D-peak in 
the pristine graphene signalled the lack of such defects 
(Fig. [25]). Moreover, they demonstrated that while trans- 



port in pristine graphene is dominated by charged im- 
purities, after ion irradiation the electron scattering off 
these vacancies appears consistent with the theory in- 
cluding midgap states (Eq. 3.36). We mention that in 
this review, we consider only the case where the disorder 
changes graphene's transport properties without modi- 



fying its fundamental chemical structure (Hwang et al. 



2007b Schedin et al. 


2007) 


graphane ( 


Elias et al. 


2009 



The subject of transport in 



Sofo et al. 2007 ) and other 



chemical derivatives of graphene is beyond the scope of 



2009a b) 



(|Bostwick et al. 


2009 


Cheianov et al. 


IRobinson et al. 


2008 


Wehling et al. 



2. Effect of strain and corrugations 

While graphene is often assumed to be an atomically 
perfect 2D sheet, in reality, graphene behaves more like 
a membrane. When placed on a substrate, graphene 
will conform to the surface roughness developing ripples. 
Even without a substrate, experiments reveal significant 



deformations (Meyer et al. 2007), although the theo- 



retical picture is still contentious (Fasolino et al. 2007 



Pereira et al. 2009 Thompson-Flagg et al. 2009). It 



is nonetheless an important theoretical question to ad- 
dress the nature of electronic scattering off such ripples. 
Ripples, by their very nature, are correlated long-range 
fluctuations across the entire sample (i.e. most experi- 
ments measuring ripples calculate a height-height corre- 
lation function). Yet, for electronic transport, one would 
like to isolate a "single ripple" and calculate its scatter- 
ing cross-section (assuming that the rest of the sample 
is flat), and then treat the problem of electrons scatter- 
ing off ripples as that of random uncorrelated impurities 
with the cross-section of a single ripple. This was the 



approach followed by Katsnelson and Geim (2008) and 



Guinea (20081; Prada et al. (2010) 



With this qualitative picture in mind, one could esti- 
mate the transport time due to ripples as 



- w 2ttD(Ef)(V„V^ 

T 



(3.37) 



where Vq is the scattering potential caused by the strain- 
fields of a single ripple. 

Introducing a height field h{r) (that measures displace- 
ments normal to the graphene sheet), one can approxi- 
mate ( Katsnelson and Geim[ [2008 ) 



hvp 

a 



E 

qi,q2 



(^q— qi ^q^— q+q2 q2 ) 



[(q-qi) -qi] [(q-q2) •q2], 



(3.38) 



where a is the lattice spacing. Following |Ishigami| 
et al. ( |2007 ) ripple correlations can be parametrized as 



{[h{r) — /i(0)]^) = r^-^, where the exponent H provides 
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information about the origin of the ripples. An expo- 
nent 2H = 1 indicates that height fluctuation domains 
have short-range correlations implying that graphene 
conforms to the morphology of the underlying substrate, 
while 2H = 2 suggests a thermally excitable membrane 
only loosely bound by Van der Waals forces to the sub- 
strate 



non- universal minimum of conductivity observed in gated 
graphene devices at the charge neutrality point with 
the typical observed minimum conductivity being much 
larger than the universal prediction. 

Unfortunately gd is ill-defined, and depending on the 
theoretical methods and approximation schemes, many 



Ishigami et al. (2007) found experimentally that different universal results have been predicted (Aleiner 



2H 1.11 ±0.013 implying that graphene mostly con- 
forms to the substrate, but with some intrinsic stiffness. 



Katsnelson and Geim (|2008|) showed that this has con- 

= 1, 
For 



sequences for transport properties, where for 2H 
a - l/ln^[fcFa]; and for 2H > 1, a{n) - n^"-^. 
the special case of 2H = 2 (flexural ripples), this scat- 
tering mimics the long-range Coulomb scattering dis- 
cussed in Sec. |III.Al For the experimentally relevant case 
of 2H > 1, electron scattering off ripples would mimic 



short-range disorder also discussed in Sec. III.A Thus, 
ripple scattering in graphene for 2H « 1 mimic surface 



roughness scattering in Si MOSFET (Ando et al. 1982). 



We should emphasize that these conclusions are at best 
qualitative, since the approximation of treating the rip- 
ples as uncorrelated single impurities is quite drastic. A 
complete theory for scattering off ripples in graphene is 
an interesting, and at present, open problem. Ripple 
scattering effects on graphene transport have a formal 
similarity to the well-studied problem of interface rough- 
ness scattering effects on carrier transport in Si-Si02 2D 
electron systems dAdam et al.[ |2008bl |Ando et al.[ |1982 



Tracy et al. 2009). 



IV. TRANSPORT AT LOW CARRIER DENSITY 
A. Graphene minimum conductivity problem 

1. Intrinsic conductivity at the Dirac point 

One of the most discussed issues in the context of fun- 
damental graphene physics has been the so-called mini- 
mum (or minimal) conductivity problem (or puzzle) for 
intrinsic graphene. In the end, the graphene minimum 
conductivity problem turns out to be an ill-posed prob- 
lem, which can only be solved if the real physical system 
underlying intrinsic (i.e. undoped) graphene is taken into 
account. An acceptable and reasonably quantitatively 
successful theoretical solution of the minimum conduc- 
tivity problem has only emerged in the last couple of 
years, where the theory has to explicitly incorporate car- 
rier transport in the highly inhomogeneous electron-hole 
landscape of extrinsic graphene, where density fluctua- 
tions completely dominate transport properties for actual 
graphene samples. 

The graphene minimum conductivity problem is the 
dichotomy between the theoretical prediction of a uni- 
versal Dirac point conductivity ajj of undoped intrinsic 
graphene and the actual experimental sample-dependent 



et a;.[|1994[ [Peres et a/.[ |2006| [Tworzydlo et al 



[and Efetovt [2 0061 |Altlaiidl [2006^ 'Bardarson et a 
Fradkin, 1986; IFritz et aLl[2008 Kashuba, 20081 ' 



2007 
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and other values. The conductivity, 

cr(T, cli, Ei?, r, A, L^^), is in general a function of many 
variables: temperature (T), frequency (lu), Fermi energy 
or chemical potential (ep), impurity scattering strength 
or broadening (F), intervalley scattering strength (A), 
system size (L). The Dirac point conductivity of clean 
graphene, cr_D(0, 0, 0, 0, 0, 0), is obtained in the limit of 
all the independent variables being zero, and the result 
depends explicitly on how and in which order these 
limits are taken. For example, w — ^ and T — ^ limit 
is not necessarily interchangeable with the T — >■ and 
a; — ^ limit! In addition, the limit of vanishing impurity 
scattering (F — )■ 0) and whether F = or F 7^ also may 
matter. In the ballistic limit (F = 0), the mesoscopic 
physics of the system size being finite {1/L ^ 0) or 
infinite {l/L — 0) seems to matter. The intervalley 
scattering being finite (A 7^ 0) or precisely zero (A = 0) 
seems to matter a great deal because the scaling theory 
of localization predicts radically different results for an, 
ctd = for A 7^ 0, (Td = 00 for A = 0, in the presence of 
any finite disorder (F 7^ 0). 

A great deal of the early discussion on the graphene 
minimum of conductivity problem has been misguided 
by the existing theoretical work which considered the 
strict T — limit and then taking the w — >■ limit. 
Many theories claim ao = 4e^/(7rft,) in this limit, but 
the typical experimentally measured value is much larger 
(and sample-dependent), leading to the so-called "prob- 
lem of the missing pi". The limit HuIi^^q a{uj,T ~ 0) 
is, in fact, experimentally irrelevant since for experimen- 
tal temperatures (even 10 mK), ksT ^ fiw, and thus 
the appropriate limiting procedure for dc conductivity is 
ImiT^o o'i^ = OjT). There is an intuitive way of study- 
ing this limit theoretically, which, however, can only 
treat the ballistic (and therefore, the completely unre- 
alistic disorder- free) limit. Let us first put w = and 
assume fi = 0, i.e. intrinsic graphene. It is then easy to 
show that at T 7^ 0, there will be a finite carrier den- 
sity Ue = nil oc thermally excited from the graphene 
valence band to the conduction band. The algebraic 
dependence of thermal carrier density, rather than the 
exponentially suppressed thermal occupancy in semicon- 
ductors, of course follows from the non-existence of a 
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band gap in graphene. Using the Drude formula for dc 
conductivity, we write ud = ne^rjra cx T'^T{T)/ra{T), 
where r, m are respectively the relaxation time and the 
effective mass. Although graphene effective mass is zero 
due to its linear dispersion, an effective definition of effec- 
tive mass follows from writing e = hvpk — {h? k'^p) / {2m) , 
which leads to m oc ^/n oc T (which vanishes as T — >■ 0) 
by using k oc ^/n. This then leads to <td ^ Tt(T). 
In the ballistic limit, the only scattering mechanism is 
the electron-hole scattering, where the thermally excited 
electrons and holes scatter from each other due to mutual 
Coulomb interaction. This inelastic electron-hole scatter- 
ing rate l/r is given by the imaginary part of the self- 
energy which, to the leading order, is given by 1/t ~ T, 
leading to (Tjj ^ T{l/T) ^ a constant in the ballistic 
limit. There are logarithmic sub-leading terms which 
indicate that (Jd{T — ^ 0) grows logarithmically at low 
temperature in the ballistic limit. The conductivity in 
this picture, where interaction effects are crucial, is non- 
universal even in the ballistic limit, depending logarith- 
mically on temperature and becoming infinite at T = 0. 
The presence of any finite impurity disorder modifies the 
whole picture completely. More details along this idea 



can be found in the literature ( Foster and Aleiner 2008 



2009 



2008 



[Fritz et aij [20081 [Kashuba[ [2008[ [Miiller et al. 



2. Localization 

A fundamental mystery in graphene transport is the 
absence of any strong localization-induced insulating 
phase at low carrier density around the Dirac point, 
where kpl <ti 1 since fcp w at the charge neutrality 
point and the transport mean free path I is finite (and 
small). This is a manifest violation of the loffe-Reggel 
criterion which predicts strong localization for kpl ^ 1- 
By contrast, 2D semiconductor systems always go insu- 
lating in the low density regime. It is conceivable, but 
does not seem likely, that graphene may go insulating due 
to strong localization at lower temperatures. Until that 
happens, the absence of any signature of strong local- 
ization in graphene is a fundamental mystery deserving 
serious experimental attention. Two noteworthy aspects 
stand out in this context. First, no evidence of strong 
localization is observed in experiments that deliberately 



break the A-B sublattice symmetry ( Chen et al. 2009c I . 



Thus, the absence of localization in graphene cannot be 
attributed to the chiral valley symmetry of the Dirac 
fermions. Second, the opening of an intrinsic spectral 
gap in the graphene band structure by using graphene 



nanoribbons (Adam et al. 2008a Han et al. 20071 or 



biased BLG (Oostinga et al. 2008 Zhang et al. 2009c) 



immediately introduces an insulating phase around the 
charge neutrality point. These two features indicate that 
the insulating behavior in graphene and 2D semiconduc- 



tors is connected more with the existence of a spectral 
gap than with the quantum localization phenomena. 



3. Zero density limit 

It is instructive to think about the intrinsic conductiv- 
ity as the zero-density limit of the extrinsic conductivity 
for gated graphene. Starting with the Boltzmann theory 



high density result of Sec. Ill we see that 



(Jd = a{n 0) 



Coulomb scattering 
Ci zero-range scattering 



,(4.1) 



where the non-universal constant Ci is proportional to 
the strength of the short-range scattering in the system. 
We note that the vanishing of the Boltzmann conductiv- 
ity in the intrinsic zero-density limit for Coulomb scat- 
tering is true for both unscreened and screened Coulomb 
impurities. The non-vanishing of graphene Boltzmann 
conductivity for zero-range (5-function scattering poten- 
tial in the zero carrier density intrinsic limit follows di- 
rectly from the gapless linear dispersion of graphene car- 
riers. We emphasize, however, that ajj is non-universal 
for zero-range scattering. 

For further insight into the zero-density Boltzmann 
limit the T — for a let's consider Eq. (3.3 1. In gen- 
eral, T~'^{E) ^ D{E) since the availability of unoccupied 
states for scattering should be proportional to the den- 
sity of states. This immediately shows that the intrinsic 
limit, Ep{n — > 0) — > 0, is extremely delicate for graphene 
because D{E — >■ 0) — >■ 0, and the product Dt becomes 
ill-defined at the Dirac point. 

We emphasize in this context, as discussed in Sec. |l] 
that as a function of carrier density (or gate voltage), 
graphene conductivity (at high carrier density) is quali- 
tatively identical to that of semiconductor-based 2DEG. 
This point needs emphasis because it seems not to be 
appreciated much in the general graphene literature. In 
particular, a{n) ^ rf for both graphene and 2DEG 
with a = 1 for graphene at intermediate density and 
a w 0.3 to 1.5 in 2DEG depending on the semiconduc- 
tor system. At very high density, a ~ (or even nega- 
tive) for both graphene and 2DEG. The precise nature 
of density dependence (i.e. value of the exponent a) de- 
pends strongly on the nature of scattering potential and 
screening, and varies in different materials with graphene 
{a « 1) falling somewhere in the middle between Si- 
MOSFETs (a « 0.3) and modulation doped 2D n-GaAs 
{a sa 1.5). Thus, from the perspective of high-density 
low-temperature transport properties, graphene is sim- 
ply a rather low-mobility (comparable to Si MOSFET, 
but much lower mobility than 2D GaAs) 2D semicon- 
ductor system. 
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4. Electron and hole puddles 

The low-density physics in both graphene and 2D semi- 
conductors is dominated by strong density inhomogene- 
ity ( "puddle" ) arising from the failure of screening. This 
inhomogeneity is mostly due to the random distribution 
of unintentional quenched charged impurity centers in 
the environment. (In graphene, ripples associated with 
either intrinsic structural wrinkles or the substrate inter- 
face roughness may also make contribution to the inho- 
mogeneity.) At low density, the inhomogeneous puddles 
control transport phenomena in graphene as well as in 
2D semiconductors. Inhomogeneous puddles would form 
also in doped 3D semiconductors at low carrier densities 



(Shklovskii and Efros 1984). 



In section IV B we discuss the details of electron-hole 
puddle formation in graphene around the charge neu- 
trality point and describe its implications for graphene 
transport properties. Here we emphasize the qualita- 
tive difference between graphene and 2D semiconductors 
with respect to the formation of inhomogeneous puddles. 
In 2D semiconductors, depending on whether the sys- 
tem is electron-doped or hole-doped, there are only just 
electron- or just hole-puddle. At low density, n « 0, 
therefore most of the macroscopic sample has little finite 
carrier density except for the puddle regime. From the 
transport perspective, the system becomes the landscape 
of mountains and lakes for a boat negotiating a hilly lake. 
When percolation becomes impossible, the system be- 
comes an insulator. In graphene, however, there is no 
gap at the Dirac point, and therefore, the electron (hole) 
lakes are hole (electron) mountains, and one can always 
have transport even at zero carrier density. This pic- 
ture breaks down when a spectral gap is introduced, and 
gapped graphene should manifest an insulating behav- 
ior around the charge neutrality point as it indeed does 
experimentally. 



5. Self-consistent theory 

The physical puddle picture discussed above enables 
one to develop a simple theory for graphene transport at 
low densities using a self-consistent approximation where 
the graphene puddle density is calculated by consider- 
ing the potential and density fluctuations induced by the 
charged impurities themselves. Such a theory was devel- 



oped by Adam et al. (2007). The basic idea is to realize 
that at low carrier density \n\ < \ni\, the self-consistent 
screening adjustment between the impurities and the car- 
riers could physically lead to an approximate pinning of 
the carrier density at n = n* « n.^. A calculation within 



the RPA approximation yields ( Adam et al. 2007 1 



2rlC^^^{rs,a = AdV^), (4.2) 



(2- 



1 + 2r, 



+ (1 + 2r,a) e^''^^{Ei[2r,a] - Ei[a{l + 2r,)]), 

where Ei{z) — t^^e^*'dt is the exponential integral 
function. This density pinning then leads to an approx- 
imately constant minimum graphene conductivity which 
can be obtained from the high-density Boltzmann the- 
ory by putting in a carrier density of n* . This simple 
intuitive self-consistent theory is found to be in surpris- 
ingly good agreement with all experimental observations 
(I Adam et al.\ 120071 IChen et all |2008a|). In the next 



section we describe a more elaborate density functional 
theory and an effective medium approximation to cal- 
culate the puddle electronic structure and the resultant 



Sarma 2008). 



transport properties (Rossi et al. 2009 Rossi and Das 



B. Quantum to classical crossover 

The starting point for the quantum transport proper- 
ties at the Dirac point discussed in Sec. |II.CTT] is the bal- 
listic universal minimum conductivity CTmin — 4e^/(7r/i) 
for clean graphene. The addition of disorder i.e. includ- 



ing potential fluctuations (given by Eq. 2.10) that are 
smooth on the scale of the lattice spacing increases the 
conductivity through weak anti-localization. This pic- 
ture is in contrast to the semi-classical picture discussed 
above where the transport properties are calculated at 
high density using the Boltzmann transport theory and 
the self-consistent theory is used to handle the inhomo- 
geneities of the carrier density around the Dirac point. 
This theory predicts that the conductivity decreases with 
increasing disorder strength. Given their vastly different 
starting points, it is perhaps not surprising that the two 
approaches disagree. 

A direct comparison between the two approaches has 
not been possible mainly because the published predic- 
tions of the Boltzmann approach include screening of the 
Coulomb disorder potential, whereas the fully quantum- 
mechanical calculations are for a non-interacting model 
using Gaussian disorder. Notwithstanding the fact that 
screening and Coulomb scattering play crucial roles in 
transport of real electrons through real graphene, the 
important question of the comparison between quantum 
and Boltzmann theories, even for Gaussian disorder, was 



addressed only recently Adam et al. (2009a I where they 



considered non-interacting Dirac electrons at zero tem- 
perature with potential fluctuations of the form shown 
in Eq. |2.10[ They numerically solved the full quantum 
problem for a sample of finite size L ^ ^ (where ^ is the 
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FIG. 26 (Color online) Resistance, R — 1/G (left) and con- 
ductivity (rigiit) obtained using a = [WdR/dL]~^ , as a func- 
tion of sample length L. The three curves shown are for 
W/^ = 200, Ko ^ 2 and 7r<^ = 0, 0.25, and 1 [from top 
to bottom (bottom to top) in left (right) panel] . Dashed lines 
in the right panel show dcr/dlnL — Ae^ /ixh. The inset in the 
left panel shows the crossover to diffusive transport (L 2> 0- 




FIG. 27 (Color online) Semiclassical conductivity o' = 
limi_>oo [a"(i) — tt"'^ ln(L/f )] versus disorder strength at the 
Dirac point (left) and at carrier density nn = Kq / (irS,)'^ , cor- 
responding to the edge of the minimum conductivity plateau 
of [Adam et al.\ ( |2007| ) (right). Data points are from the nu- 
merical calculation for L — 50^ and the (solid) dashed curves 
represent the (self-consistent) Boltzmann theory. 



correlation length of the disorder potential in Eq. 2.10 ), 
for a range of disorder strengths parameterized by Kq 
Typical results for the quantum transport are shown in 



^ Quantum effects are a small correction to the conductivity only 
if the carrier density n is increased at fixed sample size L. This 
is the experimentally relevant limit. If the limit L — ^ oo is taken 
at fixed n, quantum efi'ects dominate (see Eq. |2.8| l, where the 
semiclassical theory does not capture the logarithmic scaling of 
conductivity with system size. Here, we are not considering the 
conceptually simple question of how quantum transport becomes 
classical as the phase coherence length decreases, but the more 
interesting question of how this quantum-Boltzmann crossover 
depends on the carrier density and disorder strength. 



Fig. 26 For L < the transport is ballistic and the con- 
ductivity given by the universal value (Tmin = 4e^/(7r/i). 
For L ^ ^, one is in the difltusive transport regime. For 



the diffusive regime, Adam et al. (2009a) demonstrated 



that away from the Dirac point, both the Boltzmann the- 
ory and the full quantum theory agree to leading order 
with 



a{n) = 



Koh 



(4.3) 



While this agreement is perhaps not surprising, it vali- 
dates the assumptions of both theories and demonstrates 
that they are compatible at high carrier density. More 
interesting are the results at the Dirac point. General- 
izing the self-consistent Boltzmann theory to the case of 
a Gaussian correlated disorder potential (Eq. [2.10 ) , one 
finds 



2e^ 
nh 



exp 



-Kq 

2tt 



h 



27r 



(4.4) 



where Ii is the modified Bessel function. Shown in 



the left panel of Fig. 27 is a comparison of the nu- 
merical fully quantum Dirac point conductivity where 
the weak antilocalization correction has been subtracted 
a' = limi_j.oo[o'(i) — 7r~^ ln(L/^)] with the semiclassical 



result (Eq. 4.4) 



The right panel shows the conductivity slightly away 
from the Dirac point (i.e. at the edge of the minimum 
conductivity plateau). The numerical calculations at 
the edge of the plateau are in good quantitative agree- 
ment with the self-consistent Boltzmann theory. At the 
Dirac point, however, the quantum conductivity tIKq) 
is found to increase with Kq for the entire parameter 
range considered, which differs from the Boltzmann the- 
ory at small Kq. At large Kq the numerical data fol- 
low the trend of the self-consistent theory which predicts 
a ~ 2e^ K^^"^ / {-kK) for Kq ^ 10. This implies that even 
at the Dirac point, for large enough disorder, the trans- 
port is semi-classical and described by the self-consistent 
Boltzmann transport theory. 

For smaller ifg. Fig. [27] shows that upon reducing 
Kq below unity, the conductivity first decreases sharply 
consistent with a renormalization of the mean free path 
due to the ultraviolet logarithmic divergences discussed 
in Sec. |II.C.4| Upon reducing Kq further, the Dirac 
point conductivity satura tes at t he ballistic value CTmin = 
ie^/nh (discussed in Sec. llI.B.2|) 



In a closely related work, Lewenkopf et al. (2008) did 



numerical simulations of a tight-binding model to ob- 
tain the conductivity and shot noise of graphene at the 
Dirac point using a recursive Green's function method. 
This method was then generalized to calculate the metal- 
insulator transition in graphene nanoribbons where, as 
discussed in Sec. |II.C.2[ edge disorder can cause the An- 
derson localization of electrons (Mucciolo et al. 2009). 
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The important conclusion of this section is that it pro- 
vides the criteria for when one needs a full quantum me- 
chanical solution and when the semi-classical treatment 
is sufficient. For either sufficiently weak disorder, or when 
the source and drain electrodes are closer than the scat- 
tering mean free path, then the quantum nature of the 
carriers dominates the transport. On the other hand, for 
sufficiently large disorder, or away from the Dirac point, 
the electronic transport properties of graphene are semi- 
classical and the Boltzmann theory correctly captures the 
most of graphene's transport properties. 



C. Ground state in the presence of long-range disorder 

In the presence of long-range disorder that does not 
mix the degenerate valleys the physics of the graphene 
fermionic excitations is described by the following Hamil- 
tonian: 



dovich and Popov 1972). In the context of graphene 



j d\VD{v)^L'^ro. (4.5) 



e 

2k. 



where vp is the bare Fermi velocity, are the 

creation annihilation spinor operators for a fermionic ex- 
citation at position r and pseudospin a, cr is the 2D 
vector formed by the 2x2 Pauli matrices Gx and Gy 
acting in pseudospin space, ^ is the chemical potential, 
1 is the 2x2 identity matrix, k is the effective static 
dielectric constant equal to the average of the dielec- 
tric constants of the materials surrounding the graphene 
layer, V{\y — r'|) — l/||r — r'| is the Coulomb interaction 
and Vd{y) is the bare disorder potential. The Hamilto- 
nian (4.5) is valid as long as the energy of the fermionic 



excitations is much lower than the graphene bandwidth 



3 eV. Using (4.5 1 if we know Vd we can characterize 
the ground state carrier density probability close to the 
Dirac point. In this section we focus on the case when Vd 
is a disorder potential whose spatial autocorrelation de- 
cays algebraically, such as the disorder induced by ripples 
or charge impurities. 



1. Screening of a single charge impurity 

The problem of screening at the Dirac point of a sin- 
gle charge impurity placed in (or close to) the graphene 
layer illustrates some of the unique features of the screen- 
ing properties of massless Dirac fermions. In addition 
the problem provides a condensed matter realization of 
the QED phenomenon of "vacuum polarization" induced 
by an external charge ( Case[ 195^ 



Darwin 



1928; Gor- 



|don[ 1928 Pomeranchuk and Smorodinsky[ |1945( |Zel 



the problem was first studied by [Divincenzo and Mele 



(1984) 


and recently more in detail by several authors 


(Biswas et al.\ 2007 


Fistul and Efetov 2007 Fogler 


et al.\ 


20071 INovikov 


2007b Pereira et al. 2007 Shy- 


tov et al. 2007a|b| Terekhov et al. 2008). The parame- 



ter j5 = Zfr/lnhvf) = Zr^ quantifies the strength of the 
coupling between the Coulomb impurity and the massless 
Dirac fermions in the graphene layer. Neglecting e — e 
interactions for |/3| < 1/2 the Coulomb impurity induces 
a screening charge that is localized on length scales of the 
order of the size of the impurity itself (or its distance d 
from the graphene layer). Even in the limit < 1/2 the 
inclusion of the e— e interactions induces a long-range tail 
in the screening charge with sign equal to the sign of the 



charge impurity (Biswas et al. 2007). For \/3\ > 1/2 the 



Coulomb charge is supercritical, the induced potential is 



singular (Landau and Lifshitz 1977), and the solution of 



the problem depends on the regularization of the wave- 
function at the site of the impurity, r — > 0. By setting 
the wavefunction to be zero at r = a the induced electron 
density in addition to a localized {S{r)) term, acquires a 
long-range tail ~ 1/r^ (with sign opposite to the sign 
of the charge impurity) ( Novikovj 2007a| Pereira et al. 



2007 Shytov et al.\ 2007b) and marked resonances ap- 



pear in the spectral density ( Fistul and Efetov 2007 Shy- 



tov et al. 2007a ) that should also induce clear signatures 
in the transport coefficients. Up to now neither the os- 
cillations in the LDOS nor the prdicted signatures in the 



conductivity (Shytov et al. 2007a) have been observed 



experimentally. It is likely that in the experiments so far 
the supercritical regime |/3| > 1/2 has not been reached 
because of the low Z of the bare charge impurities and 
renormalization effects. [Fogler et al. ( |2007 ) pointed out 
however that the predicted effects for |/3| > 1/2 are in- 
trinsic to the massless Dirac fermion model that however 
is inadequate when the the small scale cut-off niin[(i, a 
smaller than ar^\/^. 



IS 



2. Density functional theory 

Assuming that the ground state does not have long 



range order ( Dahal et al. 2006 Min et al. 2008b Peres 



et al. 2005 ) a practical approach to many-body problems 



is the Density Functional Theory, DFT, (Giuliani and 



Vignalej |2005l [Hohcnberg and Kohn[ |1964[ |Kohn[ |1999 



Kohn and Sham 1965). In this approach the interaction 
term in the Hamiltonian is replaced by an effective Kohn- 
Sham potential Vks that is a functional of the ground- 
state density n(r) — J2cr '^la'^ra- 

d\'<Sjl^VKs[n{r)]^ro. (4.6) 
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The Kohn-Sham potential is given by the sum of the 
external potential, the Hartree part of the interaction, 
Vh, and an exchange-correlation potential, V^c that can 
only be known approximately. In its original form V^c is 
calculated within the local density approximation, LDA 
(Kohn and Sham 1965), i.e., Vxc is calculated for a uni- 
form liquid of electrons. The DFT-LDA approach can 
be justified and applied to the study of interacting mass- 

For graphene 



less Dirac fermions (Polini et al. 2008b I 



the LDA exchange-correlation potential, within the RPA 
approximation is given with very good accuracy by the 



following expression ( Barlas et al. 


2007 


Gonzalez et al. 


1999VHwang et aL';2007c; 


Katsnelson 


2006 


Mishchenko 


2007i|Polini et aLn2008b, 


iVafek 2007 


I: 



^ ^/^\ sgn(n) In 

4 V47rn 

\/Tr\n\ sgn(n) In -— 



(4.7) 



where g is the spin and valley degeneracy factor (g — A), 
and kc is an ultraviolet wave-vector cut-off, fixed by the 
range of energies over which the pure Dirac model is valid. 
Without loss of generality we can use kc — 1/a, where 
a is the graphene lattice constant, corresponding to an 
energy cut-off Ec « 3 eV . Equation (4.7) is valid for 
kp = ^7r|n| ^ kc. In (4.7) ^ is a constant that depends 
on rs giben by (Polini et al. 2008b): 



1 



+00 



dx 



2 Jo {l + x^)^Vl+^ + 7:gr,/8y 



(4.8) 



The terms on the r.h.s. of (4.7) are the exchange and the 
correlation potential respectively. Notice that exchange 
and correlation potentials have opposite signs. In Fig. 28 
the exchange and correlation potential and their sum, 
V^c, are plotted as a function of n for = 0.5. We 
see that in graphene the correlation potential is smaller 
than the exchange potential but, contrary to the case 
of regular parabolic band 2DEGs, it is not negligible. 



However from (4.7) we have that exchange and correla- 
tion scale with n in the same way. As a consequence 
in graphene the correlation potential can effectively be 
taken into account by simply rescaling the coefficient of 
the exchange potential. In Fig. |28] the green dotted line 
shows Vxc for a regular parabolic band 2DEG with ef- 
fective mass 0.067me in a background with k ~ A. The 
important qualitative difference is that V^c in graphene 
has the opposite sign than in regular 2DEGs: due to in- 
terlayer processes in graphene the exchange correlation 
potential penalizes density inhomogeneities contrary to 
what happens in parabolic-band electron liquids. 



Using the DFT-LDA approach [Polini et al.\ ( |2008b[ ) 
calculated the graphene ground-state carrier density for 
single disorder realizations of charge impurities and small 
samples (up to 10 x 10 nm). The size of the samples 
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FIG. 28 (Color online). Exchange, black solid line, and RPA 
correlation potentials blue dashed line, as functions of the 
density n for = 0.5. The magenta dash-dotted line shows 
the full exchange-correlation potential, Vxc- The green dotted 
line is the quantum Monte Carlo exchange-correlation poten- 



tial of a standard parabolic-band 2D electron gas ( Attaccalite 
et al.\ I2OO2 



with effective mass 0.067 me placed in back- 
ground with dielectric constant 4. Adapted from |Polini et aL\ 
( j2008b| ). 



is limited by the high computational cost of the ap- 
proach. For single disorder realizations the results of 



Polini et al. (2008b) show that, as it had been predicted 
(Hwang et al. 2007a), at the Dirac point the carrier 



density breaks up in electron-hole puddles and that the 
exchange-correlation potential suppresses the amplitude 
of the disorder induced density fluctuations. Given its 
computational cost, the DFT-LDA approach does not 
allow the calculation of disordered averaged quantities. 



3. Thomas Fermi Dirac theory 

An approach similar in spirit to the LDA-DFT is the 



Thomas-Fermi, TF, (Fermi 1927 Giuliani and Vignale 



2005[ |Spruch[ |1991[ |Thomas[ |1927[ ) theory. Like DFT 
the TF theory is a density functional theory: in the 
Thomas-Fermi theory the kinetic term is also approx- 
imated via a functional of the local density n(r). By 
Thomas-Fermi-Dirac, TFD, theory we refer to a mod- 
ification of the TF theory in which the kinetic func- 
tional has the form appropriate for Dirac electrons and 
in which exchange-correlation terms are included via the 
exchange-correlation potential proper for Dirac electron 
liquids as described above for the DFT-LDA theory. The 
TF theory relies on the fact that if the carrier density 
varies slowly in space compared to the Fermi wavelength, 
then the kinetic energy of a small volume with density 
n(r) is equal, with good approximation, to the kinetic 
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energy of the same volume of a homogeneous electron 
liquid with density n = n(r). The condition for the va- 



lidity of the TFD theory is given by the inequality ( Brey 



and Fertigl |2009b| [Giuliani and Vignalel [20051 ) 



|Vrn(r) 
n(r) 



< fcF(r). 



(4.9) 



Whenever inequality (4.9) is satisfied the TFD theory 
is a computationally efhcient alternative to the DFT- 
LDA approach to calculate the ground state properties of 
graphene in presence of disorder. The energy functional, 
E[n] in the TFD theory is given by: 



E[n\ =hvp 



2a/^ 



(frsgn{n)\n\^^'^ 

,^,n(r)n(rO^ /dVt4Jn(r)]n(r) 
F-r'l J 



+ r, / d'^rVD(r)n(r) - / d^rnir) 

' nvf ' 



(4.10) 



where the first term is the kinetic energy, the second is 
the Hartree part of the Coulomb interaction, the third 
is the term due to exchange and correlation, and the 
fourth is the term due to disorder. The expression for 



the exchange-correlation potential is given in Eq. (4.71 



The carrier ground state distribution is then calculated 



by minimizing E[n\ with respect to n. Using (4.10) the 
condition SE/6n = requires 



sgn{n)^/\Ti 
+ rsVDir) 



2 



n(r') 



= 0. 



(4.11) 



Equation (4.111 well exemplifies the nonlinear nature of 



screening in graphene close to the Dirac point: because in 
graphene, due to the linear dispersion, the kinetic energy 
per carrier, the first term in (4.11 ), scales with y/n when 
(n) = the relation between the density fluctuations 
5n and the external disorder potential is not linear even 
when exchange and correlation terms are neglected. 

We now consider the case when the disorder potential 
is due to random Coulomb impurities. In general the 
charge impurities will be a 3D distribution, C(r), how- 
ever Si02 we can assume to a very good approximation 
C(r) to be effectively 2D. The reason is that for normal 
substrates such as Si02 the charge traps migrate to the 
surface of the oxide, moreover any additional impurity 
charge introduced during the graphene fabrication will 
be located either on the graphene top surface or trapped 
between the graphene layer and the substrate. We then 
assume C(r) to be an effective 2D random distribution 
located at the average distance d from the graphene layer. 
An important advantage of this approach is that it limits 
to two the number of unknown parameters that enter the 



theory: charge impurity density, ni,„p, and d. With this 
assumption we have: 



C(r') 



r — r 



'12 



•d2]l/2 



(4.12) 



The correlation properties of the distribution C(r) of the 
charge impurities are a matter of long-standing debate 
in the semiconductor community. Because the impuri- 
ties are charged one would expect the positions of the 
impurities to have some correlation, on the other hand 
the impurities are quenched (not annealed), they are ei- 
ther imbedded in the substrate or between the substrate 
and the graphene layer or in the graphene itself. This 
fact makes very difficult to know the precise correlation 
of the charge impurity positions but it also ensures that 
to good approximation the impurity positions can be as- 
sumed to be uncorrelated: 

(C(r)) - 0; (C(ri)C(r2)) = n„„p5(r2 - ri). (4.13) 

where the angular brackets denote averaging over disor- 
der realizations. A non zero value of (C(r)) can be taken 
into account simply by a shift of the chemical potential 
/i. It is easy to generalize the theory to correlated impu- 
rities (e.g. impurity clusters) if the correlation function 
is known. 

The parameters Tiimp and d that enter the theory are 
reliably fixed by the transport results, see Sec. |III.A[ at 
high doping. Transport results at high density indicate 
that d is of the order of 1 nm whereas riimp varies depend- 
ing on the sample quality but in general is in the range 
n-imp = 10^° — 10^^ cm~2, where the lowest limit applies 
to suspended graphene. The distance d is the physical 
cutoff for the length-scale of the carrier density inhomo- 
geneities. Therefore to solve numerically Eq. (4.11) one 



can use a spatial discretization with unit step of the order 
of d. For the TFD results presented below it was assumed 
d — 1 nm and therefore a spatial step Ax = Ay = 1 nm 
was used. 

shows the TFD results for the carrier den- 



Fig. 29 



sity distribution at the Dirac point in the presence of 
charge impurity disorder for a single disorder realiza- 
tion with riinip = 10^2 cm~2 and k = 2.5 correspond- 
ing to graphene on Si02 with the top surface in vac- 
uum (or air). It is immediately clear that as was pre- 



dicted (Hwang et al. 2007a), close to the Dirac point 



the disorder induced by the charge impurities breaks up 
the carrier distribution in electron (n > 0) and hole 
(n < 0) puddles. The electron-hole puddles are sepa- 
rated by disorder-induced p — n junctions, PNJ. Apart 
from the PNJ the carrier density is locally always dif- 
ferent from zero even though the average density, (n), 
is set equal to zero. For this reason, in the presence of 
disorder it is more correct to refer to the value of the 
gate voltage for which (n) = as the charge neutral- 
ity point, CNF, rather than Dirac point: the presence of 



46 



long-range disorder prevents the probing of the physical 
properties of the Dirac point, i.e. of intrinsic graphene 
with exactly half-filling, zero-density, everywhere. The 
important qualitative results that can be observed even 
for a single disorder by comparing the results of panels 
(b) and (c) is that the exchange-correlation term sup- 
presses the amplitude of the density fluctuations. This 
fact is clearly visible in panel (d) from which we can see 
that in the presence of exchange-correlation the density 
distribution is much narrower and more peaked around 
zero. This result, also observed in the DFT-LDA results 



(Polini et al. 2008b), is a consequence of the the fact 



that as discussed in Sec. IV. C. 2 the exchange-correlation 
potential in graphene, contrary to parabolic-band Fermi 
liquids, penalizes density inhomogeneities. 
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FIG. 29 (Color online). TFD results as a function of posi- 
tion for a single disorder realization, (a) Bare disorder poten- 
tial, Vd- (b) Carrier density obtained neglecting exchange- 
correlation terms, (c) Carrier density obtained including 
exchange-correlation terms, (d) Probability density distri- 
bution at the CNF with exchange (solid line) and without 
(dashed line). Adapted from Rossi and Das Sarma (20081 



In the presence of disorder, in order to make quantita- 
tive predictions verifiable experimentally it is necessary 
to calculate disordered averaged quantities. Using TFD, 
both the disorder average {X) of a given quantity X and 
its spatial correlation function 



the TFD theory both the spatial correlation function of 
the screened potential, V^sc, and carrier density are found 
to decay at long distance as l/r'^. This is a consequence 
of the weak screening properties of graphene and was 



pointed out in (Adam et al.\ 2007 Galitski et al. 20071. 



From the spatial correlation functions rirms and ^ = ^„ 



are extracted. Fig. 30 (a), (b) show the calculated 



and ^ at the Dirac point as a function of riinip- The dis- 
order averaged results show the effect of the exchange- 
correlation potential in suppressing the amplitude of the 
density inhomogeneities, rii-ms, and in slightly increas- 
ing their correlation length. The effect of the exchange- 
correlation potential increases as riimp decreases. At the 
Dirac point, the quantity ^ can be interpreted as the ef- 
fective non-linear screening length. Fig. |30| shows that f 
depends weakly on nimp. The reason is that ^ only char- 
acterizes the spatial correlation of the regions in which 
the density is relatively high. If a puddle is defined as 
a continuous region with same sign charges then at the 
CNP the puddles have always a size of the same order 
of the system size. Inside the puddles there are small 
areas with high density and size ^ of the order of tens 
of nanometers for typical experimental conditions, much 
smaller than the system size, L. This picture is confirmed 




5X\v)^{{X{v)-{X)){X{Q)-{X))) 



(4.14) 



FIG. 30 (Color online). TFD disorder averaged results. 
The solid (dashes) lines show the results obtained includ- 
ing (neglecting) exchange and correlation terms (a) nrms 
and, (b) ^ as a function of nimp. (c) Area, Aq, over which 
|ri(r) — (n)| < n^ms /lO. (d) A verage excess charge SQ vs. 
"-imp. Adapted from Rossi and Das Sarmar ( |2008] ) 



can be efficiently calculated. For conditions typical in 
experiments 500 disorder realizations are sufficient. From 
(5X^(r) one can extract the quantities: 



X„„, = V(('5^(0))2); (,x = FWHM oi{{6X{v)f)- 

(4.15) 

respectively as the root mean square (rms) and the typ- 
ical spatial correlation of the fluctuations of X. Using 



in Fig. 30 (c) in which the disorder averaged area frac- 
tion, Ao, over which \n{v) — {n) \ < Urms/^O is plotted as 
a function of Uimp- As Uimp decreases increases reach- 
ing more than 1/3 at the lowest impurity densities. The 
fraction of area over which |n(r) — (n) | is less than 1/5 of 
Urms surpasses 50 % for riimp ^ 10^° cm~^. Fig. 30 (b) 
shows the average excess charge SQ = UrmsT^^'^ at the 
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Dirac point as a function of ni, 



Notice that as de- 



fined JQ, especially at low nimp, grossly underestimates 
the number of charges both in the electron-puddles and 
in the small regions of size ^. This is becuase in the 
regions of size ^ the density is much higher than nrms 
whereas the electron-hole puddles have a typical size 
much larger than ^. Using for the small regions the esti- 
mate |Vn(r)|/?i = 1/^ and the local value of n inside the 
regions and, for the electron-hole puddles the estimates: 



|Vn(r)| 



(4.16) 



we find that the inequality (4.9 1 is satisfied guaranteeing 



the validity of the TFD theory even at the Dirac point. 

As we move away from the Dirac point more of the area 
is covered by electron (hole) puddles. However the den- 
sity flucutuations remain large even for relatively large 
values of Vg. This is evident from Fig. 31 where the 



probability distribution P(n) of the density for different 
values of (n) is shown in panel (a) and the ratio n^ms/in) 
as a function of (n) is shown in panel (b). The probabil- 



ity distribution P{n) is non-Gaussian (Adam et al. 2007 



2009b Galitski et al. 20071. For density {n 



< 



P{n) does not exhibit a single peak around (n) but rather 
a bimodal structure with a strong and narrow peak arond 
zero. The double peak structure for finite Vg provides di- 
rect evidence for the existance of puddles over a finite 
voltage range. rirms/(?^) decreases with (n), a trend that 
is expected and that has been observed indirectly in ex- 
periments by measuring the inhomogenous broadening of 
the quasiparticle spectral function (iHong et al. 2009a). 
In the limit Ts ^ 1 it is possible to obtain analytic 



results using the TFD approach (Fogler 2009). The first 
step is to separate the inhomogeneities of the carrier den- 
sity and screened potential in slow, n, Vsc and fast com- 
ponents, Sn, SVsc- 

n{r) = n{r) + Sn{r) , V^r) = Kc(r) + (5Kc(r) , (4.17) 

where n and Vsc contain only Fourier harmonics with 
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FIG. 31 (Color online), (a) Density distribution averaged 
over disorder for different values of the applied gate voltage 



assuming k — 2.5, d = 1 nm and Uir, 



(b) 



s/{n) as a function of (n) for d — 1 nm and different 



values of rii- 



circles. 



= 1.5 X 10' 



riimp = 10 cm ; triangles, rump 
from Rossi and Das Sarma ( 2008 1 



5 X 10" cm 



squares. 
Adapted 



fc < A where 1/A is the spatial scale below which the spa- 
tial variation of n and Vsc are irrelevant for the physical 
properties measured. For imaging experiments 1 /A is the 
spatial resolution of the scanning tip and for transport 
experiments 1/A is of the order of the mean free path. 
Let the Zimp = l/2rsy/ni^p and R the non- linear screen- 
ing length. It is assumed that / < 1/A ^ i?. With these 
assumptions and neglecting exchange-correlation terms, 
from the TFD functional in the limit A kp and small 
non-linear screening terms compared to the kinetic en- 
ergy term it follows ( Fogler[ |2009[ ) : 

, Vsc \Vsc\ sgn (Vsc) 
n{Vsc)^' 



In- 



Vs, 



iVscl > fivA . 



iT{hvY 2P "' hvA 

(4.18) 

with Vsc given, in momentum space, by the equation: 

2T:rshvF 



Vsc(k) = Voik) 



-n(k) 



(4.19) 



where we have assumed for simplicity d = 0. Equation 



(4.19) can be approximated by the following asymptotic 



expressions: 



Vscih) = 



Vb(k) 



kR 
1 + kR ' 



fci?> 1, 
fci?< 1. 



(4.20) 



Equations (4.181, (4.19) define a nonlinear problem that 
must be solved self-consistently and that in general can 
only be solved numerically. However in the limit ^ 1 
an approximate solution with logarithmic accuracy can 
be found. Let Kq be the solution of the equation: 



K„ - \n{l/{ArsKo)). 



(4.21) 



Kq is the expansion parameter. To order 0{Kf^^) Vsc 
can be treated as a Gaussian random potential whose 
correlator, 5Vs^{r) can be calculated using (|4.20) to find 
(I Adam et al.\ l20071 iFoglerl [20091 [GaTitski et aimOOT I 



Kir) = SVscir) 




with R = l/{4rsKo). Using ( [4.22^ and ( |4.18[ ) we can find 
the correlation function for the carrier density ( Fogler[ 
2009|: 



'27rZ4 




K{r) 

m 



(4.23) 



The correlation functions given by Eqs. (4.22) and (4.23) 
are valid in the limit ^ 1 but are in qualitative agree- 
ment also with the numerical results obtained for r. sa 1 



( Rossi and Das Sarma 2008 ) 
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The location of the disorder induced PNJ is identified 
by the isohnes n(r) = 0, or equivalently Kc(r) = 0. The 
CNP corresponds to the "percolation" threshold in which 
exactly half of the sample is covered by electron puddles 
and half by hole puddles (notice that conventionally the 
percolation threshold is defined as the condition in which 
half of the sample has non-zero charge density and half 
is insulating and so the term percolation in the context 
of the graphene CNP has a slight different meaning and 
does not imply that the transport is percolative) . At the 
percolation threshold all but one PNJ are closed loops. 
Over length scales d such that 1/A ^ d <C R, Vsc is 



(2008) found: 



6Ef = -a^iV^h)'-' (4.24) 
where a is a constant estimated to be approximately 



equal to 9.23 eV. Recently Gibertini et al. (2010) have 



used the DFT-LDA to study the effect of ripples on the 
spatial carrier density fluctuations. Transport theories in 
the presence of topological disorder were considered by 



Cortijo and Vomediano 


( 


2007 


); 


Cortijo and Vozmediano 


(2009 


) and 


Herbut et al. 


(2008 





logarithmically rough (Eq. (4.18)) and so the PNJ loops 

3/2 ( ,Kondev, 



of diameter d have fractal dimension Dh 

2000). At larger d the spatial correlation of 



et al. 



decays rapidly (Eq. (4.18)) so that for d larger than i?, 
Dh crosses over to the standard uncorrelated percolation 
exponent of 7/4 ( Isichenkoj 1992). 



5. Imaging experiments at the Dirac point 

The first imaging experiments using scanning tunnel- 
ing microscopy, STM, were done on epitaxial graphene 



Brar et al. 


2007 


Rutter et al. 


2007 



4. Effect of ripples on carrier density distribution 

When placed on a substrate graphene has been shown 



( Geringer et al. 2009 Ishigami et al. 2007 ) to follow with 



good approximation the surface profile of the substrate 
and therefore to have a finite roughness. For graphene on 
Si02 the standard deviation of the graphene height, h, 
has been measured to be 5h « 0.19 nm with a roughness 
exponent 2H « 1. More recent experiments (Geringer 



were able to image the atomic structure of graphene and 
reveal the presence of in-plane short-range defects. So 
far one limitation of experiments on epitaxial graphene 
has been the inability to modify the graphene intrinsic 
doping that is relatively high (> 10^^ cm~^) in most of 
the samples. This fact has prevented these experiments 
to directly image the electronic structure of graphene 
close to the Dirac point. The first scanning probe exper- 
iment on exfoliated graphene on Si02 (Ishigami et al. 



et al.. 2009) have found larger roughness. Even when 



suspended, graphene is never completely flat, and it has 



2007 ) revealed the atomic structure of graphene and the 
nanoscale morphology. The first experiment that was 
able to directly image the electronic structure of exfoli- 
ated graphene close to the Dirac point was performed 



been shown theoretically to possess intrinsic ripples ( Fa- by Martin et al. ( 2008 ) using scanning single-electron 
solino et al. 2007). A local variation of the height pro- transistor, SET, Fig. |32[ The break-up of the density 



file /i(r) can induce a local change of the carrier density landscape in electron-hole puddles as predicted in (Adam 



through different mechanisms, de Juan et al. ( 2007 ) con- 
sidered the change in carrier density due to a local vari- 
ation of the Fermi velocity due to the rippling and found 
that assuming h{r) = ^exp"!""! , a variation of 1 % 
and 10 % in the carrier density was induced for ratios 
A/b of order 0.1 and 0.3 respectively. Brey and Pala- 



et al. 2007 Hwang et al. 2007a) and shown by the DFT 



cios 



(2008) observed that local Fermi velocity changes 
induced by the curvature associated with the ripples in- 
duce charge-inhomogeneities in doped graphene but can- 
not explain the existence of electron-hole puddles in un- 
doped graphene for which the particle-hole symmetry is 
preserved, and then considered the effect on the local 
carrier density of a local variation of the exchange en- 
ergy associated with the local change of the density of 
carbon atoms due to the presence of ripples. They found 
that a modulation of the out-of-plane position of the car- 
bon atoms of the order of 1 — 2 nm over a distance of 
10 — 20 nm induces a modulation in the charge density 
of the order of lO^^cm^. 



LP A ([Polini et al.\ |2008b[ ) and TED theory ( [Rossi and 
Das Sarma[ 2008 ) is clearly visible. The result shown in 
Fig. 32 (a) however does not provide a good quantitative 
characterization of the carrier density distribution due to 
the limited spatial resolution of the imaging technique: 
the diameter of the SET is 100 nm and the distance be- 
tween the SET and the sample is 50 nm and so the spa- 
tial resolution is approximately 150 nm. By analyzing 
the width in density of the incompressible bands in the 



quantum Hall regime, Martin et al. (20081 were able to 



extract the amplitude of the density fluctuations in their 
sample. By fitting with a Gaussian the broadened in- 
compressible bands Martin et al. extracted the value of 
the amplitude of the density fluctuations, identical for all 

and found 



the incompressible bands (Hani et al. 
it to be equal to 2.3 10" — 



2004) 



Kim and Castro Neto 



(2008) 



considered the effect due to the rehybridization of the tt 
and a orbital between nearest neighbor sites. For the 



local shift SEp of the Fermi level, Kim and Castro Neto 



cm~^. Taking this value to 
s using the TFD a corresponding value 
cm~^ is found consistent with typical 
values for the mobility at high density. By calculating 
the ratio between the density fluctuations amplitude ex- 
tracted from the broadening of the incompressible bands 



be equal to rip 
of riimp = 2.4 10 
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FIG. 32 (Color online), (a) Carrier density color map at the 
CNP measured with an SET. The blue regions correspond to 
holes and red regions to electrons. The black contour marks 
the zero density contour. Adapted from Martin et al. (20081. 



(b) Spatial map, on a 80 nmx80 nm region, of the energy shift 
of the CNP in BLG from STM dl/dV map. Adapted from 
Deshpande et al. (2009b 1 . 60 nmx60 nm constant current 
STM topography, (c), and simultaneous dl/dV map, (d) at 
the CNP for SLG , {Vuas = -0.225 VI = 20 pA). Adapted 
from [Zhang et al\ ( |2009b| ) . 



in the Quantum Hall regime and the amplitude extracted 
from the probability distribution of the density extracted 



from the imaging results, Martin et al. (2008) obtained 



the upper bound of 30 nm for the characteristic length of 
the density fluctuations, consistent with the TFD results 



( Rossi and Das Sarma 2008 ) 



An indirect confirmation of the existence of electron- 
hole puddles in exfoliated graphene close to the CNP 
came from the measurement of the magnetic field- 
dependent longitudinal and Hall components of the re- 
sistivity Pxx{H) and pxy{H) (Cho and Fuhrer 2008). 
Close to the Dirac point the measurements show that 
Pxx{H) is strongly enhanced and pxy{H) is suppressed, 
indicating nearly equal electron and hole contributions to 
the transport current. In addition the experimental data 
were found inconsistent with uniformly distributed elec- 
tron and hole concentrations (two-fluid model) but in ex- 
cellent agreement with the presence of inhomogeneously 
distributed electron and hole regions of equal mobility. 

The first STM experiments on exfoliated graphene 
were performed bylZhang et a/.| (|2009b|). The STM ex- 



periments provided the most direct quantitative char- 
acterization of the carrier density distribution of exfo- 
liated graphene 
a 60 



Fig. 32 (c) shows the topography of 

whereas 



60 nm area of exfoliated graphene, 



dl/dV value is directly proportional to the local density 
of states. We can see that there is no correlation between 
topography and dl/dV map. This shows that in current 
exfoliated graphene samples the rippling of graphene, ei- 
ther intrinsic or due to the roughness of the substrate 
surface, are not the dominant cause of the charge den- 
sity inhomogeneities. The dl /dV maps clearly reveal the 
presence of high density regions with characteristic length 
of w 20 nm as predicted by the TFD results. Recently 
more experiments have been performed to directly image 
the electronic structure of both exfoliated single layer 



(Deshpande et al. 2009b). In particular Deshpande et al. 



graphene ( Deshpande et al. 2009a ) and bilayer graphene 



(2009a) starting from the topographic data calculated 



the carrier density fluctuations due to the local curva- 



ture of the graphene layer using Eq. (4.24) and compared 
them to the fluctuations of the dl/dV map. The com- 
parison shows that there is no correspondence between 
the density fluctuations induced by the curvature and 
the ones measured directly. This leads to the conclusion 
that even though the curvature contributes to a variation 
in the electrochemical potential, it is not the main fac- 
tor responsible for the features in the dl/dV map. The 
results for BLG of Fig. 32 (b) show that close to the 



CNP the density inhomogeneities are very strong also in 
BLG and are in semiquantitative agreement with theo- 



retical predictions based on the TF theory (Das Sarma 



et al.\ [2010|. Using an SET [Martin et al.\ ((2009|) have 
imaged the local density of states also in the quantum 
Hall regime. The disorder carrier density landscape has 
also been indirectly observed in imaging experiments of 



coherent transport ( 


Berezovsky et al. 


2010 


and Westervelt 2010 


I- 



D. Transport in the presence of electron-hole puddles 

The previous section showed, both theoretically and 
experimentally, that close to the Dirac point, in the pres- 
ence of long-range disorder, the carrier density landscape 
breaks up in electron hole puddles. In this situation the 
transport problem becomes the problem of calculating 
transport properties of a system with strong density in- 
homogeneities. The first step is to calculate the conduc- 
tance of the puddles Gp and PNJ, Gpnj. 

We have Gp = Ta where F is a form- factor of order one 
and a is the puddle conductivity. Away from the Dirac 



point (see Sec. HI. A ) the RPA-Boltzmann transport the- 
ory for graphene in the presence of random charge impu- 
rities is accurate. From the RPA-Boltzmann theory we 
have a — e\{n)\p{{n),ni^p,rs, d,T). For the purposes of 
this section it is convenient to explicitly write the depen- 



dence of a on 



''imp 



by introducing the function 



Fig. 32 (d) shows the dl/dV map of the same area. The 



Firs,d,T)^^^^=2nn, 
2e 



Vf 



(4.25) 
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so that we can write: 



2e2|(n)| 



F(r„d,T). 



(4.26) 



Expressions for F{rs, d,T) at T = (or its inverse) were 
originally given in ( |Adam et al] [2007| (see Eq. ( [3^ ). 
We can define a local spatially varying "puddle" conduc- 
tivity (j{r) if a{r) varies on length scales that are larger 
than the mean free path I, i.e. 



Vcr(r) 



a(r) 



> I. 



(4.27) 



By substituting (n) with n(r) let us use Eq. (4.26) to 
define and calculate the local conductivity: 



a(r) 



2e^ n(r) 
h ni„ 



F(r„d,T). 



(4.28) 



''imp 

Considering that I = ha/{2e'^kp) and using Eq. (4.28) 
then inequality (4.27) takes the form: 



Vn(r) 



n(r) 



> 



Firs,d,T) 



(4.29) 



''imp 



As shown in the previous sections, at the CNP most of the 
graphene area is occupied by large electron-hole puddles 
with size of the order of the sample size L and density of 
the order of rii-ms ~ "-imp- For graphene on Si02 we have 
Tg = 0.8 for which is F = 10. Using these facts we find 



that inequality (4.291 is satisfied when 
F{r,,d,T) 



1 



(4.30) 



mip 



i.e. when the sample is much larger than the typical 
in-plane distance between charge impurities. Consider- 
ing that in experiments on bulk graphene L > 1 /im 



and Tiiyyip G 



[10 



10 



10 ] cm we see that inequality 



(4.30) is satisfied. In this discussion we have neglected 
the presence of the small regions of high density and size 
^. For these regions the inequality (4.301 is not satis- 
fied. However these regions, because of their high carrier 
density, steep carrier density gradients at the boundaries 
and small size ^ < I, are practically transparent to the 
current carrying quasiparticles and therefore, given that 
they occupy a small area fraction and are isolated (i.e. 
do not form a path spanning the whole sample), give a 
negligible contribution to the graphene resistivity. This 



fact, and the vahdity of inequality (4.30) for the large 
puddles ensure that the local conductivity cr(r) as given 
by Eq. (4.28) is well defined. In the limit rg ^ 1 one can 



use the analytical results for the density distribution to 



reach the same conclusion ( Fogler 2009 ) 



To calculate the conductance across the PNJ, quan- 
tum effects must be taken into account. In particular, 
as seen in Sec. |II.B] for Dirac fermions, we have the phe- 
nomenon of Klein tunneling (Dombey and Calogeracos 



1999 Klein 1929) i.e. the property of perfect transmis- 
sion through a steep potential barrier perpendicular to 
the direction of motion. The PNJ conductance, GpNj, 



2006a 


) and 


Zhang et al. 


(2008 



mate the steepness of the electrostatic barrier at the PNJ, 
i.e. the ratio between the length scale, D, over which the 
screened potential varies across the PNJ and the Fermi 
wavelength of the carriers at the side of the PNJ. From 
the TED theory we have that at the sides of the PNJ 
n ss TT-i-ms so that kp — y^TrnT^, and that D « 1/kp so 
that kpD M I. In this limit, the conductance per unit 



length of a PNJ, 5pnj is given by (Cheianov and Fal'ko 



2006b Zhang et al. 2008): 



5PNJ 



h 



kp. 



(4.31) 



so that the total conductance across the boundaries of 
the electron-hole puddles is GpNj = P<7pnj with p the 



perimeter of the typical puddle ( Fogler 2009 Rossi et al. 



2009 ) . Because the puddles have size comparable to the 
sample size, p ^ for typical experimental conditions 
(L > 1 ^m and n™p € [10^° - lO^^] cm^^) ^gj^^g 



and (4.31 ) we find 



;p»Gp = rF(r„d,T) 



2e^ 
h 



Frmsl 



imp 

(4.32) 

^ 1 the inequality 
(4.32) is vahd for any value of nimp (Fogler 2009). In- 



i.e. GpNj ^ Gp. In the limit 



equality (4.32) shows that in exfoliated graphene sam- 
ples, transport close to the Dirac point is not percola- 
tive: the dominant contribution to the electric resistance 
is due to scattering events inside the puddles and not 



to the resistance of the puddle boundaries (Fogler 2009 



Rossi et al. 2009). This conclusion is consistent with the 
results of Adam et al. (2009a) in which the graphene 



conductivity in the presence of Gaussian disorder ob- 
tained using a full quantum mechanical calculation was 
found to be in agreement with the semiclassical Boltz- 
mann theory even at zero doping provided the disorder 
is strong enough. Given (i) the random position of the 
electron-hole puddles, (ii) the fact that because of in- 



equality (4.30) the local conductivity is well defined, and 



(in) the fact that Gpnj ^ Gp, the effective medium the- 



ory, EMT, (Bruggeman 1935 Hori and Yonezawa 1975 



Landauer 1952) can be used to calculate the electrical 
conductivity of graphene. The problem of the minimum 
conductivity at the CNP can be expressed as the problem 
of correctly averaging the individual puddle conductivity. 



Using Eq. (4.28), given a carrier density distribution, the 
conductivity landscape can be calculated. In the EMT 
an effective medium with homogeneous transport prop- 
erties equivalent to the bulk transport properties of the 
inhomogeneous medium is introduced. Starting from the 
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local relation between current J and electric potential V: 



(a) 



J(r) = -a{r)WV{r) 



(4.33) 



the effective-medium conductivity, ctemt, is defined 
through the equation: 



(J(r)) = -aEMT(VV^(r)) 



(4.34) 



where the angle bracket denotes spatial and disorder av- 



erages. Equations (4.33), (4.34) along with the condition 
that in the effective medium the electric field —WV is uni- 
form are sufficient to calculate (Temt- The derivation of 



the relation between (Temt and a{r) using (4.33), (4.341, 
requires the solution of the electrostatic problem in which 
a homogeneous region of conductivity a(r) is embedded 
in an infinite medium of conductivity (Teut- When the 
shape of the homogeneous regions, puddles, in the real 
medium is random the shape of the homogeneous regions 
used to derive the expression of (Jemt is unimportant, 
and they can be assumed to be spheres having the same 
volume as each puddle. (Bruggeman 1935 Landauer 



1952). For a 2D system the solution of the electrostatic 



problem gives (Bruggeman 1935 Landauer 1952) 



cr(r) + (Temt 



0. 



(4.35) 



Equation (4.35) can also be viewed as an approximate 
resummation of the infinite diagrammatic series for the 
macroscopic a using the self-consistent single-site approx- 
imation (iHori and Yonezawa 1975). Disorder averaging 
Eq. (|4.35| we find: 



f^EMT 



> ^ da 



P(fT) =0, 

a(r) + aEMT / J (t + ctemt 

(4.36) 

where P(cr) is the probability for the local value of cr. 
Using the relation between the local value a and the local 



value of the carrier density n, Eq. (4.281, Eq. (4.361 can 
be rewritten in the form: 



dn44^^PW=0. (4.37) 
a[n) + (Temt 



where P[n\ is the density probability distribution that 
can be calculated using the TFD theory. Fig. [311 Using 



the TFD results and equation (4.37) the conductivity at 
the Dirac point and its vicinity can be calculated. 

Fig.j^shows a{Vg) as obtained using the TFD-fEMT 
theory ^fiossi et aL[|2009[ ). The theory correctly predicts 
a finite value of a very close to the one measured ex- 
perimentally. At high gate voltages, the theory predicts 
the linear scaling of cr as a function of Vg. The theory 
correctly describes the crossover of a from its minimum 
at Vg — to its linear behavior at high gate voltages. 



Fig. 33 (a) also shows the importance of the exchange- 
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FIG. 33 (Color online). Solid (dashed) lines show the EMT- 
TFD results with (without) exchange, (a) g as function of 
Vg for three different values of nin 



the results presented in Ref. Adam et al. ( 2007 1 for the same 
values of nimp- (b) Minimum conductivity as a function of 
interaction parameter. Solid squares show the experimental 
results of pang et al. (20081. (c) amin as a function of Uimp, 
ra — 0.8 and d = 1 nm. For comparison the results obtained 
in Ref. |Adam et al.\ (|2007|) are also shown, dot-dashed line. 

1 and Uimp = 10^^ cm~^. 
(20(39[); (c) adapted 



,m et al. (20071 are also shown, 
(d) Uniin as a function of rs for d ■ 



[a,), (b), (d) Adapted fr om Rossi et al 
from Chen et al. 



(2008a I 



Cmin on Uiaip is shown in Fig. 33 In panel (a) the theoret 



ical results are shown, crmin increases as riimp decreases, 
the dependence of a^in on nimp is weaker if exchange- 
correlation terms are taken into account. Panel (c) of 
Fig. [33] shows the dependence of (Jmin on the inverse mo- 
bifity 1/^ cx nimp 



. In 



as measured by Chen et al. ( 2008a ) 
this experiment the amount of charge impurities is con- 
trolled by potassium doping. Fig. 33 (b) shows the re- 
sults for cr„iin as a function of r^. The solid (dashed) line 
shows the calculated values of o-min including (neglect- 
ing) exchange. fTmin has a non-monotonic behavior due 
to the fact that affects both the carrier density spatial 
distribution by controlling the strength of the disorder 
potential, screening and exchange as well as the scat- 
tering time T. The dependence of (Tmin on Vg has been 



measured in two recent experiments (Jang et al. 2008 



Ponomarenko et al. 2009). In these experiments the fine 
structure constant of graphene, r^, is modified by placing 
the graphene on substrates with different k and/or by us- 
ing materials with k 7^ 1 as top dielectric layers. In [Jang] 
et al. (2008) graphene was placed on Si02 and Vs was 



reduced from 0.8 (no top dielectric layer) to 0.56 by plac- 
ing ice in vacuum as a top dielectric layer. The resulting 
change of cTmin is shown in Fig. [33] (b) . As predicted by 
the theory, when V^c is included, for this range of values 
of r^, CTmin is Unaffected by the variation of r^. Overall 



correlation term at low gate voltages. The dependence of the results presented in (Jang et al. 2008) are consistent 
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with charge impurity being the main source of scatter- 
ing in graphene. In Ponomarenko et al. (20091, rs was 



varied by placing graphene on substrates with different 
dielectric constants and by using glycerol, ethanol and 



water as a top dielectric layer. Ponomarenko et al. ( 2009 ) 
found very minor differences in the transport properties 
of graphene with different dielectric layers, thus conclud- 
ing that charge impurities are not the dominant source of 
scattering. At the time of this writing the reasons for the 



discrepancy between the results of Jang et al. { 2008 ) and 



Ponomarenko et al. ( 2009 ) are not well understood. The 
experiments are quite different. It must be noted that 
changing the substrate and the top dielectric layer, in 
addition to modifying rg, is likely to modify the amount 
of disorder seen by the carriers in the graphene layer. By 
not modifying the substrate and by placing ice in ultra- 



high- vacuum, Jang et al. ( 2008 1 minimized the change of 



disorder induced by modifying the top dielectric layer. 

The approach presented above based on TFD and 
EMT theories can be used to calculate other transport 
properties of SLG and BLG close to the Dirac point. 



In Hwang et al. (20091 the same approach was used to 



calculate the thermopower of SLG with results in good 



agreement with experiments 


Checkelsky et al. 


2009 


Wei 


et a/.r2009; Zuev et al. 


2009 


. In 


Das Sarma et al. 


(2010) 



the TFD-I-EMT approach was used to calculate the elec- 



trical conductivity in BLG. In Tiwari and Stroud (2009) 

the EMT based on a simple two-fluid model was used to 
calculate the magnetoresistance at low magnetic fields of 
SLG close to the CNP. 

It is interesting to consider the case when Gp » Gpn j • 
This limit is relevant for example when a gap in the 
graphene spectrum is opened. The limit Gp ^ GpNj 
was studied theoretically in ( Cheianov et aZ.[|2007b ) by 
considering a random resistor network (RRN) model on 
a square lattice in which only nearest neighbors and next 
nearest neighbor are connected directly. Mathematically 
the model is expressed by the following equations for 
the conductance between the sites and with 
\^-^'\ < 1, \j-f\ < 1: 



G 



j) 

j+i) 



9[l + i-^y+^v^,,] A 



g [1 - i-iy+'v^,,] /2; 

7 < 1. 



G 



19, 



where rjij is a random variable, 



(4.38) 
(4.39) 
(4.40) 

(4.41) 



and p is proportional to the doping (n). For 7 — 0,p = 
we have percolation. Finite p and 7 are relevant per- 
tirbations for the percolation leading to a finite corre- 
lation length ^ (p, 7) . On scales much bigger than ^ 
the RRN is not critical, and consists of independent re- 
gions of size ^ so that a is well defined with scaling 



(Cheianov et al. 



2007bl ) a(p,7) = [a/aP,l)rg with 
p* = 7'^/" and V = 4/3, 
: 7/4 and x sa 0.97 are re- 



f(p,7) ^ aj-yPip/p*) 
fi — l/{h + x) where h 
spectively the fractal dimension of the boundaries be- 
tween the electron-hole puddles and the conductance ex- 
ponent, (G(L)) = {a/LYg, at the percolation threshold 
(iCheianov et aLl |2007bl llsichenkol 119921). Fig. Isll shows 



the results obtained solving numerically the RRN defined 
by Eq. |4.40[ The numerical results are well fitted using 
for F{p/p*) the function F{z) = (l +z^)('^/^). Estimating 
" \ akF and 73 ^ j- ^akpY^^ (Cheianov and Fal'ko 



2006b), Cheianov et al. (2007b) estimate that 



n 

From the TFD results one gets 6n 



''imp 



(4.42) 
and so 

Eq. (4.42) predicts that that fTmm should increase with 
f^imp, a trend that is not observed in experiments. The 
reason for the discrepnacy is due to the fact that for cur- 
rent experiments the relevant regime is expected to be 
the one for which GpNj 3> Gp. The minimum conduc- 
tivity was also calculated for bilayer graphene ("Adam and 

2010) 



Das Sarma 2008a Das Sarma et al. 



Other works 

calculated a^n using different models and approximations 



for regimes less relevant for current experiments (Adam 



Trushin and Schliemann 



2007) 



and Das Sarma, ,2008a, ,Cserti 


2007 


Cserti et al.\ 


2007 


Groth et al. 2008 


Katsnelson 


20061 


Tushin et al.\ 


2010 



Although the semiclassical approach presented in this 
section is justified for most of the current experimental 
conditions for exfoliated graphene its precise range of va- 
lidity and level of accuracy close to the Dirac point can 
only be determined by a full quantum transport calcula- 
tion that takes into account the presence of charge impu- 
rities. This is still an active area of research and work is 
in progress to obtain the transport properties of graphene 



using a full quantum transport treatment (Rossi et al. 



2010b) 
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FIG. 34 Collapse of the conductivity data obtained for RNNs 
with various — l/2<p<l/2 and values of the parameter 7 
onto a single curve. Points: numerical results; line: best fit 
obtai ned using the equations presented in the text. Adapted 
from ICheianov et al.\ (|2007b I 
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V. QUANTUM HALL EFFECTS 
A. Monolayer graphene 

1. Integer quantum Hall effect 

The unique properties of the quantum Hall effect in 
graphene are among the most striking consequences of 
the Dirac nature of the massless low energy fermionic 
excitations in graphene. In the presence of a perpendic- 
ular magnetic field, B, electrons (holes) confined in two 
dimensions are constrained to move in close cyclotron or- 
bits that in quantum mechanics are quantized. The quan- 
tization of the cyclotron orbits is reflected in the quan- 
tization of the energy levels: at finite B the B = dis- 
persion is replaced by a discrete set of energy levels, the 
Landau levels (LL). For any LL there are iV^ — BA/(j)o 
degenerate orbital states, where A is the area of the sam- 
ple and is the magnetic quantum flux. Quantum Hall 



chiral nature of the quasiparticles in graphene. The fac- 



effects ( 


Das Sarma and Pinczuk 


1996 


MacDonald 


1990 


Prange and Girvin 


1990 


) appear when TV is comparable 



to the total number of quasiparticles present in the sys- 
tem. In the quantum Hall regime the Hall conductivity 
axy exhibits well developed plateaus as a function of car- 
rier density (or correspondingly magnetic field) at which 
it takes quantized values. At the same time, for the range 
of densities for which axy is quantized, the longitudinal 



conductivity axx is zero (Halperin, 1982 Laughlin 1981 1. 



For standard parabolic 2DEG, (such as the ones created 
in GaAs and Si quantum wells), the LL have energies en- 
ergy hujcin + 1/2) where n = 0, 1, 2, ... and ojc = eB/mc, 
m being the effective mass, is the cyclotron frequency. 
Because the low energy fermions in graphene are mass- 
less it is immediately obvious that for graphene we cannot 
apply the results valid for standard 2DEG (wc would ap- 
pear to be infinite). In order to find the energy levels, 
En, for the LL the 2D Dirac equation must be solved in 



the presence of a magnetic field ( Gusynin and Sharapov 



2005 


Haldane 


1988 


Jackiw 


1984| Peres et al. 


2006 



The result is given by Eq. (1.13) (a). Differently from 
parabolic 2DEG, in graphene we have a LL at zero en- 
ergy. In addition we have the unconventional Hall quan- 



tization rule for axy ( 


Gusynin and Sharapov 


2005 


Peres 


et al. 2006 


Zheng and Ando 2002 


I: 



ax 



e 
~h 



compared to the one valid for regular 2DEGs: 

c2 



gn 



h 



(5.1) 



(5.2) 



shown in Fig. 35 where g is the spin and valley degen- 
eracy. Because in graphene the band dispersion has two 
inequivalent valleys valleys (7 = 4 (for GaAs quantum 
wells we only have the spin degeneracy so that g = 2). 



tor 1/2 in (5.1) can be understood as the term induced 
by the additional Berry phase that the electrons, due to 
their chiral nature, acquire when completing a close orbit 



(Luk'yanchuk and Kopelevich 2004 Mikitik and Shar- 



lai 1999). Another way to understand its presence is 
by considering the analogy to the relativistic Dirac equa- 
tion (iGeim and MacDonaldl 120071 lYangl 120071). From 



this equation two main predictions ensue: (i) the elec- 
trons have spin 1/2, (ii) the magnetic g— factor is exactly 
equal to 2 for the "spin" in the non-relativistic limit. As a 
consequence the Zeeman splitting is exactly equal to the 
orbital-splitting. In graphene the pseudospin plays the 
role of the spin and instead of Zeeman splitting we have 
"pseudospin-splitting" but the same holds true: the pseu- 
dospin splitting is exactly equal to the orbital splitting. 
As a consequence the n — th LL can be thought as com- 
posed of the degenerate pseudospin-up states of LL n and 
the pseudospin-down states of LL n — 1. For zero mass 
Dirac fermions the first LL in the conduction band and 
the highest LL in the valence band merge contributing 
equally to the joint level at E = resulting in the half- 



odd-integer quantum Hall effect described by Eq. (5.1 1. 
For the E — LL, because half of the degenerate states 
are already filled by hole-like (electron-like) particles, we 
only need 1/2N^ electron- like (hole- like) particles to fill 
the level. 
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FIG. 35 (Color online) Illustration of the integer QHE found 
in 2D semiconductor systems (a), incorporated from MacDon- 
aldl (pOOf; [Prange and G"Ir^([l990|), SLG (b), BLG (c) The 



sequences of Landau levels as a function of carrier concentra- 
tions n are shown as dark and light peaks for electrons and 
holes, respectively. Adapted from Novoselov et al. (20061. 
(d) a^y (red) and p^x (blue) of SLG as a function of carrier 
density measured experimentally at T = 4K and B = 14T. 



Adapted from Novoselov et al. { 2005a I 



The remarkable quantization rule for axy has been ob- 



The additional 1/2 in Eq. (5.1) is the hallmark of the et al. 2005) as shown in Fig. 35 (d). The experimental 



served experimentally (Novoselov et al. 2005a Zhang 
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observation of Eq. ( 5.1 ) shows clearly the chiral nature of 
the massless quasiparticles in graphene. There is another 
important experimental consequence of the Dirac nature 
of the fermions in graph ene. Because in graphene i?„ 
scales as ^/nB (Eq. (1.13) (a)) rather than linearly as in 
regular 2DEG (Eq. (1.131 (c)) at low energies (n) the en- 
ergy spacing A„ = En+i — En between LL can be rather 
large. Because the observation of the quantization of a^y 
relies on the condition A„ ^ kgT {T being the temper- 
ature) it follows that in graphene the quantization of the 
LL should be observable at temperatures higher than in 
regular parabolic 2DEG. This fact has been confirmed 
by the observation in graphene of the QH effect at room 



temperature (Novoselov et al.i 20071 Graphene is the only 



known material whose quantum Hall effect has been ob- 
served at ambient temperature (albeit at high magnetic 
fields). 

By applying a top gate, p-n junctions, PNJ, can be cre- 
ated in graphene. In the presence of strong perpendicular 
fields graphene PNJ exhibit unusual fractional plateaus 
for the conductance that hav e be en studied experimen- 
tally bylWilhams et al.\t007\ andlOzyilmaz et al.H2007\ 



and theoretically by Abanin and Levitov ( 2007 1 . Numer- 



ical studies in the presence of disorder have been per- 



Low (2009) 



formed by Long et al. (2008), Li and Shen (2008), and 



2. Broken-symmetry states 



The sequence of plateaus for given by Eq. (5.1 1 de- 
scribes the QH effect due to fully occupied Landau levels 
including the spin and valley degeneracy. In graphene 
for the fully occupied LL we have the "filling factors" 
ly = gN/N^ = 4(n -h 1/2) = ±2, ±6, ±10, .... In this sec- 
tion we study the situation in which the spin or valley, or 
both, degeneracies are lifted. In this situation QH effects 
are observable at intermediate filling factors 1^ = 0, ±1 
for the lowest LL and v — ±3, ±4, ±5 for n ~ ±1 LL. 
The difficulty in observing these intermediate QH effects 
is the lower value of the energy gap between successive 
splitted LL. If the gap between successive Landau levels 
is comparable or smaller than the disorder strength, the 
disorder mixes adjacent LL preventing the formation of 
well defined QH plateaus for the Hall conductivity. For 
the most part of this section we neglect the Zeeman cou- 
pling that turns out to be the lowest energy scale in most 
of the experimentally relevant conditions. 



Koshino and Ando (2007) showed that randomness in 



the bond couplings and on-site potential can lift the val- 
ley degeneracy and cause the appearance of intermediate 



Landau Levels. Fuchs and Lederer (2007) considered the 



electron-phonon coupling as the possible mechanism for 
the lifting of the degeneracy. However in most of the 
theories the spin/valley degeneracy is lifted due to inter- 



2006 


Ezawa 2007: '2008| |Goerbig et al. 


'2006'; 'Gusyninl 


et al. 


2006 Herbut,|2007 Nomura and MacDonald„2006| 


Yang et al. 2006), in particular, electron- electron inter- 



actions. When electron-electron interactions are taken 
into account, the quasiparticles filling a LL can polarize 
in order to minimize the exchange energy (maximize it in 
absolute value). In this case, given the S'J7(4) invariance 
of the Hamiltonian, the states 



(5.3) 



i<i<j\/ fc 



where i is the index of the internal states that runs from 
1 to M = V - 4(n - 1/2) < 4, and |0) is the vacuum, 
are exact eigenstates of the Hamiltonian. For a broad 
class of repulsive interaction \'^q) is expected to be the 



exact ground state (Yang, 2007 Yang et al. 2006). The 



state described j^'g) is a "ferromagnet" , sometimes called 
a QH ferromagnet, in which either the real spin or the 
pseudospin associated with the valley degree of freedom 
is polarized. The problem of broken symmetry states 
in the QH regime of graphene is analogous to the prob- 
lem of "quantum Hall ferromagnetism" studied in regu- 
lar 2DEG in which, however, normally only the SU{2) 
symmetry associated with the spin can be spontaneously 
broken (notice however that for silicon quantum wells the 
valley degeneracy is also present so that in this case the 
Hamiltonian is SU{N) {N > 2) symmetric). Because in 
the QH regime the kinetic energy is completely quenched, 
the formation of polarized states depends on the rela- 
tive strength of interaction and disorder. For graphene. 



Nomura and MacDonald ( 2006 ) , using the Hartree-Fock 



approximation derived a "Stoner criterion" for the exis- 
tence of polarized states, i.e. QH ferromagnetism, for a 
given strength of the disorder. |Chakraborty and Pieti-| 
lainen (2007) have numerically verified that QH ferro- 



magnetic states with large gaps are realized in graphene. 



Sheng et al. (2007) using exact diagonalization studied 
the interplay of long-range Coulomb interaction and lat- 
tice effects in determining the robustness of the z/ = ±1 
and u = ±3 states with respect to disorder. Nomura 



et al. ( 2008 ) studied the effect of strong long-range disor- 



der. Wang et al. (2008) have performed numerical studies 
that show that various CDW phases can be realized in 
the partially-filled v = ±3 LL. 

Experimentally the existence of broken-symmetry 
states has been verified by Zhang et al. (2006), 



Fig. 36 (a), that showed the existence of the interme- 

0,±1 for the n = LL 
= ±4 for the n = 1 LL 



action effects (Abanin et al. 2007b Alicea and Fisher 



diate Landau levels with = 
and the intermediate level v 
that is therefore only partially resolved. Given that the 
magnetic field by itself does not lift the valley degener- 
acy, interaction effects are likely the cause for the full 
resolution of the n = LL. On the other hand a careful 
analysis of the data as a function of the tilting angle of 
the magnetic field suggests that the partial resolution of 
the n = 1 LL is due to Zeeman splitting ( Zhang et al. 
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2006) 
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FIG. 36 (Color online) (a) tj^y, as a function of gate voltage 
at different magnetic fields: 9 T (circle), 25 T (square), 30 T 
(diamond), 37 T (up triangle), 42 T (down triangle), and 45 
T (star). All the data sets are taken at T=f .4 K, except for 
the B=9 T curve, which is taken at T=30 mK. Left upper 
inset: R^x and R^y for the same device measured at B=25 T. 
Right inset: detailed a^y data near the Dirac point for B=9 
T (circle), 11.5 T (p entagon), and 17.5 T (hexagon) at T=30 



mK. Adapted from Zhang et al 
resistance 7?; 



(20061 



(b) Longitudinal 
as a function of gate voltage Vg = — Vo at 
0.3 K and several values of the magnetic field: 8, 11 and 14 T. 
The inset shows in false color a graphene crystal (dark red) 
with An leads deposited (yellow regions). The bar indicates 
5 iim. At Vg — 0, the peak in R^^ gro ws to 190 kfJ at 14 T. 



Adapted from Checkelsky et al. ( 2008 1 



Using Eq. (5.4), a possible resolution of the v = Q 



"anomaly" is obvious: for any finite Oxy the vanishing 
of a XX corresponds to the vanishing of pxx^ however for 
'^xy = we have Pxx — ^I'^xx so that Oxx ~^ implies 
Pxx oo- This is very similar to the Hall insulator phase 
in ordinary 2D parabolic band electron gases. This sim- 
ple argument shows that the fact that pxx seem to diverge 
for T — for the v — Q state is not surprising. How- 
ever this argument may not be enough to explain the 
details of the dependence of Pxxi^^ = 0) on temperature 



and magnetic field. In particular Checkelsky et al. (2009) 



found evidence for a field-induced transition to a strongly 
insulating state at a finite value of B. These observa- 
tions suggest that the v = ground state might differ 



from the 5J7(4) eigenstates (5.3 1 and theoretical calcu- 
lations proposed that it could be a spin-density-wave or 



charge-density- wave (Herbut 2007 Jung and MacDon- 



aid, 2009). It has also been argued that the divergence 
of Pxx{v = 0) might be the signature of Kekule instability 



(Hou et al. 2010 Nomura et al. 2009) 



Giesbers et al. (2009) have interpreted their experi- 



mental data using a simple model involving the open- 



ing of a field dependent spin gap. Zhang et al. (2009a 



have observed a cusp in the longitudinal resistance pxx 
for v 1/2 and have interpreted this as the signature of 
a transition from a Hall insulating state for v > 1/2 to 
a collective insulator, like a Wigner crystal (|Zhang and 
Joglekar 2007), for < 1/2. No consensus has been 



reached so far, and much more work is needed to under- 
stand the v = state in graphene. 



3. The V = Q state 



In the previous section we have seen that in strong 
magnetic fields the lowest LL can be completely resolved 
and the spin and valley degeneracies may be lifted. In 
particular, an approximate plateau for axy appear for 
V = Q. This state has been experimentally studied in a 



series of works by Checkelsky et al. (2008 2009); Gies- 



bers et aL| ( |2009[ ); |Jiang et al.\\2mi\ , see Fig.|36|(b). The 
state is unique in that the plateau of Uxy corresponds to 
a maximum of the longitudinal resistivity pxx in con- 
trast to what happens for z/ 7^ where a plateau of Gxy 
corresponds to zero longitudinal resistivity. In addition, 
the f = edge states are supposed to not carry any 



charge current, but only spin currents (Abanin et al. 
2006 2007a|c ). As pointed out by Das Sarma and Yang 
( 2009 ) however the situation is not surprising if we re- 
call the relations between the resistivity tensor and the 
conductivity tensor: 



Pxx 



Pxy 



'xy 



(5.4) 



and the fact that the quantization of <Jxy is associated 
with the vanishing of Gxx ■ This can be seen from Laugh- 



lin's gauge argument (Halperin 1982 Laughlin 1981). 



4. Fractional quantum Hall effect 



In addition to QH ferromagnetism, the electron- 
electron interaction is responsible for the fractional quan- 
tum Hall effect, FQHE. For the FQHE the energy gaps 
are even smaller than for the QH ferromagnetic states. 
For graphene the FQHE gaps have been calculated by 
Apalkov and Chakraborty|(|2006|);|Toke et al.\\2QQQ\. For 



the V — l/i the gap has been estimated to be of the 
order of Q.Qbe^ / kIb, where Ib = {hc/eBY/'^ is the mag- 
netic length. Because of the small size of the gaps the 
experimental observation of FQHE requires high quality 
samples. For graphene very low amount of disorder can 
be achieved in suspended samples and in these suspended 



samples two groups (Bolotin et al. 2009 Du et al.\ 2009) 



have recently observed signatures of the = 1/3 frac- 
tional quantum Hall state in two-terminal measurements, 
see Fig. |37] A great deal of work remains to be done in 
graphene FQHE. 
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(a) 



FIG. 37 (Color online) Graphene fractional quantum Hall 
data taken from (a) Du et al. { 2009 1 and (b) Bolotin et al 



(20091, observed on two probe suspended graphene samples. 



B. Bilayer graphene 

1. Integer quantum Hall effect 

In bilayer graphene the low energy fermionic excita- 
tions are massive, i.e. with good approximation the 
bands are parabolic. This fact would suggest that the 
bilayer QH effect in graphene might be similar to the one 
observed in regular parabolic 2DEG. There are however 
two important differences: the band structure of bilayer 
graphene is gapless and the fermions in BLG, as in SLG, 
are also chiral but with a Berry phase equal to 2tt instead 
of TT (iMcCann and Fal'kol 120061). As a consequence, as 



shown in Eq. (1.13) (c), the energy levels have a dif- 
ferent sequence from both regular 2DEGs and SLG. In 
particular BLG also has a LL at zero energy, however, 
because the Berry phase associated with the chiral na- 
ture of the quasiparticles in BLG is 2-n, the step between 
the plateaus of (Jxy across the CNP is twice as big as in 
SLG (as shown schematically in Fig. 35) (c). One way to 
understand the step across the GNP is to consider that 
in BLG the n = and n = 1 orbital LL are degener- 
ate. The spin and valley degeneracy factor, g, in BLG is 
equal to 4 as in MLG. In BLG the valley degree of free- 
dom can also be regarded as a layer degree of freedom 
considering that without loss of generality we can use a 
pseudospin representation in which the K valley states 
are localized in the top layer and the K' states in the 
bottom layer. The QH effect has been measured exper- 
imentally. Fig. 38 (a) shows the original data obtained 
by Novoselov et al. ( 2006 1 . In agreement with the theory 
the data show a double size step, compared to SLG, for 
a^y across the CNP. 



2. Broken-symmetry states 

As discussed in the previous section, the = LL 
in BLG has an 8-fold degeneracy due to spin degeneracy. 
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FIG. 38 (Color online) (a) Measured Hall conductivity Oxy 
in BLG as a function of carrier density for i? = 12 T, in blue, 
and B = 20 T, in red, at r=4 K. (b) Measured longitudinal 
resistivity in BLG at r=4 K and B = 12 T. The inset shows 
the calculated BLG bands close to the CNP. Adapted from 
Novoselov et aZ.|(2006). (c) Two terminal conductance, G, as 



a function of carrier density at T^lOO mK for different values 
of the magnetic field in suspended BLG. Broken symmetry 
states in BLG. Adapted from |Feldman et aT] ( |2009| ). 



valley (layer) degeneracy and n = 0, n = 1 orbital LL 
degeneracy. The En ^ LL have only a 4-fold degener- 
acy due to spin and valley degeneracy. As discussed for 
SLG, it is natural to expect that the degeneracy of the 
full LL will be lifted by external perturbations and/or 
interactions. Similar considerations to the ones made 
in section |V.A.2| for SLG apply here: the splitting can 



be due to Zeeman effect (Giesbers et al. 2009), strain- 



induced lifting of valley degeneracy ( Abanin et al. 2007b ) 



or Coulomb interactions (Barlas et al. 2008 Ezawa 



2007). Barlas et al. (2008) considered the splitting of the 
En, = LL in BLG due to electron-electron interactions 
and calculated the corresponding charge-gaps and filling 
sequence. As in SLG, the charge gaps of the splitted LLs 
will be smaller than the charge gap, hujc, for the fully 
occupied LLs and so the observation of QH plateaus due 
to the resolution of the LL requires higher quality sam- 
ples. This has recently been achieved in suspended BLG 



samples (Feldman et al. 2009 Zhao et al.\ 20101 in which 



the full resolution of the eightfold degeneracy of the zero 
energy LL has been observed. Fig. [38] (b). By analyzing 
the dependence of the maximum resistance at the CNP 



on B and T Feldman et al. ( 2009 ) concluded that the the 



observed splitting of the ii^„ = LL cannot be attributed 
to Zeeman effect. Moreover the order in magnetic fields 
in which the broken-symmetry states appear is consistent 



with the theoretical predictions of (Barlas et al. 2008). 
These facts suggest that in BLG the resolution of the 
octet zero energy LL is due to electron-electron interac- 
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tions. 



VI. CONCLUSION AND SUMMARY 

In roughly five years, research in graphene physics 
has made spectacular advance starting from the fabri- 
cation of gated variable-density 2D graphene monolay- 
ers to the observations of fractional quantum Hall effect 
and Klein tunneling. The massless chiral Dirac spectrum 
leads to novel integer quantum Hall effect in graphene 
with the existence of a n = quantized Landau level 
shared equally between electrons and holes. The non- 
existence of a gap in the graphene carrier dispersion leads 
to a direct transition between electron-like metallic trans- 
port to hole-like metallic transport as the gate voltage is 
tuned through the charge neutral Dirac point. By con- 
trast, 2D semiconductors invariably become insulating at 
low enough carrier densities. In MLG nanoribbons and 
in BLG structures in the presence of an electric field, 
graphene carrier transport manifests a transport gap be- 
cause there is an intrinsic spectral gap induced by the 
confinement and the bias field, respectively. The precise 
relationship between the transport and the spectral gap 
is, however, not well understood at this stage and is a 
subject of much current activity. Since back-scattering 
processes are suppressed, graphene exhibits weak anti- 
localization behavior in contrast to the weak localization 
behavior of ordinary 2D systems. The presence of any 
short-range scattering, however, introduces inter-valley 
coupling, which leads to the eventual restoration of weak 
localization. Since short-range scattering, arising from 
lattice point defects, is weak in graphene, the weak anti- 
localization behavior is expected to cross over to weak 
localization behavior only at very low temperatures al- 
though a direct experimental observation of such a local- 
ization crossover is still lacking and may be difficult. 

The observed sequence of graphene integer quantized 
Hall conductance follows the expected formula, Uxy — 
{Ae'^ /h){n+l/2), indicating the Berry phase contribution 
and the n — Landau level shared between electrons and 
holes. For example, the complete lifting of spin and valley 
splitting leads to the observation of the following quan- 
tized Hall conductance sequence i/ = 0, ±1,±2, ... with 
Cxy = ve^ /h whereas in the presence of spin and valley 
degeneracy (i.e., with the factor of 4 in the front) one gets 
the sequence v = ±2, ±6,... The precise nature of the 
V — Q IQHE, which seems to manifest a highly resistive 
{Pxx ~^ oo) state in some experiments but not in others, 
is still an open question as is the issue of the physical 
mechanism or the quantum phase transition associated 
with the possible spontaneous symmetry breaking that 
leads to the lifting of the degeneracy. We do mention, 
however, that similar, but not identical, physics arises in 
the context of ordinary IQHE in 2D semiconductor struc- 
tures. For example, the 4-fold spin and valley degeneracy. 



partially lifted by the applied magnetic field, occurs in 2D 
Si-(IOO) based QHE, as was already apparent in the orig- 



inal discovery of IQHE by Klaus von Klitzing (Klitzing 



et al. 19801. The issue of spin and valley degeneracy 



lifting in the QHE phenomena is thus generic to both 
graphene and 2D semiconductor systems although the 



origin of valley degeneracy is qualitatively different ( Eng 



et al. 2007 McFarland et al. 2009) in the two cases 



The other similarity between graphene and 2DEG QHE 
is that both systems tend to manifest strongly insulating 
phases at very high magnetic field when v <^ 1. In semi- 
conductor based high mobility 2DEG, typically such a 



strongly insulating phase occurs ( Jiang et al. 1991 1990 ) 
for i/< 1/5 — 1/7 whereas in graphene the effect man- 
ifest near the charge neutral Dirac point around v 0. 
Whether the same physics controls both insulating phe- 
nomena or not is an open question. 

Very recent experimental observations of = 1/3 
FQHE in graphene has created a great deal of excitement. 
These preliminary experiments involve 2-probe measure- 
ments on suspended graphene samples where no distinc- 
tion between p^x and Pxy can really be made. Further 
advance in the field would necessitate the observation of 
quantized plateaus in Pxy with pxx ~ 0. Since FQHE in- 
volves electron-electron interaction effects, with the non- 
interacting part of the Hamiltonian playing a rather mi- 
nor role, we should not perhaps expect any dramatic dif- 
ference between 2DEG and graphene FQHE since both 
systems manifest the standard 1/r Coulomb repulsion 
between electrons. Two possible quantitative effects dis- 
tinguishing FQHE in graphene and 2DEG, which should 
be studied theoretically and numerically, are the differ- 
ent Coulomb pseudopotentials and Landau level coupling 
in the two systems. Since the stability of various FQH 
states depends crucially on the minute details of Coulomb 
pseudopontentials and inter-LL coupling, it is conceivable 
that graphene may manifest novel FQHE not feasible in 
2D semiconductors. 

We also note that many properties reviewed in this ar- 



ticle should also apply to topological insulators (Hasan 



and Kane 2010), which have only single Dirac cones on 



their surfaces. Although we now have a reasonable the- 
oretical understanding of the broad aspects of transport 
in monolayer graphene, much work remains to be done 
in bilayer and nanoribbons graphene systems. 
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